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INFLUENCE OF THE REACTOR TEMPERATURE CHARACTERISTICS 
UPON THE CHOICE OF THE OPTIMUM THERMODYNAMIC CYCLE 
OF AN ATOMIC-ELECTRICAL GENERATING STATION 


D. D. Kalafati 


Translated from Atomnaya Energiya Vol. 8, No. 1, pp. 5-14, January, 1960 


Original article submitted November 18, 1958 


In this article we investigate the possible temperature changes in a nuclear energy station as 
a function of the thermal power of the reactor when there are two limiting temperatures: for 
the shell and for the center of the heat-emitting elements. We find the changes in the al- 
lowed thermal and electrical power of the unit as a function of the parameters of the thermo- 
dynamic cycle. We give the reader an understanding of the boundary thermal power of the 


reactor and the efficiency of the generator. 


We also give the conditions under which the formulas which we have derived may be 
used for a preliminary calculation of the optimum parameters of the thermodynamic cycle. 
Our analysis gives the curves which show the increase in the parameters and the efficiency 
of the electrical generating station as a function of the material of which the shell is con- 


structed and of the type of nuclear fuel used. 


Changes in the initial parameters of the thermody- 
namic cycle of atomic-electrical generating stations and, 
consequently, of the temperature of the heat carriers leads 
to changes in the thermal power of the reactor as well as 
to changes in the efficiency of the cycle; therefore, the 
optimum parameters of the thermodynamic cycle and of 
the heat catriers are determined only by a simul- 
taneous analysis of the operating conditions of the 
cycle and the reactor. Besides, in order to calculate the 
thermal operation of the reactor we must first know the 
optimum initial parameters of the cycle and of the heat 
carriers. Thisis difficult to do merely by means of varia- 
tional calculations. 

The following formula (1) was obtained on the ba- 
sis of using the condition of minimum cost of electrical 
energy in a preliminary calculation of the optimum 
mean temperature of the heat supply for the heat cycles 
in the steam turbines of the atomic-electrical generating 
stations 


Tr. To¢ 
eK, (1) 


7SBy” oa )/ 1—z 

where “ie is the limiting temperature of the shell ie or 
of the center Th of the heat-emitting elements (HEEL) 
of the reactor; Tracy is the temperature in the condenser; 
Z = mm + Cy (m isthe thermal efficiency of the cycles, cp 
is the thermal component of the cost of the electrical 
energy). 


For a small thermal component where we can as- 
sume that c-p = 0, formula (1) corresponds to the condi- 
tion of maximum electrical power of the generating 
station [2]: 

TyPy™ =V Ti TacyK. (2) 

Formula (2) gives the value of the optimum mean 
temperature not of the heat carriers, as in the formula 
obtained by G. Melets [3, 4], but of the thermodynamic 
cycle. As a result of the completely identical results for 
limiting temperatures of the center and periphery [1, 2], 
in formulas (1) and (2) the temperature 7 must corres - 
pond to the limiting shell temperature or the limiting 
center temperature and which of the temperatures is 
the limiting one ina givenreactor. Since the optimum tem- 
perature of the cycle differs basically, it ismecessary to estab 
lish which optimum is preferable for nuclear energy in- 
stallations; that determined on the basis of the limiting 
temperature of the shell or that determined on the basis 
of the limiting temperature of the center. 

In addition to this we often assume that the optimum 
temperatures exist simultaneously at the periphery and at 
the center, The question then arises as to how we can,in 
such a case, determine the optimum cycle parameters, 


Possible Temperature Characteristics of 
the Nuclear Energy Station 


Let us examine the possible temperature changes in 
a given nuclear energy installation for changes in the 


thermal power of the reactor, due to changes in the ini- 
tial parameters and to the thermodynamic cycle. Let the 
form, the surface, the number of HEEL be given and let 
us choose the velocity and consumption of heat carriers. 
Under these conditions, the thermal power of the reactor 
is proportional to the difference in the limiting tempera - 
tures of the shell (or center) of the HEEL and the average 
temperature of the heat supply pipes in the cycle [1]: 


Qr = k, nk; F,w (Tt yo Tiey) (3) 


where k, is the coefficient of nonuniformity of heat em- 
ission with respect to the reactor radius, n is the number 
of HEEL, ky is the coefficient of heat transmission, Pg is 

the surface of the HEEL, and w is the coefficient of utili- 
zation of the possible reactor power, 


The electric power of the station is equal to 


= Ome PI —AP.y = ee = 


q- ee es 


where Nj Nm Ng 1s the product of three efficiencies: 
the relative inner, the mechanical and the generator ef- 
ficiency; and AP,, is the power required for the opera- 
tion of the station. If we change the thermal power Q, 
for the given constant thermal resistances of the HEEL, 
constant heat transmission to the heat carriers and the 
steam generator, then all the temperature drops in the in- 
dicated elements change proportionally to the thermal 
power of the reactor, However,the temperature changes 
of these elements cannot be the same (Fig. 1) since they 
depend upon the conditions under which the changes oc- 
cur in the thermal power of the reactor, 

For asmall thermal power(Zone I) anincrease in the 
thermal flux leads to an increase in all the temperatures 
of the unit except the condenser temperature Tacy = 
const. In Zone I an increase in the thermal power of 
the reactor leads to an increase in the average tempera - 
ture of the cycle, i.e., in this zone there is no maximum 
electrical power, In the indicated zone we may operate 
reactors that are used for study purposes, but as far as the 
operation of power reactors is concerned, it is always ad- 
vantageous to increase the thermal power until the maxi- 
mum peripheral power T!** up to the limiting value. 

A further Epa reese in the thermal power of the reactor 
must occur for Te = const. This can only be accomp- 
lished by lowering the mean temperature of the heat 
carriers in the cycle,while at the same time increasing 
the maximum temperature of the center HEEL TH14*, 
Since in this process the thermal power increases and 

the cycle efficiency drops, then in Zone II, according to 
equation(4),we get the maximum electrical power of the 
station which corresponds to the optimum mean tempera - 
ture of aS eke quation (2), if we substitute in the for- 
mula a for = Then the point A,corresponding to the 
athena thermat reactor power, is determined by the 


~ APsn (4) 
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Fig. 1, Possible temperature changes in a nuc- 
lear energy installation with TPRy" in Zone I. 


PeQy 
100% 


80 
60 
40 


20 


0 fj 
Tacy 100 = 7 °P™M ——s200 


Fig. 2: Dependence of the thermal and elec- 
trical HOEY of a reactor with TPP. Cy. in Zone 
II upon ange cy" 


point of intersection of the lines of the optimum and al- 
lowed temperature cycles (see Fig. 1). 

The dependence of the thermal and electrical pow- 
er of the station upon the average temperature of the 
cycle for one limiting temperature of the shell or of the 
HEEL of the center is discussed in papers [1, 2]. For two 
limiting temperatures,the allowed changes in the ther- 
mal and electrical power have been computed from equa- 
tions (3) and (4), for the case of an optimum in the av- 
erage power Zone II, as a function of the mean tempera- 
ture of the heat conductor for the cycle Ticy and are 
shown in Fig, 2, We can see from Figs. 1 and 2 that an 


increase in the thermal power in Zone II for T! = const 
is possible until we reach the limiting temperature of 
the center of the HEEL, after which possible changes in 
thermal and electrical power can occur only for Ths 
const, i.e., in the high thermal power Zone III, in this 
case we must substitute an for T! in equations(3) and (4), 

Since the limiting shell temperature is bounded by 
the quality of the materials used or the properties of the 
heat carriers and often does not exceed 300 - 400°C, 
then for a condenser pressure py = 0.04 atm (t, = 28.6°C) 
and tT = 300°C we get from equation (2): 


yop m 


OPM. V573-302 = 416 °K = 143°C, 


this corresponds (for the cited temperature range of the 
shell) to a saturated steam cycle with an initial pressure 
of 5 - 10 atm. Such steam parameters would lead to low 
efficiencies (of the order of 20% and less),therefore the 
operation of atomic-electrical generating stations in 
Zone II would not be economical. For the same thermal 
reactor power an increase in the efficiency of the station 
{the factor ne in equation(4)] increases proportionately, 
or lowers the specific capital losses. For a given electri- 
cal power for the station using nuclear fuel, i.e., 
for the same depth of burning (degree of fuel utilization), 
the duration of the reactor charge is directly proportion- 
al to the efficiency of the station. Therefore an increase 
in the efficiency of atomic generating stations is signi- 
ficant even for the small fuel cost component of elect- 
rical generating stations. 

Since the limiting temperatures of the center of 
the HEEL are significantly higher than those of the shell 
(from 3 = 660°C for metallic uranium to T. = 2800°C 
for uranium dioxide) then the optimum mean tempera- 
ture of the heat conductor during the cycle,as given by 
formula (2) according to the limiting temperature of the 
center of the HEEL, increases discontinuously as for 


200 


Fig. 3, Possible temperature changes in a nuclear en- 
ergy installation with TPR in Zone Il. 


example, for Tl = 1000°C, Tidy = (1275-302 = 623°K 
= 350°C. This mean temperature necessitates the use of 
steam with ultra-high parameters,with station efficiency 
of the order of 40% which is twice as great as the effi- 
ciency of the station during operation in Zone II. In 
connection with this it is expedient to take such possible 
temperature characteristics of the nuclear energy instal - 
lation so that the lines of theoptimum and possible temp- 
erature cycles cross in Zone II and in Zone III; in Zone 
III both the thermal power of the reactor and the opti- 
mum cycle efficiency are greater (Fig. 3). 

Upon reaching the limiting temperature of the 
shell and of the center of the HEEL simultaneously, the 
thermal power, which we shall henceforth call the 
boundary thermal power of the reactor, oe does not as 
yet reach its maximum possible value, as has sometimes 
been indicated in the literature, but lies only at the be- 
ginning of Zone III. A further increase in the thermal 
power of the given reactor in Zone III for Tr = const is 
possible if we decrease the maximum shell temperature 
as well as the mean temperatures of the cycle and of the 
heat carriers. 

Possible changes in the thermal and electrical po- 
wers of the nuclear energy installation in the presence 
of an optimum in Zone III as a function of the changes 
in the mean temperature of the cycle are given in Fig. 
4, As we see from the figure,if we increase the mean 
temperature of the cycle, when operating in Zone III, 
the electrical power of the installation increases at first, 
reaches a maximum Tr Ry, (determ ined from equation 
(2) by substituting Th for T!),and then decreases, After 
we reach the boundary thermal power of the reactor, 
changes in the thermal and electrical power are possible 
only for the conditions of Zone II. 

Thus the economy in the operation of atomic elec- 
trical generating stations operating in Zone III at a high- 
er thermal power and greater efficiency is significantly 
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Fig. 4. Dependence of the thermal and electrical po- 
wer of a reactor with TPRy" in Zone IM upon TP cy: 


greater than that for electrical generating stations oper- 
ating in Zone II where the optimum cycle parameters 
are determined by the limiting temperature of the shell, 
Therefore the best conditions for the operation of atom- 
ic-electrical generating stations are the conditions 
prevailing in Zone III where the optimum cycle para- 
meters are determined by equation(2) according to the 
limiting temperature of the center of the HEEL. Since 
in Zone III the temperatures of the shell and cycle de- 
crease relatively slowly, it is more probable that the 
operation in Zone III will be in the vicinity of point M 
(Fig. 5). 

Due to the independence of the coefficients of con- 
ductivity and thermal emission of the mean temper~ 
ature, the possible temperature characteristics take the 
form of line segments, Fig. 3. If the indicated coeffi- 
cients change somewhat with temperature, then these 
lines will have an insignificant curvature,but the shape 
of the diagram and the division of the characteristics 
into three characteristic zones will not change. We 
wish to underline the fact that the diagram we have 
shown is not a diagram of the changing mode of opera~ 
tion of the reactor during power regulation in the ex- 
ploitation process (for example, at average heat carrier 
temperatures T;" = const, see dotted line in Fig. 1) but 
the characteristic temperature changes in the unit due 
to the choice of thermal power of the reactor, 
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Fig. 5. Possible temperature changes in a nuclear 
eae), installation for a reactor operating with 
Q,. 
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Limiting Thermal Power of the Reactor 


and Greatest Electrical Power of the 
Station . 

Depending upon the relationships between the limit- 
ing temperatures of the shell and the center, the rela - 
tionships between the thermal resistances of the HEEL, 


4 


the processes of heat transfer to the thermal carriers and 
the resistances of the steam generator, an intersection of 
the lines giving the possible changes in the mean tem- 
peratures of the cycle with the lines of optimum mean 
temperature in Zone III is not always possible. In par- 
ticular this is impossible when the limiting temperature 
of the shell of a given reactor in Zone III is lower than 
the optimum mean temperature of the cycle according 
to (2) with respect to the limiting temperature of the 
center, i.¢€., if 


Vey (>) 


The indicated temperature relationship occurs, for 
example,in uranium dioxide water-water reactors. 
Since, for WWER (water-water element reactors) aus = 
2200°C and due to the absence of the boiling away of 
the heat carriers for py = 100 atm,it follows that ql = 
309°C [5]. Then according to the equation (5) 


V 2473-302 = 860 °K > 582 °K, 


i.e., for the indicated type of reactor itis impossible to 
obtain the optimum mean temperature of the cycle in 
Zone III,due to the low limiting temperature of the 
HEEL of the shell, For this case, the mean temperature 
of the cycle must be as near as possible to the limiting 
temperature of the shell, since it is entirely expedient 
to raise the possible temperature characteristics,and the 
extent to which they are raised is limited only by the 
necessary optimum temperature drop in the steam gen- 
erator Atse . The temperature characteristics of such a 
unit are shown in Fig. 5, while the possible changes in 
thermal and electrical power,as a function of the mean 
cycle temperature or of the heat carriers,are shown in 
Fig. 6. For the case under consideration,there is no 
maximum electrical power either in Zone II or in Zone 
III, since the limiting temperatures of the HEEL of the 
shell and center are reached simultaneously for the 
point when, in one of these zones, the maximum is ob- 
tained. However on the boundary of Zones II and III, 
for the limiting thermal power of the reactor, the great- 
est electrical power is obtained at the point of intersec- 
tion of the curves which show the changes in electrical 
power in Zones II and III (Fig. 6). 

The optimum mean temperature of the cycle, cor- 
responding to the value of greatest electrical power, is 
always located between the two values of optimum 
mean temperature of the cycle which are determined 
from equation(2) if we substitute first the value of the 
limiting shell temperature and then the limiting temp- 
erature of the center of the HEEL. The value of the 
mean temperature of the cycle corresponding to the 
maximum electrical power cannot be found from equa- 
tion (2), Therefore,the attempt [6] to find PH* py in- 
troducing, in equation(2), the limiting temperature of 
the shell plus a correction factor cannot be successful. 


However,for the case of the limiting thermal power we 
get from the condition of the simultaneous existence of 
the limiting temperatures of the shell and of the center 
of the HEEL a single-valued determination of the mean 
temperature of the heat carriers. This solves the problem. 
For this case the optimum mean temperature of the 
cycle is determined from the expression 


ine =e = Agem . (6) 

If,for the maximum electrical power there was a 
flat part of the power curve for which it would be pos- 
sible that a small deviation of the mean temperatures of 
the cycle and of the heat carriers from the optimum 
would not result in a drop in the electrical power,then 
the derivative dP, /dTy cy would have a discontinuity. 
At this point E on the electric power curve we geta sharp, 
peak (see Fig. 6) Therefore, even a small deviation of 
the mean temperature of the heat carriers from the op- 
timum can lead to a noticeable lowering of the electri- 
cal power of the installation. For example, the maxi- 
mum temperature of the center of the HEEL, 1200°C, 
for the atomic generating station in Shippingport [7], 
had not yet reached the limiting value for uranium dio- 
xide, i.e., the steam pressure p, = 39.2 atmospheres, 
which was too high for the given installation, corresponds 
to point III on Fig. 6, and does not insure the maximum 
possible electrical power of the station. 

Thus in the case under study, when the limiting 
temperature of the center of the HEEL is reached, it is 
expedient to raise the temperature characteristics for the 
mean cycle temperature, but if for the given installation 
To < To, it is expedient to lower the mean tempera- 
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Fig. 6. Dependence of the thermal and electrical 
power of a reactor with QP upon TY'cy: 


ture of the cycle so that the limiting temperature of the 
center is not reached. 


ici of the Steam Generator 


For the greatest electrical power output of the in- 
stallation there is a boundary value for the thermal power 
of the reactor,but for a constant consumption of heat 
carriers there exists a mean temperature to which the 
heat carriers are raised. For this case,changes in the pa- 
rameters and in the efficiency per cycle of the installa- 
tion can occur only as the result of changes in the temp- 
erature drop within the steam generator due to changes 
in its surface heating or in the coefficient of heat trans- 
mission kp¢. 

If, during alternations in the thermal power of the 
reactor,changes in the steam generator temperature drop 
influence the coefficient of utilization of the thermal 
power of the reactor [1, 2);then for a limiting thermal 
power, changes in the mean temperature drop in the 
steam generator, due to changes in its surface heat, in- 
fluence only the changes in the thermal efficiency of the 
cycle. 

Evaluations made in the literature regarding the 
quality of the steam generator sometimes make use of a 
value of the steam generator which takes into account 
the heat losses. These can usually be neglected. 

For a given reactor,thermal power ope const due to 
the fact that the heat exchange process in the steam gen- 
erator is not reversible, the index of thermodynamic per- 
fection of the steam generator must be set up not in 
terms of the heat but in terms of the power efficiency of 
the steam generator and be equal to the ratio of the max- 
imum work capability of the heat,before and after the 
heat exchange in the steam generator: 


AL uarey ang TO ie ra (7) 
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SALT (Th — Trey) Troy 


where AL®Y represents the operation of an ideal steam 
cycle (Fig. 7) and AL{°Y represents the maximum pos- 
sible work capability of the heat received by the heat 
carrier, which is equivalent to the operation of the re- 
versible cycle,abcda,with an alternating temperature in 
the heat carrier heat supply. The work ratio in equation 
(7) is equal to the ratio of the efficiencies of the steam 
and reversible cycles expressed in terms of the tempera - 
ture,since they represent the same amount of supplied 
heat, 

In order to evaluate the thermodynamic perfection 
of the steam generators in atomic-electrical generating 
stations, we wish to point out, for example, that the pow- 
er efficiency of the steam generator for a station with 
(WWER) computed from equation (7) is fgg = 0.92,that 
for the atomic icebreaker “Lenin” is se =. 0.835, i.e., 
the power efficiency of the steam generators is close to 
the thermal efficiency of the boilers, The optimum 


Fig. 7- Loss in the working capability of 
the heat during nonreversible heat exchange 
in the steam generator. 


temperature drop in a steam generator is uniquely de- 
termined by comparing the surface heating of the steam 
generator and the energy losses during heat exchange, 
however this requires a special derivation. 


Electrical Energy Upon the Optimum 
Cycle Parameters | 

Previously we have, for simplicity,examined in- 
stallations in which the cost of the fuel components of 
the electrical energy was so small that it could be neg- 
lected, i.e., for which the optimum parameters were de- 
termined by the condition of maximum electrical pow- 
er. 

When the cost of the fuel component reaches 0.1 
- 0.2 and higher,the optimum mean temperature of the 
cycle must be determined from the condition of the 
minimum electrical energy cost, Changes in the cost 
of the electrical energy may be expressed as a function 
of the electrical power and the efficiency of the station 


[1]: 


C= Cx ai Cr = a ai 7 cents/kw hour, (8) 
where K denotes the capital outlay,"amortization, sal- 
aries, and other exploitation expenses which are propor- 
tional to time and computed per hour of operation of 
the atomic power generating station,and V is the cost of 
nuclear fuel used per kilowatt hour of heat generated. 
The optimum mean temperature of the cycle, tak- 
ing the cost of the nuclear fuel into account for K = 
const ,is determined from equation(1). This preserves all 
the previously obtainedconclusions regarding the influ- 
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ence of the possible temperature characteristics of the 
installation upon the optimam mean temperature of the 
cycle. 

The optimum mean temperature and efficiency of 
the cycle corresponding to a minimum cost of electri- 
cal energy for either Zones II or III (determined res- 
pectively from the limiting temperature of the shell and 
the central HEEL) will be somewhat higher than the 
temperature and the efficiency corresponding to maxi- 
mum electrical power (2), 

However, if the unit operates at the limiting ther- 
mal reactor power, then as a result of the discontinuity 
in the power derivative dP./dTf"cy at the point E (see 
Fig. 6) then the derivative of the costs dC, /dT}'cy will 
also have a discontinuity at this point. Instead of ob- 
taining the minimum cost for the limiting thermal reac- 
tor power, we will get the minimum value of the cost of 
the electrical energy (except for the boundary portions 
at points M and P, Fig. 5). Thus,for the limiting ther- 
mal power of the unit the maximum electrical power 
usually coincides with the condition of minimum cost 
of electrical energy. 

For uranium, for Tl= 650°, cp = 0.40 and ny = 
0.45, the coefficient z is0.18. Under these conditions, 
the optimum mean temperature of the heat supply is 
equal to 

opm = Yi agg = 585° K = 312°C. 

At the same time for the cycle we have used for py 
= 90 atm, ty = 500°C and t = 215°C (temperature to 
which the water is heated), the mean temperature of the 
heat supply is given by: 

m iy—in 810.1—220,4 _ 
Ti cy 4 poe og 004s 0a SaGae 

It follows from the example cited, that the optimum 
initial cycle parameters for atomic-electric generating 
stations,using metallic uranium with ee = 650°C as de- 
termined from the conditions for the minimum electri- 
cal energy cost,consist,for example,of 90 atm pressure 
and 500°C, For the type of fuel considered a further in- 
crease in the initial steam parameters would be inad- 
visable. This conclusion coincides with the data of the 
project for the possible improvement in the cycle of an 
atomic-electric generating station using uranium with a 
sodium graphite reactor SGR [8]; in this project an in- 
crease in the initial steam parameters from 56 atm 
(440°C) to 88.5 atm (510°C) was foreseen. In this pro- 
cess the electrical power of the station does not increase, 
and an increase in the steam parameters results only in 
a decrease in the cost of the electrical energy. A fur- 
ther increase in the cycle parameters was advisable only 
for an increase in the temperature of the central HEEL, 
i,e., in the case where another nuclear fuel or an alloy 


585° K=312°C. 


with a greater limiting temperature was used. 


“The amortization time of an electrical generating 
station is taken as eight years, 


The Line of Optimum Temperatures for 
the Heat Supply During the Operating 
Cycle of an Atomic-Electric Generating 
Station 

All the positions which we have examined for the 
line giving the optimum mean cycle temperature in the 
various zones,can be connected with one line of opti- 
mum mean cycle temperatures in TQ, coordinates, 
this is the heavy line in Fig. 8. 

The point of intersection of the line of possible 
temperature characteristics of the cycle, which can be 
drawn during the planning of the unit, with the line of 
optimum cycle temperatures will correspond to the op- 
timum operating conditions of the given atomic-elec- 
tric generating station. For the case of operation at the 
limiting thermal reactor power, the point of intersection 
of the possible changes in the mean temperature of the 
heat carriers with the vertical portion of the line giving 
the optimum conditions corresponds.also to the mean 
optimum temperature of the heat carriers. 

If we take into account the component giving the 
fuel costs of the electrical generating stations,the line 
showing the optimum cycle temperature in Zones II and 
Ill is raised,depending upon the coefficients z = c-pn, in 
formula (1); this corresponds to the dotted lines in Fig. 
8. We wish to underline however,that all the points on 
the optimum temperature line of the cycle are of equal 
weight. The best conditions, which insure simultaneous - 
ly the maximum thermal power and efficiency of the 
unit,lie in Zone III near point M. 

The optimum mean temperatures of the cycle or of 
the heat carriers,for a given nuclear energy installation 
having a given form of, number of and type of surface 
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Fig. 8- Line of optimum temperatures for ey” 
of an atomic-electric generating station. 


HEEL, and with a given velocity and utilization of heat 
carriers and with a constant surface heating of the steam 
generator, were considered previously. 

If the characteristics of the possible temperature 
changes in a nuclear energy installation do not insure 
an optimum mean cycle temperature in Zone III (de- 
termined by equation (2) from the limiting center tem- 
perature), then we should raise the line of the possible 
temperature characteristics of the heat carriers in the 
cycle by changing the ratios of the thermal resistances 
of the HEEL, the heat transmission to the heat carriers , 
and the steam generator resistance,in order to insure the 
choice of the optimum parameters. Such an increase in 
the temperature characteristics of the heat carriers and 
of the cycle is possible if we change the diameter and 
heat conduction of the HEEL, the surface temperature of 
the steam generator,and the coefficient of heat trans- 
mission, depending upon the speed and parameters of 
the heat carriers and the working substance. It is pos- 
sible to entirely eliminate Zone II and obtain only Zone 
III and its left boundary. 

The method of obtaining more favorable thermody- 
namic operating conditions for the nuclear energy in- 
stallation consists in raising the possible temperature 
characteristics, The expediency of such an increase in 
the indicated characteristics must be verified by cal- 
culations, 

In order to evaluate the applicability of the pre- 
viously derived conclusions, we show in Fig. 9 curves of 
the optimum mean temperatures of the heat supply dur- 
ing the cycle computed fromequation (1) as a function 
of the limiting temperature of the shell and of the cen- 
tral HEEL and the fuel cost component of the electrical 
energy. Points 1 - 17, Fig. 9, correspond to given 
mean temperatures of the thermodynamic cycles of the 
constructed or proposed electric generating stations 
listed in the table. 

The computations for points 1 - 4 (Fig. 9), which 
we have carried out either as examples or to be used for 
double -purpose stations agree well with the theoretical 
curves computed using the limiting shell temperature, 
and points 5 - 11 characterize stations using uranium as 
a fuel,computed using the optimum cycle temperature 
for the limiting temperature of the central HEEL and for 
a fuel cost component of electrical energy from 0 — 
30%, 

When we considered the case of uranium dioxide 
as the nuclear fuel, we assumed that the temperature of 
the center was 2200 - 2760°C,and that the optimum 
mean temperature of the operating cycle of the atomic- 
electrical generating station (points 12 - 17, Fig. 9) 
was higher than the limiting temperature of the shell. 
Therefore such stations operate at the limiting thermal 
reactor power. 

For a limiting temperature of the central HEEL, of 
the order of 1000 - 1200°C, the optimum parameters 


Mean Temperature of the Heat Supply over the Heat Cycles in the Installation 


: Point No. 1 o Tia e 
Condition (see Fig.9) | Name of Reactor or Station | toh | Trey» © 
l * po ‘ : a 
nie graphite reactor (calculations in) 930 118 
For a limiting tempera 2 
ture of the shell HEEL Uranium- graphite reactor (calculations in 300 137 
3 [4]) i 
r Two-purpose reactor at Marcoule (France) ae We 
Reactor at Calder Hall (England) | 
5 First atomic generating station (USSR) 370 485% 
6 Station at Henderstone (England) 569 227 
7 | Heavy water reactor with organic heat ca 
riers (switzerland) 600 224 
§ ' Uranium- graphite reactor using super- 
For a limiting temper heated steam (USSR) 550 312 
erect heicentral 9 Sodium- graphite reactor at Santa Susana sigs sigs 
HEEL 10 Reactor SGR (USA), variation 1[8] 650 269 
: H Reactor SGR (USA), variation 3[8] 650 297 
Reactor “Enrico Fermi" (USA), variation 1 713 262 
dls 
| Reactor “Enrico Fermi” (USA), variation 2 713 287 
3 1 1 
12 eee station at seh viet suas) ae we 
Ae 13 oiling reactor at Kale (FRG) 176 247 
For a limiting thermal 14 |WWER reactor (USSR) [5] 2200 235 
power of the reactor 15 Dresden reactor (US A) 2520 298 
16 “yankee Atomic” station (USA) 2600 243 
17 Reactor SGR (USA) 2760 297 


"For points 5, 8, 10, 12, and 14 in paper [6] the values given for Tey are not correct: 40-60°C lower 
without taking regeneration into account 
** Point 6, Fig. 3, paper [1] for the station in Shippingport is erroneous. 
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Fig. 9. Comparison of age for stations in operation, being planned or being built,with the theoritical 
curves for TPRy" for various fuel components. 


will become so high that they will be unattainable for ample ,the critical thermal loading, the limiting temp- 

cycles utilizing ‘ultrahigh parameter steam. Then the erature of the moderator, the radiational instability, 

choice of the cycles for steam turbine installations may ete., which affect the temperature characteristics 

be based on the parameters attainable for steam and the which we have presented in TQ, coordinates and which 

determination of the optimum parameters will be nece- require a special analysis for each case. 

ssary only for gas-turbine installations. On the basis of the data presented and also of the 
We have so far considered the temperature limita - typical computations shown in papers [1, 2], we can con- 

tions of HEEL. There are also other limitations, for ex- clude that equations (1) and (2) may be en to obtain a 


preliminary estimate of the optimum parameters for the 
thermodynamic cycle of the steam turbine installations 
of atomic electric generating stations. A theoretical 
analysis using the indicated formulas enables us to find 
the factors which influence the optimum cycle para- 
meters and also shows us how to develop and improve 
the thermodynamic operating cycles of atomic-electric 
generating stations. 
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A large detector filled with a liquid organic scintillator containing cadmium has been used to 


. . Cael . 235 ce. 
measure the number of neutrons emitted by individual fragments inU™ fission 


by thermal 


neutrons, Using 4m geometry the dependence of the umber of neutrons emitted by Uae 
pairs on the mass ratio has been measured, The excitation energy used in neutron Pies. 
is determined on the basis of the semiempirical Weizsacker formula. A sharp asymmetry in the 
distribution of excitation energy between the heavy and light fragments is noted. The data 
which are obtained are found to be in disagreement with the statistical theory of fission pro- 


posed by Fong. 


Introduction 

The distribution of excitation energy between frag- 
ments is an important characteristic of the fission pro- 
cess in heavy nuclei. In fission by thermal neutrons, the 
initial excitation energy of the nucleus is equal to the 
binding energy of the the neutron in the intermediate 
nucleus which is formed and is approximately 6 Mev, 

a figure which is many times smaller than the mean 
value of the excitation energy of the fragment pair 
which, in fission of uranium isotopes, is approximately 
30 Mev. Considerations based on a quasi-static process 
[1] indicate that the fragments remain in an unexcited 
state up to the moment directly preceding the division 
of the neck which connects them. 

Bohr and Wheeler [2] have indicated that the exci- 
tation of the fragments is basically connected with the 
fact that the shape of the fragments directly after for- 
mation does not correspond to equilibrium; the energy 
of deformation is consequently transformed into energy 
of internal degrees of freedom. B. T. Geilikman [3] 
has considered this question in detail, relating the ex- 
citation energy of each fragment with the deformation 
parameters and the shape of the nucleus at the instant 
preceding the actual division into fragments. The ex- 
citation energy E, is, as is well known, dissipated in the 
emission of “prompt” neutrons and y rays: 


E (M), = v(M)e(M) + ey (M), 


where M is the mass of the fragment; D is the mean 
number of neutrons emitted by the fragments and 

is the energy carried away by the y photons; e(M) is 
the mean energy required for the evaporation of a sin- 
gle neutron, The energy carried away by the y pho- 
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tons is a very weak function of the ratio of fragment 
masses [4] and their excitation energies [5]. On the 
other hand, direct experiments show [6] that an increase 
in the excitation energy of the nucleus leads to an in- 
crease in the number of neutrons emitted in fission and 
that on the average the emission of a single neutron re- 
quires an energy of approximately 7 Mev. Hence,an 
investigation of the neutron yield from the separate 
fragments makes it possible to obtain data on the dis- 
tribution of excitation energy from the individual frag - 
ments and thereby, information on the deformation of 
the fragments at the time of formation. The depen- 
dence of excitation energy of a fragment pair on mass 
ratio is of interest from the point of view of checking 
the statistical theory of fission proposed by Fong [7]; 
this theory predicts an increase in the excitation energy 
of the fragment pair characterized by the most probable 
mass ratio, 

The neutron yield from individual fragments in 
U™ fission by thermal neutrons was first measured by 
Fraser and Milton [8] in 1954, The results obtained by 
these authors indicate that in fission,which is almost 
symmetrical in mass,the excitation energy is distributed 
between the fragments in a manner which is far from 
symmetrical; for a fragment mass ratio of 1:1 the light 
fragment emits approximately four times more neutrons 
that the heavy fragment, The total neutron yield is a 
weak function of the mass ratio, In the region of the 
most probable mass ratio there is apparently a small in- 
crease in the number of emitted neutrons. This last 
feature has been interpreted by Fong [7] as evidence of 
the validity of his theory, From the point of view of 
fission physics it is interesting to determine whether the 
pattern observed by Fraser and Milton is an accidental 


one, and is characteristic only of the fissile nucleus 

u™4 or whether the excitation energy distribution, which 
has been observed,is a general feature which pertains 
under certain conditions for all nuclei (similarly, for 
the fission asymmetry). 

In the present work we have measured the average 
number of neutrons emitted by the individual fragments 
and the fragment pairs, in fissionof U® by thermal neu- 
trons. 


Experimental Apparatus 


The measurement of the neutron yieldv(M) from 
individual fragments is based on the fact that the an- 
gular distribution of the neutrons with respect to the di- 
rection of motion of the fragment is highly anisotropic: 
approximately seven times more neutrons are emitted 
in the direction of motion than at an angle of 90°, 
while the emission of neutrons in the backward direc- 
tion is essentially negligible [9]. A schematic arrange- 
ment of the apparatus is shown in Fig. 1. A layer of 
u", 30u g/cm? thickis deposited by evaporation on a 
colloidal film 20 g/cm? thick which is placed on the 
central electrode of a double grid ionization chamber. 
The chamber is filled with a mixture of carbon dioxide 
and argon (partial pressures of 50 and 900 mm Hg res- 
pectively). The angle of emission of the fragments 
with respect to the normal to the layer is defined by 
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Fig. 1. Diagram of the apparatus: 1) double grid ioni- 
zation chamber; 2) scintillation tank for neutron de- 
tection; 3) photomultipliers; 4) collimated neutron 
beam, 


the collimator; this is a plate of phosphor bronze 0.4 
mm thick in which an aperture 0.4 mm in diameter has 
been drilled. The most probable emission angle for the 
fragments is approximately 25°. The uranium layer is 
located at a distance from the neutron detector such 
that the most probable angle at which neutrons are de- 
tected is also approximately 25°, 

The neutron detector is a hexagonal tank with a 
volume of 200 liters which is filled with a liquid scin- 
tillator, The scintillator is a solution of 2,5-diphenyl- 
oxazole (PPO) in dioxane witha concentration of 4 g/ 
liter which contains an almost saturated water solution 
of cadmium nitrate in amounts such that there is one 
cadmium atom for each 400 hydrogen atoms. The scin- 
tillation flashes ,produced by the captured y photons, are 
recorded by 32 FEU-24 photomultipliers with cathodes 
78 mm in diameter, At the center of the detector there 
is an aperture 113 mm in diameter through which the 
collimated neutron beam passes. The lifetime of the 
neutrons in the scintillator is llpsec. The measured ef- 
ficiency for detection of neutrons from U™® fission in the 
detector is 62%, The efficiency is determined by locat- 
ing the layer of fissile material at the center of the de- 
tector. With the exception of the scintillator,the detec- 
tor is essentially the same as that described in detail by 
Reines, Cowan, et al. [10]. 

The operation of the electronic system is controlled 
by the coincidence pulse from the fragments. This 
pulse opens for 25 psec a gate which transmits pulses 
from the detector to a high-speed counting system and 
pulses from the double ionization chamber to the input 
of the ratio analyzer, The gate is opened after a delay 
of approximately 0.6 sec from the fission time in or- 
der to avoid detection of the prompt y rays, The re- 
solution time of the neutron counting channel is approx- 
imately 0.4 usec. Information for each event (the mag- 
nitude of the ratio and the number of neutrons) is recor- 
ded on a recording system. The ratio range from 2,2 
to 1/ 2,2 is covered by 30 channels of the analyzer. 
Cases in which the light fragment is emitted in the di- 
rection of the detector for ratios from 1/2.2 to 1 are an- 
alyzed by the first 15 channels while cases in which the 
heavy fragment is emitted in the direction of the detec- 
tor for ratios from 1 to 2.2 are recorded by the remaining 
15 channels, Under operating conditions approximately 
80 fissions/min are recorded, 

We measured 500,000 fission events in the course 
of these measurements; this corresponds to approximate - 
ly 215,000 neutron pulses of which 85,000 were caused 
by fission neutrons,with the remainder due to the back- 
ground of scattered y rays and neutrons from the beam. 

The data on the total yield of neutrons from frag - 
ment pairs can be obtained by measuring the neutron 
yield from the individual fragments v(M). However, in 
order to obtain the quantity v(M), it is necessary to in- 
troduce a correction in the experimental results; this 
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correction takes account of the dependence of neutron 
detection efficiency on angular distribution with respect 
to the fragments, To determine the total yield indepen- 
dent measurements are made in which the ionization 
chamber containing the layer of fissile material is 
placed at the center of the detector and the neutrons are 
detected in a 4m geometry. The apparatus is operated 
in the same way as in the measurements of neutron yield 
from the individual fragments with the only difference 
being that the ratios x and1/x, corresponding to the 
same mass ratio for the heavy and light fragments,are 
entered in the same channel of the analyzer. In these 
measurements the counting rate for coincidences be- 
tween fragments is15/min while the background is equal 
to one count per fission, In all approximately 70,000 
fission events were recorded, 


Corrections for Absorption in the Layer 
and lonization Defect 

In the ionization chamber used to measure the neu- 
tron yield from individual fragments, one of the frag - 
ments passes through the colloidal film and the colli- 
mator whereas the other loses energy only in the uran- 
ium layer. The average energy loss in the uranium lay- 
er is less than 0.4% and no corrections for this effect 
have been made. The energy loss in the film and the 
collimator are determined from the displacement of the 
point corresponding to symmetric fission with respect to 
the analyzer channel which records equal-amplitude 
pulses. The energy losses of the fragment in symmetric 
fission are found to be approximately 5 Mev. The ap- 
propriate corrections are introduced in the experimen- 
tal distribution of fragment mass ratios, In introducing 
the corrections account has been taken of the depend - 
ence of initial ionization of the fragments on velocity 

In the measurement of the total neutron yield,one 
of the fragments passes through the collimator while the 
other passes through the film. The energy losses of the 


fragments are approximately the same in magnitude and 
the displacement of the mass-ratio distribution curve is 
less than 1%. Hence,we have not deemed it necessary 
to introduce corrections for absorption,and record events 
characterized by the ratio x and 1/x in the same chan- 
nel, 

The yield curve for the fragments is also corrected 
for the ionization defect, In introducing this correction 
we use the values 5.7 and 6,7 Mev [11] for the ioniza- 
tion defect of the most probable light and heavy frag - 
ments respectively. It is assumed that the defect is a 
linear function of fragment mass; AE(M) = (4 + 0.019M) 
Mev [8]. 


Correction for Neutron Detection 
Efficiency ai Hed 

This correction must be introduced only in the 
measurements of neutron yield from the individual frag - 
ments. The larger the fragment velocity, the greater the 
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fraction of the neutrons emitted by it which enter the 
detector and are recorded, 

The angular distribution of isotropically emitted 
neutrons from a fragment converted to the laboratory 
coordinate system has been given, for example, in [9]. 
In introducing corrections, we assume that the neutron 
detection efficiency n(M,q) for neutrons which are emit- 
ted in the coordinate system of a fragment with energy 


q is proportional to the probability for emission of a 


neutron at an angle 0° in the laboratory coordinate sys- 
tem. Under these conditions 


1 (M, q)=const (1+7)?, raV 4) iis > 4) 


where E(M) is the kinetic energy of a fragment with 
mass M, For the most probable energies q and E the 
value of r is slightly smaller than unity. In going from 
symmetric to the most probable fission the efficiency 

n varies by approximately 20%, It can easily be shown 
that the magnitude of this variation is a very weak func- 
tion of the energy q in the range 0.5 - 2,0 Mev, which 
covers most of the neutrons. Hence, although the neu- 


. tron spectrum is continuous,it may be assumed that the 


detector efficiency is determined by the mean neutron 
energy q(M) emitted by a fragment of a given mass, 
The mean value of q(M) for the most probable mass de- 
pends on : the excitation energy and is given by 
0.62 0+1)2 Mev, where V is the mean number of neu- 
trons .which appear per fission event (Terrel [12]). For 
as qis 1.15 Mev. Since,experimentally,one observes 
a strong dependence of fission excitation energy on 
mass, the corresponding change in mean neutron energy 
is taken into account. It is assumed that 


q(M) = 1.5 + [ vimy—s | k. 


The factor k can be found either from the relation 

q = 0.62(v +1)? Mev or from the experimental results 
by studying the change in fragment temperature as a 
function of excitation energy [13]. The quantities 
found by these two methods differ somewhat and the 
average value is taken: k = 0,25 Mev/neutron, Thus 


q(M) =0.84-+- v (M) 0.25. 


Substituting the value of q(M) in Eq. (1) and expressing 
the fragment energy in terms of its mass and the total 
kinetic energy of the fragment pair Ep (x) with respect 
to the mass ratio x, finally we have 

const > 


n(M)= 
Ei ( a) 


x M)—M 2 
x [44 V5 V momen |: i 


/ 


where Mg is the mass number of the fissile nucleus. 
The value of the constant is found from the normaliza- 
tion condition: the results must be normalized to the 
known average value of v for U* which is taken to be 
2.45, 


Discussion of Results 

The results of the measurements are shown in Figs. 
2and 3, The mass ratio of the fragments is plotted 
along the abcissa axis, The fragment yield per unit 
mass ratio, corrected for the ionization defect and ab- 
sorption, is shown in Figs, 2a and 3a. For purposes of 
comparison we also show the results of radiochemical 
and mass-spectrometer studies of fission products [14] 
(Fig. 3a, dashed curve). 

In Fig. 2b is shown the neutron yield for individual 
fragments corrected for the efficiency of the neutron 
detector. Analysis shows that on the average the light 
fragment emits about 17% more neutrons than the 
heavy fragment. In Fig. 3b the points indicate the 
measured neutron yield from the fragment pair while the 
crosses denote values of vy (x) computed from the results 
of the measurements of v from individual fragments. A 
comparison of the calculated and measured values shows 
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that although the correction for detector efficiency is 
based on a large number of simplifying assumptions, the 
error is less than approximately 5% when these correc 
tions are taken into account, 

The following remarks should be made with regard 
to these results. The values of v obtained experiment- 
ally are averaged over some range of ratios, This aver- 
aging takes place primarily because of the finite resolv- 
ing power of the ionization chamber. The resolving 
power Ax/ x (Ax is the dispersion in the value of X) can 
be determined from a comparison of the obtained frag- 
ment yield curve with the accurate data obtained by 
radiochemical and mass-spectrometer analysis, An an» 
alysis such as that carried out by Leachman [15] shows 
that in the present experiment Ax/x = 0,08, It is as- 
sumed that for all values of x the error distribution is 
given in a Gaussian form with dispersion A x which depends 
linearly on the ratio: A x= const 

The effect of the finite resolving power is specially 
noticeable where the fragment yield falls off sharply. 
Hence, in the region close to symmetric fission ( x<1. 20) 
the behavior of the quantity v(M) is not well defined. 

In Figs, 2c and 3c we also show the calculated val- 
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Experimental data for the individual fragments: a) yield distribution for fission 


fragments; b) neutron yield; c) energy carried away by the neutrons. 


13 


19 


~ 
lo ~) 


€(x} v(x), Mev 
Aas 


v (x) 


mR NNRASN MANNS 
OS =AwWKAMNYDW™SO 


fragment yield 


oof] Lae ee 
LOI ICR IMIS O67 CMS 20 me 


mass ratio 


Fig. 3. Experimental data for fragment pairs: a) yield 
distribution of fission fragments; b) neutron yield; 
c) energy carried away by the neutrons, 


ues of v(M) €(M) and v(x) €(x),the energy carried 
away by the neutrons from the individual fragments and 
the fragment pair. Stein and Whetstone [16] have attemp- 
ted to determine experimentally the mean value of the 
energy €(x) dissipated in evaporation of a single neu- 


tron from a fragment pair. However, because of effects 
associated with the finite resolving power of the appara - 
tus with which the velocity of the fragments was deter- 
mined the value obtained by these authors is approxi- 
mately 14 Mev/neutron, which is twice as large as the 
known value of approximately 7,0 Mev/neutron, The 
results contain a correction which amounts to approxi- 
mately 100% of the quantity being measured. In the 
present work the resolving power is smaller than that of 
[16]eand we find it more convenient to compute the 
value of €(M). 

The quantity €(M) is equal to the sum of the neu- 
tron binding energy and its mean kinetic energy q(M) 
in the fragment system, The Weizsdcker semiempirical 
formula can be used to compute the neutron binding 
energy. However, to carry out this procedure one must 
know the charge distribution between the fragments. 
The charge Z) (M,) is found from the relation 


z= Darah + AZ(M1), 
AZ(M,)=5.9—0.09M , (3) 
where the subscripts "1* and "h" denote the light and 
heavy fragment,respectively. The coefficients in the 
expression for AZ are chosen to give equal lengths for 

the B-decay chains and to make AZ vanish for sym- 
metric fission, unity for the most probable fission and 
equal to two for a mass ratio of two, Equation (3) gives de- 
tailed values of Z. For each mass the value of the bind- 
ing energy is smoothed for the odd-even effect in the 
number of protons and interpolated to the exact value. 
The dependence of binding energy on fragment mass 
obtained in this way is shown in Fig, 4. The values of 
the binding energy for integral values of Z are taken 
from the table given by Cameron [17]. 
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Fig. 4, Computed values of the neutron binding energy. 
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Fig. 5. Neutron yield from individual fragments vj /vp and total yield from the fragment pair v (X)/v 
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The values of € (x) are given by the weighted av- 
erage values: 


6 (M)v(M)+e(Mo—M)v(Mo—™M) . 


€ (zt) = v (M)--v (My—M) 
Mo—M _,, 
et Nd reaeh ie 


Several months before the present work was finished 
data were published in [16] on the neutron yield from 
fragment pairs in spontaneous fission of Cite At the 
Second International Conference on the Peaceful Uses 
of Atomic Energy, Leachman [12] reported on the re- 
sults of measurements of neutron yield from Cf® frag- 
ments obtained by Whetstone. 

In Fig. 5 we show the results for iow [8], cf® and 
u?% (data from the present work); the neutron yield 
from the individual fragments is converted into the 
ratio ¥, /v;, while the total yield is given in units of 
v(yx/D The arrows on the abcissa axis indicate the 
‘most probable values for the fragment mass ratios, 

As is apparent from the figure, the general charac- 
teristic for all fissile nuclei is the extremely sharp de- 
pendence of the number of emitted neutrons on frag - 
ment mass, For mass ratios of 1.20-1.25 practically all 
the neutrons are evaporated from the light fragment, 
while for strong asymmetry the neutrons are practically 
all emitted from the heavy fragment. At the same 
time the total number of neutrons emitted by the frag - 
ment pair is a weak function of fragment mass ratio, 

In the region of ratios which covers approximately 90% 
of all the fission cases,the value of v (x) remains con- 
stant to + 10%, For U”® and Cf?™ no increase in v(x) is 
found in the region of the most probable ratio; an in- 


crease of this’kind is predicted by the Fong theory. The 
results for U8? are obtained by combining the yields 
from the individual fragments and it is possible that the 
small increase in the yield is a consequence of a combi- 
nation of experimental errors, As is apparent from Figs. 
2, 3, and 4 the excitation energy of the fragments is 
approximately proportional to neutron yield. While the 
excitation energy of the fragments is basically asso- 
ciated with deformation, the nature of the excitation 
energy distribution between the fragments indicates that 
in fission which is almost symmetrical (in mass),the 
light fragment is deformed much more strongly than the 
heavy fragment. For the most probable mass ratio the 
deformation of the fragments is comparable while for 
high fission asymmetry the deformation of the heavy 
fragment is somewhat greater, It is possible that this 
pattern is connected with the oscillatory nature of the 
rupture of the neck of the fissioning nucleus." 

The authors wish to express their gratitude to B. G. 
Erozolimskii, who was one of the people responsible for 
initiating this work; he also helped in its execution by 
a number of valuable comments and discussions, The 
authors also wish to thank A. A, Markov and A, A, 
Borodin for the design and construction of the multi- 
channel ratio analyzer and K, S. Mikhailov for advice 
and help in the preparation of the scintillator for the 
neutron detector, 


* The suggestion as to the oscillatory nature of the rup- 
ture of the neck of the fissioning nucleus was proposed 
by A. Bohr at the end of 1958 (private communication 
to V. M. Strutinskii). 
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In this paper we consider an approximation method for determining the critical parameters of 
a body of arbitrary shape (or a group of bodies) made from fissionable material. In contrast 

with the variational method, the method described here is extremely simple and gives an er- 
ror which is favorable from the point of view of safety as applied to fissionable materials; 


hence it is useful for engineering calculations. 


The method has been applied to fast neutron systems where the single-velocity approxi- 


mation gives good accuracy. 


Derivation of the Basic Equation 

The critical parameters of a body made from a 
material which is fissionable under neutron bombard- 
ment are described by the Peierls integral equation [1] 


e@= layer) K@,r)av’, a 


where ¢ (f) is an eigenfunction which coincides with the 
neutron density distribution function in the critical state; 
a(f) is the reciprocal mean free path length for the 
neutrons; y is the eigenvalue of the equation. 

The kernel K(r, r') is defined by the expression 


exp (— \ adl) 


Ker re inl? : 


l=|r—r’'|. 
The integration is taken over the entire volume V. 

Equation (1) is valid for a body or a material of 
uniform properties; the density need not be constant 
over the volume of the body. 

In the critical state the first eigenvalue y coin- 
cides with the material constant, which is equal to the 
reciprocal mean neutron multiplication factor for one 
collision with a nucleus: 


OF +95 + Ge 


In VO;+Os Y 


a 


where Of, Os and dg are the elementary cross sections 
for fission, scattering and capture; v is the number of 
neutrons per fission. The difference between the eigen- 
value and the material constant Ay = ym ~- y charac- 
terizes the departure of the system from the critical 


state. The excess value of the quantity being sought.y. 
is favorable from the point of view of safety for applica- 
tion to fissionable materials because in this case the ac- 
tual subcriticality of the system Ay is too low. 

An upper estimate for the eigenvalue y can be ob- 
tained from Eq. (1) if we know the local maximum for 
the neutron density rp. In this case, replacing ¢ (r) by 
the constant value we find the following expression for 
the upper limit of y 


ye \ a(r)K (ro, r) dV. 


However, in practice this method is inconvenient 
because of the fact that it gives a large error and be- 
cause the location of the neutron density maximum in 
a body of arbitrary shape is not, as a rule, known pre- 
cisely. Generally we can only determine the region in 
which this maximum is located, 

In order to obtain a more feasible method of mak- 
ing an estimate we consider region Vj, which includes 
the point which contains the neutron density maximum. 
For brevity, in what follows we will call the region Vj 
the primary inner region while the remaining volume of 
the body will be called the supplementary external re- 
gion; we use the subscripts i and e to denote quantities 
referring to these regions, 

We assume that y(r)= const, and integrate Eq. (1) 
over the volume of the primary region: 


4 , , h? 
v=> NC pial A 


We transform the inner integral, We make use of 
the fact that the region of integration in it can be ex- 


it 


tended over the entire infinite space because outside 
the body a = 0 (this factor is important for extension of 
the results to a system of interacting bodies), The in- 
tegration is carried out in two steps: first we integrate 
over the direction Q along a ray which originates at the 
point fj, and then over the total solid angle: 


(a (e')K (yr) av’ = 


exp (— f adt) & 
“oer yraty wee daQ = 


U-»g 
eS \ ca \ e “du, 
4n 
4am 0 

where u= 


U gp is the optical 


agl us = eet 


Sr 


path length along the line from the point rj to infinity. : 
Integrating along the ray we obtain 


ees —u., @Q 
we 
We now change the order of integration and divide 
the optical length of the ray ug into two parts: 


Woo = Ln} t, 


where L is the optical distance along the ray from the 
point rj to the surface of the primary region; t is the op- 
tical path length along the same ray from the surface of 
the primary region to infinity. Then we have 


ape *e 
Y= V, \ aa 


( e-Le-' dV. (2) 


i 

We transform the inner integral which is a function 
of the direction 9. We consider the case in which the 
density of matter is constant in the primary region. Ex- 
pressing all linear dimensions in the primary region in 
free-path lengths,we set dV = dLdS, where dS is an ele- 
ment of area normal to the direction Q. We integrate 
along dL: 


where Lo is the “optical thickness” of the primary region 
along a ray which passes through the elementary area 

dS in the direction of Q(Fig. 1). We replace the optical 
length t on the average value over the cross section S, 
by the value t (Q): 


1 = ¢ . 
t—y= 7 \ eWo).c8 \ (1-e-to) ds. (8) 
Ss 
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The relation in Eq. (3) becomes simpler if we take 
the primary region to bea sphere. In this case, because 
of the spherical symmetry, the integral over S in Eq. (3) 
becomes independent of direction. 

For an isolated primary sphere the critical parame~ 
ter y is determined in our approximation by Eq. (3) with 
t = 0: 


1 BL 
1—yo=7-\ (t-¢ o) dS. (4) 
Ss 
Dividing Eq. (3) by Eq. (4),term by term, we obtain 
the relation between the quantity being sought,y, and the 
value of yo for an isolated primary sphere: 


i =| e (0) 2 (5) 
4m 

In what follows, byy we will imply the eigenvalue 
of y for a sphere determined by the exact equations, so 
that we satisfy the limiting case for the vanishing of the 
supplementary region. 

In order to simplify Eq. (5),we make use of the fact 
that the average values of t are multiplied by the weight- 
ing function (1 - e7 Lo), that is,the largest contribution in 
the average value t(Q) comes from the central region of 
the cross section S. We assume that t(Q) coincides with 
the value ty (Q) for the ray which passes through the cen- 
ter of the primary sphere. There is some error associated 
with this approximation, but in the general case there 
will be different signs for the different directions of Q and 
the error will average out in the integration over the total 
solid angle. It is shown later that this error does not 
change the sign of the total error. 

The final relation assumes the following form: 


1—y 
{= 


te (Q) €2 
Sa4h GE (6) 


4n 


Fig. 1, Diagram for the integration pro- 
cess, The primary region is shown cross- 
hatched, 


* For brevity, we will call the distances expressed in 
neutron mean -free path lengths the optical distance, 


The dependence of y on the optical radius of an is- 
olated primary sphere is shown in Fig. 2. The branch of 
the function yg(R) for dimensions of a sphere large com- 


pared with the neutron free path length,is obtained by the 


asymptotic diffusion method, which gives good accuracy 


in this region [2]. For R 1, the following relation from 


[1] applies 
Yo = & (0,78 — 0,39). 


Since both branches almost intersect each other, inter- 


polation between them in the small intermediate region 


gives reliable results. 
Equation (6) also applies when the density of the 
matter in the primary sphere is an arbitrary function of 


distance from the center of the sphere. We can convince 


ourselves of this if Eq. (6) is obtained directly from Eq. 
(2). This approach can be useful if the distribution of 
density in the primary sphere allows us to determine 
the quantity y) by exact methods (for example, if there 


is a spherical cavity in the center of the primary sphere). 


We may recall that Eq. (6) applies for an arbitrary 
distribution of density in the supplementary region. 
However, this relation can also be used for an arbitrary 


distribution of density over the total volume of the body. 


For this purpose it is sufficient to take as yp a value 
which is too high, corresponding, for example, to the 
maximum density in the primary region; we can also 
reduce the primary region to dimensions such that the 
density becomes essentially constant. For this purpose 
we must use the requirement that there be a maximum 


critical parameter, yp 


N 
wl 


optical radius, R 


Fig, 2. The critical parameter yo as a func- 
tion of the optical radius R for an isolated 


sphere. 


in the neutron density within the primary region. In the 
first approximation,the location of the maximum in the 
neutron density corresponds to a minimum in the quan- 
tity 
\ pr tteo IQ 
4x 
Errors of the Method 

The errors in this method come from two sources. 
The basic error yy , is caused by replacing the true 
neutron density distribution by a fixed quantity; this 
procedure leads to a value of y which is too high. The 
second error 6,7 is of the opposite sign and is due to 
the inexact averaging of the quantity t in Eq.(6). In 
order to determine the sign of the total error,we con- 
sider the behavior of these two errors as the dimensions 
of the primary region are changed. 

In the limiting case, when the primary region 
shrinks to a point about the maximum in the neutron 
density, the error 52y vanishes so that the averaging 
region for t (cross section of the primary sphere) ap- 
proaches zero while 6yy assumes its maximum value. 
This is the case for which the error in this method is a 
maximum. In the other limiting case, in which 
X= v,V = 1, both errors vanish. Close to this limiting 
case the error Say is of greatest importance because of 
the smallness of the first error. In order to show that 
the absolute magnitude of the error 5,y is smaller than 
537 when x “1, we find the derivative £05, “v) 

Ox Vv /x=1 
by comparison with the exact solution in the limiting 
case (1 - x)= AV/V<« 1, which we obtain by means 
of perturbation theory 3jt 

According to the perturbation theory analysis 


Ay AV 
ae eee 
where 
; = 1 
Ay =Y— Yo; k= = \oav (7) 


We find A;y from Eq. (3) which does not include 
the error yy. Considering the limiting case for the 
small supplementary region we obtain 


AV ( dQ 
AW = — Yo= \ Ue ati tad) regi tal) 


4n 


where the integration is carried out over the rays which 
emanate from the volume AV. For a primary region in 
the form of a sphere of radius R, as an approximation 
we obtain + 


+ This perturbation theory approach was developed in- 
dependently by N. A. Dmitriev in 1948 (cf. [4]). 


+ The accuracy of this approximation is discussed in the 
"two-sphere” section. 
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Taking account of Eqs. (4) and (9), we find that Eq. 
(8) assumes the following form: 


) oo 
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Similarly, from Eq. (6), including both errors we 
obtain 
Ay 2V 
Vi" 


AV, eee 
=hy ea Kemah y ° 


In Fig. 3 we show the dependence of the derivative 


a7 Ay 
ca ( ae 
on the constant y for the exact solution, the approxi- 
mate solution neglecting the error §,y, and for the so- 
lution which contains both errors. In the limiting case 
y « 1, the values of the derivatives are determined by 


the method which has been described in[1]. In the re- 
gion 0.3 < y< 1 we use the asymptotic diffusion method 


{2}. 
Comparing the derivatives of the first and second 
errors at xX = 1 


Om eOsVR Nes FOR GON Bre on oc 
peg) Ec himag oy )— BT he 


we see that the following relation always holds 


ees orn d ( Sev a bY , o2¥ 
a y DE fl ay ands, ( we av) <0, 


as was to be demonstrated. 


Fig. 3, The error derivatives for x = 1 as a function 

of y for the exact solution (curve k) the approximate 
solution without the error 5,y (curve k,) and for the 

solution which contains both errors (ky,»). 
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The total error increases as the dimensions of the 
primary region are reduced; the total error reaches a 
maximum when the primary region shrinks to a point, in 
which case the error due to averaging the neutron density 
becomes greatest while the error of opposite sign asso- 
ciated with the averaging of the quantity t vanishes. 
Hence, in taking solutions of actual cases it is convenient 
to choose the primary region to be a sphere of the largest 
possible dimensions, i.e., a sphere inscribed in the geo- 
metric shape of the body. In carrying out this procedure 
we must satisfy the basic requirement that the maximum 
neutron density is found within the primary sphere. 


Comparison of the Theoretical and Experi- 
mental Results. 


We compare the results obtained by the method de- 
scribed above with the experimental data for a system of 
different geometric shape of Oralloy - 93.5,which has 
been published in (5) ‘ 

The parameter ym is 0.74 for a sphere of Oy (93.5), 
in agreement with the experimental measurements of the 
mass of a critical sphere of Oy (93.5) of 51.3 kg. [5]. 

The neutron mean free path isa~!= 4 cm. 

Infinite Plate. We inscribe the primary sphere in the 
plate as shown in Fig, 4. Carrying out the integration in 
the right-hand side of Eq. (6) we obtain the critical equa- 
tion for the plate 


dV RES 
hasnay =1-+ Re® Ei (— R). (10) 


For the value y = ym = 0.74, from Eq. (10) we find 
that R= 0.77, i.e., the critical thickness of the plate 
H = 2Ra7!= 6.16 cm. 

The critical thickness of the plate found by the 
asymptotic diffusion method with the material constants 
in agreement with the experimental values of the criti- 
cal mass of a sphere of Oy (93.5) is 6.57 cm. The criti- 
cal thickness of the plate,obtained by extrapolation of 
the experimental data is H~6.1 cm [5]. 

Infinite cylinder. In similar fashion, for an infi- 
nite cylinder of Oy (93.5), from Eq. (6) we find R= 1.45; 
D = 2Ra*= 11.6 cm. 

The experimental value of the critical diameter D 
for an infinite cylinder of Oy (93.5) is 12.4 cm[5]. A 
calculation of the critical diameter of a cylinder by the 
asymptotic diffusion method also gives the value D = 
12.4 cm. 

Sphere. We estimate the radius a for a critical 
sphere of Oy (93.5) using Eq. (6). For this purpose the 
primary sphere of radius R is located eccentrically as 
shown in Fig, 5. Carrying out the integration in Eq. (6) 
for the supplementary region,we obtain; 


** Oralloy - 93.5 or Oy (93.5) is metallic uranium 
which contains 93.5% of the U?® isotope. 


for R==037 50 
for R= 0,048 


Cini 
a= 1,95, 


= 00; "= 7.0. ci; 


r= 7,8 cm. 


Experimentally, for a density of 18.8 g/cm’, the 
value of the radius r for a critical sphere of Oy (93.5) 
is 8.7 cm [5]. 

With a further reduction in R the error is reduced 
and the error changes sign when R < 0.37a. This feature 
serves to illustrate the fact that the requirement that the 
maximum in the neutron density be included in the pri- 
mary sphere is a sufficient condition but not a neces- 
sary condition for obtaining an error whose sign is favor- 
able in the sense of safety. 

Two spheres. For a system consisting of two iden- 
tical spheres, one of which is assumed to be the primary 
sphere, Eq. (6) assumes the following form: 


fat at-Z(t-g), gag) o'ae, ay 
Q 


where 02 = - bees Vithren( say) is the solid 


angle subtended by the second sphere at the center of the 
first sphere (b is the distance between the centers of the 
spheres). The function g characterizes the transparency 
of the sphere. ip 


Fig. 4. Position of the primary 
sphere in the plate. 


With fairly high accuracy we may assume that g 
is equal to the transparency factor go of a sphere for a 
plane neutron beam: 

4_e-28 (4. 9R 
goa \ eds = aes ) 
Ss 

Comparison with the results of the exact integration 
in Eq. (11) shows that replacing g by gy in the most un- 
favorable case for b = 2R leads to a negligibly small 
error which is less than 1% in the quantity 1 - git This 
result allows us to extend the analysis to the case of in- 
teraction of a large number of spheres. 

In Fig. (6) we show the dependence of the critical 
radius of the spheres of Oy (93.5) (in neutron mean free 
paths) on the relative distance between the centers of 
the spheres b/R determined by Eqs. (11) and (12) (curve 
a). For comparison,we also show the critical radius R 
as a function of b/R obtained by the method described in 
[6] (curve b). 


us ult (12) 


Fig. 5. Sphere with the ec- 
centric primary sphere, 


The relation given in Eq. (12) was verified experi- 
mentally by B. D. Stsiborskii and M. L Kuvshinov for 
two spheres of Oy (93.5) located close to each other. 
The results of the measurements are in agreement with 
those obtained by Eq. (11) to within the limits of the ex- 
perimental errors. The good agreement between the 
theoretical and experimental results is explained by the 
fact that the mean values of the neutron density in the 
primary and supplementary regions coincide because the 
spheres are identical. 

Equation (11) also applies for a system consisting of 
several spheres, In this case, we choose one of the 
spheres in the center of symmetry where the neutron 
density is a maximum as the primary sphere while the 
solid angle Q in Eq. (11) is summed over the remaining 
spheres (if they do not shield each other), 

System of a large number of spheres. We apply Eq. 
(6) to a system consisting of a large number of spheres 
with radii R which are uniformly distributed over the voi- 
ume of a sphere of radius a with a mean density of n 
spheres per unit volume. We take one of the spheres in 
the center of the system to be the primary sphere. 


critical radius, R 


20 24 28 G2 36 40 
relative distance between sphere centers, Ye 


Fig. 6. The critical radius R as a function of the 
relative distance between centers b/R for two in- 
teracting spheres: a) according to Eq. (11); 

b) using the method of [6]. 


++ The sign of this error is favorable as far as safety is 
concerned, 


ail 


Each sphere of the system can be characterized by 
an effective cross section for neutron interaction: 


= nwR?, 


= (1— gp) 5, (13) 


where gy is determined from Eq.(12), In this case the 
neutron mean free path length is 


pu ((lesso) |? (14) 


From Eqs. (13) and (14) we obtain the critical re- 
lation for a system of spheres: 


LY _ g=(1 a0) Bra, 
1—Yo 


or, in more convenient form 


1— Yo 


l—g)NV(>) =n, 


where N is the total number of spheres in the system. 

Because the primary sphere is located at the neu- 
tron density maximum, Eq. (15) contains an error which 
is favorable from the point of view of safety, that is to 
say it either reduces the total number of spheres in the 
system N as compared with the critical value Neyjz or 
reduces the radius of the system a for a given N. 

In order to estimate the magnitude of the excess we 
may note that the distribution of an average neutron den~ 
sity in the system as a function of radius is, in the first 
approximation, the same as the neutron distribution in 
an isolated sphere of fissionable material with a mater- 
ial constant y’ = 1/Q, where Q is the multiplication fac- 
tor by each sphere in the system for an incident neutron 
flux.ttHence, having determined the critical radius of 
this sphere R' by the same method as that used to deter- 
mine the critical radius of a system of spheres a, we ob- 
tain the same error in these two cases with respect to 
the critical radius, 

We locate the primary sphere concentrically in a 
sphere of matter characterized by a constant y' and al- 
low the radius of the primary sphere to approach zero. 
From Eq. (6) we find the lowest value of the critical 
radius of the sphere R’: 


(15) 


R'=|In(i—y’)]. 


Comparing the magnitude of R' with the exact va- 
lue of the critical radius of the sphere and introducing 
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the same relative correction in Eq. (15) we obtain the 
corrected critical relation for our system of spheres 


3 gyn (2) = ne 


ie iad—vT | ——s 
where y' = 1/Q, while R(y’) is determined, for example, 
from Fig, 2. 

We may emphasize, that because of the absence of 
experimental data,we must take Eq. (15) as a working 
relation and this provides a reliable excess. We can use 
Eq. (16) to estimate the magnitude of this excess. 

CONCLUSIONS 

The analysis which has been made of the above er- 
rors which arise and the comparison with experimental 
data verify the stability of the sign of the error as being 
independent of any given problem. This factor is of 
great value for solution of problems associated with 


” (16) 


safe handling of fissionable materials, 

The method described above is particularly effec- 
tive in calculations of the critical parameters for groups 
of bodies made from fissionable material, in which case 
other methods become very difficult. The method is ap- 
plicable for an arbitrary distribution of the density of 
the fissionable material. Thus, this approximation me- 
thod can be used for systems with a fissionable material, 
which is uniform in its properties, without moderators or 
reflecting shells. 

We wish to express our gratitude to B. D, Stsiborskii 
and M. L Kuvshinov for calling our attention to the ex- 
periment with the two interacting spheres; we also wish 
to thank V, A. Davidenko, Yu. N. Babaev and N. A. 
Dmitriev for discussion of the present work and a number 
of valuable comments, 
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REMOVAL OF OXIDES FROM SODIUM AND TESTS 
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Oxides present in sodium ,used as a heat carrier, accelerate the corrosion of the ducts and may 
cause blockage of the coldest sections of the circuit, The article gives the results of tests on 
cold oxide traps and work on an apparatus for determining the oxygen content in sodium. The 
results of the work can be used for planning experimental and industrial apparatuses in which 
sodium or sodium-potassium alloys are used as heat carriers. 


Until recently, the use of sodium as a heat carrier 
was difficult because of the absence of methods of re- 


moving oxides from sodium and also tests for the oxide con- 


tent, The presence of oxides in liquid sodium intensi- 

fies corrosion of structural materials (for example, stain- 

less steel — see Table 1[1]). When oxides are deposited 

from solution in the cold parts of the circuit, heat trans- 

fer to the liquid sodium or sodium-potassium alloy is 
lessened and the pipelines may become blocked [2]. 

The tests of apparatuses for determining the oxide 

content in sodium (indicators, oxides) and of oxide 
traps were carried out on a special testing installation 
(Fig. 1) with the following characteristics: 


Volume of metal in the circuit, ye 0.100-0,450 
Oxygen CONTENt, Wt. Ur w.5cccseseescseercees 0.1-0.001 
Working temperatures, “GC. ....svsseseenes. 110-550 
Pressure head of the pump, kg/ cm’ ,.,. 6-10 
Length of the pipelines, M ............++. ~50 
Diameters of the pipelines, mm.......... 25-48 


For preliminary heating of the metal all the lengths 
of the circuit were equipped with electric heaters, Dur- 
ing the operation of the testing installation the metal 
was heated by electric or gas furances, 

Apparatus for determining the oxygen content in 
liquid sodium. Methods for the chemical determination 
of the oxygen content in liquid sodium — mercury [3], 
butyl bromide [4], distillation [5], vacuum fusion [6] — 
are rather complicated, very time consuming and re- 
quire special equipment, When operations are carried 
out with radioactive sodium the chemical analysis is 
markedly complicated. The most difficult operations 
are sodium sampling and the transfer and introduction of 
the samples into the apparatus without contamination of 
the metal to be analyzed by oxygen. 


A stopper indicator method was used to determine 
the oxygen content. The idea is far from new [7, 8]. 


TABLE 1. Rate of Leaching of the Components of 
Stainless Steel in Sodium at 496°C 


! 


Rate a leaching,mg/ em*| 


Components of stain- |*™ont 


i n sodium 
less steel In os ee wae 
0,003 wt%o O wt. To O> 
italien ce OS ced | nS) 40 
Gobaltieeaeeac arse Pages OR? 4 
Tealita tiie ate vee fe | 1,0 4 
TD 26 


Piiangatiese: ie.) sre | 


It is based on the fact that if the temperature of the 
sodium current passing through a finely perforated dia- 
phragm is lower than the temperature of its saturation 
by oxides,the latter are deposited and block the holes in 
the diaphragm. This method does not require highly ac- 
curate apparatus and ensures the required accuracy with- 
in a wide range of oxide concentrations, 

Employing the relation between the solubility of 
oxygen in sodium and temperature 


W=2.7-104( a) (1) 


(where W is the solubility of oxygen, in wt%, t is the 
temperature, in °C) established from the data of a num- 
ber of works [9 - 14], the oxygen content in the metal 
can be found if the temperature at which the oxides are 
deposited from solution is determined, The deposition 
of the oxides is recorded at the moment the flow of me- 
tal through the stopper indicator, shown in Fig. 2, falls. 
Valve 1 with several radial grooves is installed so 
as to facilitate washing out of the oxides from the holes 
after the determination, The metal is cooled by air 
from a fan. A somewhat altered EPP-09 electronic 
self-+ecording potentiometer with a path/time scale of 
8 sec was used for simultaneous determination of tem- 
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Fig. 1. Setup of the testing installation for tests 
of apparatuses for determining the oxide content 
of sodium: 

1) pump; 2, 8) oxide indicators; 3) sampling de- 
vice; 4) heat exchanger (sodium alloy, sodium- 
potassium); 5) electric heater; 6) gas heater; 

7) sodium-water heat exchanger; 9) measuring 
tank; 10) cold trap; 11) pump tank; 12) magne- 
tic flowmeters, 


perature and flow of sodium (Fig. 3). The most charac- 
teristic parameters of different modifications of the ap- 
paratus are given in Table 2. 

The relations between the apparatus readings and 
the rate of flow of metal through the holes, the dimen- 
sions of the holes, the number of holes and rate of cool- 
ing of the metal were determined in order to select the 
optimum design, During these tests the oxygen content 
of the sodium varied from 0.002 to 0.1 wt%, and the 
temperature of the sodium from 110 to 550°C. The re- 


Unit A 


Wi 


Gir 
XH; 16 
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lations between the readings of the apparatus and the 
rate of flow of metal through the holes at two different 
oxygen contents were determined for each modification 
of the apparatus (for purposes of comparison, two de- 
signs were generally tested), The temperature at which 
the holes became blocked did not vary atflow rates of 
2.5 to 14 m/sec (Fig. 4). At lower rates of flow the 
readings are too low, the cause of this being unknown as 
yet. 


TABLE 2. Characterisrics of the Oxide Indicators 


Modification 


Design characteristics 


A B 


umber of grooves in the 
valve. ; 16 
imensions sof the hates: minsleto<t: 
iameter of the heat-ex- 


16 34 
1x4 |0,5x0,5 


Changer tube aes ays a ee 25 49 49 
Total area of the heat-ex- 
Changerin* soe le oes 0,454. O50 2th Orne 


When the size of the holes was changed from 0,5 x 
0.5 to 1 X 1 mm the readings of the apparatus were re-’ 
duced by 2.5%. The number of grooves must be suffi- 
ciently high (10-15) to eliminate the effect of chance 
blockage of the holes, The readings of the apparatus 
are independent of the rate of cooling of sodium under 
the following conditions of the experiment: constant 
oxygen concentration (0,008-0,02 wt%), variation in the 
flow of metal from 2.5 to 13 m/sec, rate of fall of tem- 
perature in the valve from 0.3 to 37° C/min. A cooling 
rate of more than 10° C/min is undesirable because it 


ee 
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Fig. 2. Design of the stopper indicator for oxides: 1) main valve with radial grooves in the valve, re- 
taining the oxides; 2) sodium-air heat exchanger; 3) flowmeter; 4) regulator valve; 5) thermocouple for 


measuring the temperature at the blocked point, 
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Time, min, 


Fig. 3, Examples of recording of the curves of the flow and tempera~ 
ture of sodium in the secondary apparatus of the oxide indicator: 
1) emf of the magnetic flowmeter; 2) temperature in the valve of the 


main valve. 


may lead to an error in the determination of the temp- 
erature at which blockage of the holes occurs, The ac- 
curacy of the indicator readings was checked by analy- 
sis of the metal accumulated after a specific period of 
operation of the trap. The results of chemical analysis 
and the determinations of the stopper indicator are com- 
pared in Table 3, which shows that they agree fairly 
closely. 

The majority of the components of the apparatus 
were tested for several thousand hours at a temperature 
of 300-500° C and showed quite satisfactory strength. 

Cold traps for removal of oxides from sodium, The 
working principle of cold traps[{7, 8, 15] is that the me- 
tal from the hot sections of the circuit is cooled and the 


Fig. 4, Relation between the temperature of block- 
age of the holes in the oxide indicator and the rate 
of flow in the holes (approximate oxygen content: 
O-~ 0.031%; O — 0.013%). 


oxides deposited from the solution are then filtered, 
Since the solubility of sodium oxides in sodium falls 
with a reduction in temperature,the degree of removal 
of the oxides is greater the lower the temperature of 
the metal in the trap. 

The trap (Fig. 5) is formed by a cylindrical contain- 
er 2, consisting of two shells — an outer shell O and an 


inner I — packed with stainless steel cuttings or wire, 
The container is equipped with a jacket 1, containing 
toluene for withdrawing heat from the sodium, The 
toluene boils and its vapor condenses in the coil 5, 
cooled by mains water. This complex cooling system 


TABLE 3. Oxide Content in the Trap, Determined 


Trapnumber| Data of the oxide] Gas- analysis data 
indicator 
{ “890+ 100 1000-+-500 
10 4750-700 6200900 


was used for reasons of safety in the event of leakage of 
sodium. The lower part of the trap had a settling tank 
6, in which the oxides are accumulated. The cone 7 
reduces the rate at which the metal flows through the 
settling tank, The trap is heated preliminarily by a ni- 
chrome heater located in the shell 4, 


The reduction in the oxygen content of the sodium 
passing through the trap can be calculated from the 
equation of the material balance: 


Vde= c—c’)dt, 
\ vQ ( a 


where V is the volume of sodium in the circuit, m?, c is 
the oxygen content in the sodium, wt%; c’ is the solu- 
bility of oxygen in the metal at a temperature t’, wt% 
(t’ is the minimum temperature of the metal in the trap); 
Q is the flow of metal through the trap, m? /hr; yis the 
specific gravity of the metal at the temperature of the 
circuit, kg/m’; 7 is the working time of the trap, hr. 
Equation (2) was derived from the following pre- 
mises: the volume of the trap is less than the volume of 
the circuit; the flow of sodium through the trap is low 
compared with the flow through the main circuit; all 
the oxides deposted from the solution remain in the trap, 
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i.e., the oxygen content in the metal leaving the trap 
corresponds to its solubility at the minimum temperature 


any : 
in the trap. a 
Integrating equation (2) we obtain the following ise Ls 
equation: 5 : : 


cme + (q—e)e* () 


where co is the oxygen content in the sodium before fil- 
tration; nis a number indicating the number of times the 
sodium present in the circuit has passed through the trap 
in time 7; n= Qr/V. 


Tate fata ye 
q 


——_ 


The following factors had to be determined when 
the trap designs were tested: 1) the maximum degree 
of purification provided by the trap; 2) the maximum 
amount of accumulated oxides; 3) the suitability of the 
design as a whole, its thermotechnical and mechanical 
characteristics, The tests were carried out as follows. 
The initial oxygen content in the sodium was deter 
mined before filtration, The metal was then circulated 
through the trap at a predetermined rate for a specific 
time. The oxygen content in the metal was then re- 
determined, This procedure was repeated until the oxy- 
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gen content had fallen to the required value, The re- Fig. 5, Design of cold trap with 
sults obtained were compared with the theoretical curve a capacity of 32 liters of sod- 
(3). If the oxides were removed from the trap the ex- ium, 
perimental curve appeared above the theoretical curve. 
By reducing the flow of metal through the trap,a Oxygen was added to the sodium either by mixing the 
position was reached where the experimental points co- pure metal with the oxidized metal or by introducing 


incided with the theoretical or were fairly near to them, sodium peroxide into the circuit. 


TABLE 4, Characteristics of the Traps Tested 


Metal capacity of the trap 


Conditions and results of the tests 


32 liter trap* 47 liter trap 

eStat MGs Mig wera oe. <caligeae Beak tat aetna @.o <6 ae 624 720 1900 
Clreilation time or the metal; hr.) ayedee. ss ano 105 130 336 
Flow ol metal pm ihe... Pain sss Osa codivbess amv 2 * the, fBOkAS 70 tO0L | Oe Sbes() 20nd Dee ate 
Rate of the metal, mm/sec; 

EDR EDETOULEL TSC LINN NT ee MR ES, on, SS arbaciae Voie 2°82—24 ee 

INRERMEN SWC LL ocr Aussie lesa « aPole-s. tei eu clocetteu lee an eee 2,74—2 ,0 2,9-0,6 
Temperature of the metal, °C; ; 

EARS ee CNSR EAS) Gee Ren Sty ER Aan Aer 2 ar 310—420 450—270 370—284 

ISTO MORN Dt ake. Sean tacct suena shiacs) concede) «4 Guha 140—217 200—120 226—160 

AlDemrAn-OUet eas wt ae ete ere SMewwad chert 148—280 300—180 480—125 
Pressure of the toluene vapor, atm: 

WY OLIN Dabo was egve-nhejery ropa rs Nehari antvane wae ss Way 0,5—3 0,0—2 0,0—2 

DIAMAR gen Reuss eee ery hoes Boch ott Hig neparae Meds 11 10 4,0 
Oxygen content in the sodium before operation of the trap, | 

wt."/o: 


UMC 8 Lesa Bees Wace», Le are yma ON Sap OSU Oe teers 0,044—0 ,02 0,04—0 ,02 0,03—0,015 


MAXIM ae Mean ecaly Gua atert et ee ore na bala ee 83 
Amount of oxides accumulated by the trap, g........ 3010 450 4850" 
Amount of oxid 
REET LEU CTS THET  rme leleRn 10 4 10 


* 13,2 wt. vol. % of oxides of metal in the trap was found in the settling tank, The determination was 


made with respect to the density of the metal (the specific gravity of sodium oxide is more than twice 
that of the metal itself), 
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As a result of tests of different variants of the traps 
(Table 4) the following facts were established: 

1) The cold trap can reduce the oxygen content in 
sodium to 0,002 wt. 

2) The required degree of purification can be 
easily obtained by selecting appropriate operating con- 
ditions. 

3) The more oxides accumulated the better is the 
filtration capacity of the trap. At higher flow rates of 
metal there is a considerable entrainment of oxides, 

4) The reduction of the temperature of the metal 
at the trap inlet, equivalent to a reduction in the mean 
temperature inside the trap, increases the filtration ca- 
pacity of the trap. 

5) The oxide capacity of the trap depends on its 
design, the type and density of the packing and the 
operating conditions, 

6) For the same flow rates of metal, packing in 
the form of cuttings retains oxides more effectively 
than 0.5 mm diameter wire. 

Only the inner shell was equipped with a packing 
(density 100 kg/fm*), The volume of the settling tank 
was increased and the residence time of the metal in it 
was also increased. These measures made it possible to 
increase the oxide capacity of the trap to 10%. 
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Fig. 6, Variation in the oxygen content in 
sodium during operation of the cold trap: 

) curve calculated from equation (3); 
the experiments were carried out at the fol- 
lowing flow rates of sodium through the trap 
(mYhbr): @) 0.168; A, J) 0.250; mg, O,x,@, ©, 
A) 0,285; @, ©) 0.300; O,@) 0.335, 


When the oxygen content of the sodium was 0.1% 
or more the traps were blocked very rapidly. In such 
cases the metal was first pumped through the trap with- 
out toluene, i.e., the trap was used as a mechanical fil- 
ter. After the oxygen content had been reduced to 0,06- 
0.05 wt%, toluene was added and filtration was contin- 
ued, The typical characteristics of the trap operation 
are given in Fig. 6. 

It should be noted that the operating conditions of 
the trap depend greatly on its oxide capacity. For ex- 


Height, mm 


Fig. 7. Temperature distribution a- 
long the height of the trap at differ- 
ent sodium flow rates (m/hr): D) 
0.140; O) 0.080; A) 0.060, 


ample, during a test with a flow rate of sodium through 
the trap of 0.080-0,117 m*/hr its oxide capacity was on- 
ly 4%; at flow rates of’ 0.270-0,500 mhr its capacity in- 
creased to 10%, This can be explained by the more uni- 
form distribution of oxides in the volume of the trap at 
high rates. 

To determine where exactly the oxides were accu- 
mulated, the inner and outer shells of the traps were 
flushed with gas after the tests and in some of them the 
content of the settling tank was analyzed. A compari- 
son of the results of the tests of different traps and ana- 
lyses of their contents and the contents of the settling 
tank after completion of the tests allows the following 
conclusion to be drawn with regard to the role of the 
different parts of the apparatus. The oxides are depos- 
ited from solution in the outer shell. The temperature 
distribution along the height of the trap was measured 
(Fig. 7). It was found that at a distance of 140 mm 
from the inlet the sodium temperature falls from 300- 
450° C to 180-200° C, i.e., nearly all the oxides are de- 
posited immediately after entering the trap. This ex- 
plains the fact, noted in [13], that the hydraulic resis- 
tance of traps increases as a result of blockage of the 
outer shell. But if the outer shell is packed nearly all 
the oxides are retained in it, The packing was, there~- 
fore, removed from the outer shell, 

The cooled metal passes from the outer shell to the 
upper part of the settling tank, forming 10-40% of the 
volume of the whole trap, The results of the tests 
showed that the amount of oxides retained in the volume 
of the settling tank is negligible. 

If dead-end and discharge lines are heated slightly 
they operate as diffusion traps, which can accumulate 
oxides or recycle them according to the oxygen content 


eal 


in the sodium and the ratio between the temperatures of 
the dead ends and the temperature of the circulating 
metal, 
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ON THE CHANGE IN THE COLOR AND TRANSPARENCY OF 
GLASSES WHEN BOMBARDED BY GAMMA RAYS 
FROM A Co% SOURCE AND IN A NUCLEAR REACTOR 


S. M. Brekhovskikh 


Translated from Atomnaya Energiya, Vol, 8, No. 1, pp. 37-43, January, 1960 
Original article submitted November 17, 1958 


Glasses of various types are widely used in radiation shielding. The problems involved in the 
effect of different types of radiation on glasses and the choice of glass with greatest radiation 
stability become important in this area, The present article deals with these questions,Results 
of experimental investigations where glasses of different types were subjected to bombardment 
from a Co’ source and to bombardment in a nuclear reactor are discussed, Glasses of different 
composition, products of several different glassworks in the USSR, were tested. Techniques for 
enhancing the radiation stability of the glasses are recommended in this article. Those charac- 
teristics of glasses which provide a fairly complete picture of the radiation response of the ma- 
terials, e. g., index of stability, index of saturation, minimum transparency, coefficient of rate 


of darkening, are introduced. 


Glasses darken or color upon exposure to radioac- 
tive radiations. The rate of darkening depends on the 
composition of the glass and the radiation dose; when 
radiation doses are large, transparency drops rapidly and 
the color characteristics of the glass are modified. Some 
data on the stability of industrial glasses to irradiation 
by y rays may be found in [1]. 

The present article presents the results of addition- 
al investigations of color and transparency changes in 
glasses upon irradiation using a Co® source in doses up 
to 10%, and upon in-pile irradiation with doses reaching 
pelOs I, 

Polished glass specimens 30 mm in diameter and 
5+ 0.1 mm thick were used in the tests. Spectral absorb- 
tion was determined in a quartz SF-4 spectrophotometer, 
and integral light transmission was found by using a tri- 
color VEI colorimeter. 


Change in Color of the Glass 


Figure 1 shows the curves of the spectral transmit- 
tance of various glasses found prior to irradiation and 
following irradiation with a dose of 10®r from a Co™ 
source. Glass No. 10 from the "Krasnyi luch" works 
(Fig. 1a), used for the manufacture of extruded items of 
optical glass; polished (rolled) glass PG from the Dzer- 
zhinskii glass works at Gusev (Fig. 1b); and green glass 
No. 18 (Fig. 1c) from the "Krasnyi luch" works \ were 
used, Both specimens of colorless glass darkened intense- 
ly and took on a dark brown coloration, 

Figure 2a shows the curves of spectral transmittance 
prior to and following irradiation of red glasses colored 
with selenium and cadmium sulfate (selenium ruby 


glass), Fig. 2b applies to orange cadmium OS-6 glasses, 

while Fig. 2c refers to blue (cobalt) glasses, No. 17. 
Table 1 gives the color characteristics of red glass- 

es of the "copper ruby" type prior to and following irra- 


=) 
S 


S 


Transmittance, % 


i>) 


Transmittance, %o 
ae 
es8SSBsss 


Transmittance, % 


308 410 500 600 700 600 900 4000 
Wavelength, mm 


Fig. 1. Curves of spectral transmittance of 
glasses: (----) prior to irradiation; ( ) 
following irradiation with a dose of 108 r, 
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diation with a dose of 10°r. We see from the Table that 
the color characteristics of the glasses suffered neglig- 
ible change. 
Change in Transparency of the Glass 
Each glass specimen was exposed to various 

gamma-ray doses ranging from 10? to 108 r, emanating 
from a Co® source, as well as irradiation in a nuclear 
reactor, in doses from 108 to 5:10! r. 

Figure 3 shows the variationinthe transparency of 
conventional window glass, OG from the Gorky glass~ 
works (1), polished PG glass from the Gusev works (2), 
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and colorless No. 10 glass, for extruded items, from the 
"Krasnyi luch" works (3), as a function of radiation dose. 
The shape of the curves in Fig, 3 demonstrates that soda- 
lime glasses of compositions generally accepted in in- 
dustry are stable to gamma-rays for doses of up to 104 r. 
Only at that dose limit is appreciable darkening of the 
glass observed, Further irradiation results in rapid dark- 
ening; within the 10°10’ r dose range, the slope of the 
curve dips, and at 10*° r the limit of color saturation sets 
in. For polished glass, color saturation sets in at lesser 
doses, at about 108 r, 


ie ; 
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3 b c 
Fig. 2. Curves of spectral transmittance of glasses: (----) prior to irradiation; 


) after irradiation, dose 10 f. 


( 


The integral transparency of window glass irradiated 
with a maximum dose of 6-10!" r is 7.1%, and for polish- 
ed glass it is 12.8%, while it is 6.2% for glass No. 10. 

Figures 4 to 6 show the curves characterizing the 
change in transparency of the glasses tested in response 
to irradiation with y rays. Inspection of these curves 
shows clearly that there is an inflection point for each 
glass corresponding to the gamma dose at which appreci- 
able darkening or an appreciable decrease in the trans- 
parency of the glasses begins. 

The results of the test show that when a cation of 
larger ionic radius in the glass is replaced by a cation of 
smaller ionic radius, the stability of the glass to radia- 
tion bombardment is enhanced, 


_ It is plausible to assume that the substitution of the 
Mg?* ion with its radius of 0.78 A for the Na‘ ion of 
radius 0.98 A and the Kt ion of radius 1.33 A explains 
the considerable stability to radiation of the rolled glass, 
compared to No, 10 glass. 


Terminology 


Hitherto, there has existed no unified terminology 
with respect to the determination of stability of glass to 
the effects of radioactive radiations. For the stability 
index (Ky) of the glass, it would be feasible to adopt the 
value of the logarithm of the dose, Dy, at which appreci- 
able darkening of the glass sets in, i.e., Ky=log Dgas- 

A decrease in transparency to visible light of 3-5% will 


Table 1. Color Characteristics of Glasses of the "Copper Ruby” Type Prior to and Following Irradiation with 


a Dose of 10° r 


] ; } 4 
. riietnereminatee ony Color coordinates Color sat-| 
Glass species of glass g |length X,\in visible . uration, %o, 
mm THT range 
Unglazed No, 13 copper ruby, from 
"Krasnyi luch"; 
PRlOr torisradia ions. «cts steen a teeke es G 620 Dae 0,695 0,305 100 
followingurradiationn Ss)... 4.0.7. = 0,9 0,700 | 0,295 100 
SKSG copper ruby, from the Chernyatinsk works 
PLLOLMOMITLAG AON aes ees ene ney ee 5 4 638 20) Ondo 0,284 100 
following irradiation........... cee ‘ ‘i 0,8 0,705 0,285 100 
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Fig. 3. Variation in the transparency of colorless 
glasses upon irradiation. 


LEPTIN hee whe oT aad 
Logarithm of gamma dose 

Fig. 4. Variation in the transparency of 
glasses in response to gamma irradia- 
tion. 1) Opal No. 14 glass from “Kras- 
nyi luch"; 2) same type glass from 
Chernyatinsk works; 3) blue No. 17 
(cobalt) glass from "Krasnyi luch"; 
4) red No. 13 (copper ruby) glass from 
“Krasnyi luch”, 


Logarithm of gamma dose 


Fig. 5. Variation in the transparency of colored 
glasses in response to irradiation: 1) Blue-green 
SCh glass from Chernyatinsk works; 2) green 
5A/Z glass from “Kransyi luch"; 3) red SKSG 
glass from Chernyatinsk works; 4) red K4/2 
(selenium ruby) glass from *Krasnyi luch™, 


be adopted as criterion for the process of appreciable 
darkening. 

As is clear from the graphs inFigs, 4-6, the rate of 
coloration of the glass increases continuously as the irra- 
diation dose is stepped up. However, there is a certain 
limiting value to the irradiation dose for each glass spec- 
ies, any increase above which will fail to produce any 
further intensification in coloration. 

This limiting value D,,, will not be the same for 
each glass species, and it would be advisable to introduce 
still another term to evaluate the response of the glass to 
radioactive bombardment; the saturation index (Kyat) 
which is of the form K,,, =log D,a,. Furthermore, the 
minimum value of the total light transmittance T_ at 
the saturation dose Deat is highly significant in the case 
of glass, 


Transmittance, % 


Fig. 6. Variation in transparency of colored glass- 
es upon irradiation: 1) Orange OS-6 (cadmium) 
glass from Chernyatinsk works (unglazed); 2) same 
glass (glazed); 3) experimental US-4 glass, sulfur 
colorant; 4) red K 3/2 (selenium ruby) glass from 
"Krasnyi luch” works; 5) red Ks-1 (selenium ruby) 
glass from Chernyatinsk works. 


Still one more criterion, the darkening rate or co- 
efficient of rate of darkening Q, will be found useful in 
rounding out the description of the glass response, this 
quantity being accepted as the ratio of the transparency 
factor Ty prior to irradiation to the transparency factor 
T following irradiation, at a dose of 10¢r, i.e., 


The irradiation dose of 106 r is fully adequate for prac- 
tical purposes. 

These four quantities : the stability index Ky, the 
saturation index K,,;, minimum transmittance T,,, and 
coefficient of rate of darkening Q, give a fairly com- 
plete picture of the response of glass to irradiation. The 
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TABLE 2. Values of the Coefficients Characterizing the Response of Glass to Radioactive Bombardment 


Glass species 


Polished (rolled) PG glass, Gusev glassworks 


No, 14 opal glass, "Krasnyi luch” 
OPLS opal glass, Chernyatinsk glassworks 
Blue No. 17 (cobalt) glass, "Krasnyi luch" 
Blue-green SCh glass, Chernyatinsk works 


Red (selenium ruby) K 4/2 glass, "Krasnyi luch" 


values ofthese quantities are entered in Table 2, 

It may be ascertained from the data in Table 2 that 
colored glasses with such colorants as selenium, cadmium 
sulfate, cobalt oxide, copper oxide and bichromates 
possess a high stability index and a low coefficient of 
rate of darkening, However, for irradiation doses of the 
order of 9-1010 r, they lose virtually all of their trans- 
parency in the visible portion of the spectrum. 

The investigation presented in[ 2] contains an inaccu~ 
rate statement that, on receiving doses of radiation close to 
10'°r, the glass shatters. Not the least sign of shattering of 
the glass was observed in our experiments. 


On Enhancing the Stability of Glass to Radiation. 

The character and intensity of changes in the color- 
ation of glass are dependent on the composition of the 
glass and the irradiation dose, and not on the form of 
the radiation. Industrial glasses, with the exception of 
special shielding glasses, absorb y rays and neutrons 
weakly, and the effects of these forms of radiation there- 
fore result in coloration of the glass through the entire 
thickness of the glass. The depth of penetration of 
particles and 6 rays into the glass is negligible, and 
changes in glass coloration are observed only in the 
superficial layers. All forms of radiation lead to the 
formation of free electrons and excited atoms in the 
glass. 

In the presence of ions free to change valency, the 
mechanism involved in the change in glass coloration 
may be explained by the attachment or ejection of 
electrons, as for example: 


V** (green). -| e7 —>V*" (pink), 


In the case of some ions, as ce4t, Ass sb5* , and 
Pb2*, the change in valency fails to bring about any 
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change in the coloration of the glass: 


Ce’ (colorless) + e— -——> Ce** (colorless). 


Transition from Ce** to Ce3* does not result in any 
change in the absorption curve over the visible portion 
of the spectrum, since the absorption band for CeO, and 
Ce,O3 lies in the ultraviolet region. 

Polyvalent ions in glass assume the role of electron 
acceptors, accepting electrons set free in response to 
radiation. 

In glasses containing polyvalent ions, the energy of 
radiation is spent primarily in changing the valency, 
rather than in the formation of color centers. Such glass 
species accordingly exhibit high radiation stability. 

The property of increasing the radiation stability of 
glasses is also possessed by cr3+, Mn*, As5* and Sb 
and Fe3* ions. Ceric oxide is the only additive which has 
found practical use in industrial glasses as a material for 
sharply increasing the stability of glass to radiation. The 
intense effect of ceric oxide on the stability of glass to 
radiation was first disclosed as far back as 1926 [3], and 
in subsequent work [4] and [5]. Some very important 
work, in both theoretical and practical areas, has been 
done in this field by V. V. Vargin and his associates at 
the State Optics Institute. 


At the present time, ceric oxide is used in amounts 
of 0.4-1.2 wt% in the manufacture of glasses possessing 
increased radiation stability. 

The status of shielding glasses with enhanced stab- 
ility to y rays is a more complex one. High-density 
glasses with a high lead oxide content are used as trans- 
parent shielding. Such glasses possess increased radia- 
tion stability in and of themselves. For example, the 
absorption coefficient of glass containing 69.5 wt% lead 
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Glass species 


Oxides 

m10] PG K3/2 | S-1 we 18 | 5P/3 | SKSG | mi4 | M47 |: 413 
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Ore) ap are a 8,4 O70 — os. 9,31 8,58 7,94 Ooo 0,34 
MgO 0,15 3,3 a te 0,24 0,14 3,76 0,19 | 0,25 0,24 
Na,O £7563 1) 1392 2,64 16,0 11,14 135.76 “lod 7 Jd 46.3 14,7 
Ke Oe ee 1,24 = 12,44 3 9,27 = Gls nt 0,92 6,5 
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Cds ES me 432 0, 2 = as c= 28 = 
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oxide is 1.48 cm ! following an irradiation of 106 r, 
i.e,, it is significantly less than that found for window 
and polished glasses. However, the ensuing darkening 
is inadmissible for a number of cases. 

Ceric oxide additive in lead glass does not produce 
such noticeable results as in silicate glasses. In order to 
improve the y stability of high-lead content glasses, 
arsenic and antimony are used as additives. 

The transition of the As* ion to As® and the Sb*'to 
Sb® transition calls for a high expenditure of energy 
above that required for the Ce* to Ce® transition. It is 
however possible that, when arsenic trioxide is added in 
combination with other polyvalent oxides, highly stable 
glasses may be successfully produced. 

Cognizance must be taken of the desirability of 
replacing ions of larger dimensions by smaller ions, in 
working out the compositions of y-stable glasses, The 
replacement of SiO, by PyO;; and B,O, in glass improves 
the stability of the glass. Replacement of K,O and Na,O 
by LizO and MgO also has a salutary effect. 

Our investigations resulted in the elaboration of a 
new technique for improving the y stability of glass by 
addition of 0.5-1.5wt% elemental sulfur. Electron at- 
tachment with the formation of polysulfides is postulat- 
ed as an explanation. 

Figure 7 shows the spectral transmission curves for the 
experimental glass US-4. The composition of this spec- 
ies is as follows (in wt.%): SiO, (71.5); B,O3 (1.5); CaO 
(8.5); MgO (1.5); NagO (17.0); sulfur above 100 wt,% 
(1.2). 


Experiments have shown that the color of US-4 
glass does not change in response to irradiation deses of 
3-108 r. In this species of glass, the stability index Ky= 


=5, the saturation index K,,,;>10.7, the darkening rate 
Q@1.0. The transparency of the glass following irradia- 
tion with a dose of 5-10" r is 27.7%, i.e., much higher 
than in the case of all other glasses investigated. 

The processes involved in the change in coloration 
of the glass in response to irradiation have a reversible 
character, Coloration of the glass may be partially or 
completely reversed by heating the glass to temperatures 
close to the softening point, by exposing the glass to sun- 
light, or by exposing the glass to flourescent lamps. 

The first process proceeds very rapidly; a prolonged ex- 
posure is needed for sunlight to work its effect on the 
glass. 

Polished glass from the Gusev works acquired, 
following y irradiation with a dose of 5:10*r, a faint 
brownish-gray coloration, which disappeared completely 
when the specimen was left standing for four to five days 
of exposure to ordinary sunlight. The same glass, irradi- 
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Fig. 7. Spectral transmittance of experimental glass 
US-4: - — —-—) prior to irradiation; ) after 
irradiation, dose 10° r; — x— ) after irradiation, 
dose 3 + 10° r, 
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ated with a dose of 106 r, was successfully bleached only 
after 40 days of exposure. 

Under the service conditions appropriate to shield- 
ing glasses, in the presence of low source activity, the 
rate of darkening of the glass is equal to the rate of fad- 
ing, so that the transparency of the glass suffers no 
change. 

In the manufacture of shielding glasses, it is some- 
times necessary to glue separate plates of glass together 
into a sandwich, When the refractive indices of the 
glass and the film in the sandwich are approximately 
equal, reflection losses (4-6% on each surface) are ab- 
sent, so that the transparency of the glued laminates is 
much higher than would be the case in a window made 
up of separate plates of glass with air layered in be- 
tween. 
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Glass specimens glued to polyvinyl butyral film were 
subjected to y irradiation in a nuclear reactor to study 
the y stability of adhesives. As the experiments showed, 
polyvinyl butyral film exhibits excellent ae) to 
y rays, withstanding irradiation doses of 107 r without 
alteration. Bubble formation ensues when the film is 
exposed to a dose of 10° r, 

In conclusion, I would like to express my acknow- 
ledgment to G. Ya. Vasil'ev and Yu. Ya. Zhelazova for 
their kind assistance in performing the experimental re- 
search discussed herein. 
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The production of currents having a maximum re- 
lative content of either protons or ions of molecular hy- 
drogen, from hydrogen ion sources, is of significant in- 
terest for a number of physical problems. Preliminary 
mass~spectrometric studiesof ions emerging from the slit 
source of the 1.5 m cyclotron of the Academy of Scien- 
ces of the USSR have shown that ion currents of the order 
of 60 ma, containing 95% protons or 80% molecular hy- 
drogen ions, can be obtained from an area 20X2 mm [1]. 
Later, by employing optical-spectroscopic techniques, a 
successful attempt was made to trace not only the ratio 
of H* to Ht ions in the ion source, but at the same time 
to monitor variations in the relative content of neutral 
H and Hp particles in the plasma of the discharge. An 
ion source with its associated electrical supplies are di- 
agrammed in Fig. 1. A three-prism glass Zeiss spectro- 
graph, with camera objective of focal length 840 mm, 
was used in investigating the optical spectrum of the hy- 
drogen. The reciprocal linear dispersion of this instru- 
ment in the region of 6500 A was 38 A/mm. 

Light from the ion source emerges from a rectangu- 
lar opening 15 mm high and 2 mm wide, oriented along 
the discharge channel. The first line of the Balmer ser- 
ies (6562.79 A) was used in studying the emission of 


Gas admitted to 
charge channe 
Gas admitted to cathode region. 
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Fig. 1, Diagram of ion source. 


monatomic hydrogen. The half-sum of the intensities of 
lines 6031.90 A and 6018.29 A were used as measures of 
the emission of molecular hydrogen, Spectral-line in- 
tensity may be used to advantage to estimate the con- 
centration of neutral particles. It should however be 
borne in mind that this intensity relates to the number 
of excited particles at a given energy level [2]. Under 
conditions of constant excitation, variation in the rela- 
tive intensity of the spectral lines of monatomic and 
molecular hydrogen may be used to monitor variations in 
the relative concentration of neutral particles Nq/ NH: 
When the conditions of excitation are varied, the rela- 
tive intensity Jp;/ Ju, characterizes the ratio Nyyo4;/ 

/ Ny oH, where 6;, and 6;,, are the effective excitation 
cross sections for H and Hy. Since the value of 9,, var- 
ies only slightly over the range of electron energies per- 
tinent to the investigation [3], the variation in the rela- 
tive intensity Jjy/ Jy, will in this case reflect the varia- 
tion occurring in the relative concentration of neutral 
particles in the hydrogen discharge. 

During the experimentation, a sharp difference be- 
tween the yields of H* and Hj ions from the source, de- 


Dy 


S 


Aa) 


AR 
rn 
N 


MLS 
A 
BEBe tia 


9 
é 
ri 
6 
5 
5 
3 
2 
f 
0 


Be 
2 a 


g, cm’/ min 


Fig. 2, Ratios of ion currents 1 and of line in- 
tensities 2 plotted against rate of gas influx q, 


at In= 3 amp, Ub = 160 v. 
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Fig. 3. Variation in the absolute values of 
ion currents and spectral-line intensities as 
a function of gas influx, at I= 3 amp, ' 
Up= 160 v. 


Fig. 4, Ratios of ion currents 1 and line inten- 
sities 2 versus discharge current, for the case 
where gas is admitted to the discharge channel 
(q=4 cm3/min , Up=180 v ), and where gas is 
admitted to the cathode region (q=2 cm3 /min, 


Up = 160 v). 


pending on the point at which the gas was admitted, was 
disclosed. Figure 2 shows the ratios of ion currents ipgt/ bee 
and line intensities J,;/J,,, plotted as functions of gas in* 
flux where two different methods were used for introduc- 
ing the gas. The points plotted in all cases correspond 
to admittance of gas to the cathode region (hollow cir- 
cles), while curves with crosses correspond to admittance 
of gas into the gas-discharge channel. The shape of the 
curves of relative intensity J/ Ju, characterizing the dis- 
sociation process show that dissociation is greater when 
the gas is admitted to the gas-discharge channel. To 
shed light on the behavior of curves igyt/ ip;+ and Jpy/ Jyz,, 
the absolute values of the ion currents eth, and spec- 
tral-line intensities J, and Jup are plotted in Fig. 3 as 
functions of gas influx. For greater clarity, the intensity 
of the lines in the first spectrum is assigned a value of 
unity, and the intensity of lines in subsequent spectra is 
referred to it. Plotted in Fig. 4, for comparison purposes, 
are curves expressing the relation of ratios ipgt/ tut and 
Ty/ Ju, to the discharge current when the gas is admitted 
to the gas-discharge channel or to the cathode region, as 
the case may be, at the point of peak difference between 
the yields of Ht and Hj ions. 

In the case where gas is admitted to the gas-dischar- 
ge channel, the curve for the line-intensities ratio shows 
that the relative concentration of monatomic hydrogen 
increase quite steeply as discharge current rises, Like- 
wise, the relative concentration of H+ and Hf ions as a 
function of discharge current also shows this behavior. 

When the gas is admitted to the cathode region, dis- 
sociation of hydrogen molecules and the relative concen- 
tration of protons at first decline in response to increased 
discharge current, and then later rise. The variation in 


Fig. 5, Dependence of ratio of ion currents 
and line intensities on discharge voltage, 
when gas is admitted to discharge channel 
(q=4 cm3/min , Iy=2 amp ) and when gas 

is admitted to cathode region (q=2 cm3/min, 
Ip=2 amp ). 


the absolute values of the spectral line intensities and 
ion currents demonstrates that the minimum in Fig, 4 is 
due to a sharp drop in the concentration of monatomic 
hydrogen, resulting in a relatively high yield of molecul- 
ar hydrogen ions. Note that a decline in H and Ht with 
increased discharge current is characteristic for the case 
when gas is admitted to the cathode region. The reasons 
for this effect are obscure. 

Figure 5 shows the dependence of the ratios mention- 
ed on discharge voltage. Curves giving the ratios of the 
intensities of spectral lines are similar in shape to the 
ion-current curves; the curve of relative intensity for 
gas admitted to the discharge channel reaches saturation 
at about 200 v, while the curve for the ion currents con- 
tinues upward. 

In the case where gas is admitted to the cathode re- 
gion, the curve of relative intensity of spectral lines 
shifts by about 40 v in the direction of lower voltages. 
This is an indication that the most complete dissocia- 
tion sets in earlier in the plasma than does ionization, 
as the voltage increases, Spectral-line intensities and 
absolute ion-current values show that the variation in 
the degree of dissociation of molecular hydrogen and in 
the relative concentration of protons as a function of dis- 
charge voltage are associated with the increase or de- 
crease in the concentration of monatomic hydrogen in 
the plasma. 

The experimental material reviewed here shows 
that the degree of dissociation in the ion source, as like- 


wise the proton yield and yield of molecular hydrogen 
ions, are affected by the point of admittance of the gas, 
An important role is assigned to such discharge constants 
as gas influx rate, discharge voltage, and discharge cur- 
rent [4]. The similar behavior of the curves characteriz- 
ing the dissociation and ionization processes accompany- 
ing variations in the discharge constants bears evidence 
that the process of proton formation has a two-stage 
character: dissociation of Hz takes place first, followed 
by ionization of the hydrogen atoms. The large proton 
yield in the case where gas is admitted to the gas-dis- 
charge channel is explained by the fact that a larger 
number of molecules engages in interaction with the 
electrons in that case than when gas is admitted to the 
cathode region. Another possibility is that, when the gas 
is admitted to the discharge channel, the gas particles 
interact with slower electrons; the process of dissociation 
of the molecules would be more effective in that case 


[5]. 
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NEW ISOTOPES OF HOLMIUM AND ERBIUM 


N. S. Dneprovskii 


Translated from Atomnaya Energiya, Vol, 8, No. 1, pp. 46-47, January, 1960 


Original article submitted August 13, 1959 


The conversion-electron spectra of isotopes of the 
erbium fraction, formed upon bombardment of tantalum 
by protons of 660 Mev energy on the Dubna synchro~ 
cyclotron, were investigated, 

The spectra were studied on a BPP-1 magnetic beta- 
spectrograph with a focusing angle of mV2 [1]. When a 
source 1X30 mm was used, with solid angle 0.4% of 4 7, 
the resolution of the instrument at the K-661.6 line of 
Bal? was 0,2-0.25%.. The magnetic field in the spectro- 
meter was measured by means of a special magneto- 
meter, thus ensuring efficient performance and making 
it possible to study conversion lines with a half-life of 


15 min or longer. Anaccuracy of 0.1-0,3% was achieved in 


determining the energy of the conversion lines 

A group of lines with half-life Tj, = 2.4 hr was 
discovered from among the host of conversion lines (over 
100) belonging to the neutron-deficient isotopes of erbium 
and holmium [2]. The investigation of these lines made 
it possible to identify several transitions occurring both 
in the holmium nucleus (Z= 67) and in dysprosium nucleus 
(Z= 66) [cf Table]. 

Upon isolation of the daughter holmium isotope 
from the erbium fraction 2 hr. after the latter was isolat- 
ed from tantalum, the intensity of lines corresponding to 
y transitions of the dysprosium nucleus had dropped one- 
half within (27+ 2) min. 

From the facts accumulated, we may ascertain the 
presence of the decay chain 


2,4hr 
Ere oho 


27 min 
pace WORE 


The ratio of energies 98.6 and 218.2 kev and of intensi- 
ties of the transitions support the inference that these 
transitions take place between levels F,=98.6 kev (2") 
and E9= 316.8 kev (4°) of the first rotational band of the 
even-even dysprosium nucleus. The position of the first 
excitation level as a function of neutron number [3] sug- 
gests that the mass number of the nuclei belonging to the 
above decay chain is A= 158. 
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Gamma Transitions in Erbium and Holmium Nuclei. 


Charge on 
hy, nucleus Multiple 
kev hee under going order 
transition 
98,6 | 2,4+0,4 66 E2 
BAS. Vir 27548054 66 es 
320,59 2,4+0,1 66 £2 or M3 
357 te eed €6 = 
387,3 | 2,4-.0,1 67 Ja 
945,9 -| 2,5-+-0,4 a = 
848,5 | 2,5+0,4 is “= 


The values of the constants A and B in the equation 
for energy levels E = AI(I+ 1)— BI’(1 + 1) determined 
according to the position of the first two levels, are 
equal, respectively, to 16.7 + 0.35 and 0.042 + 0.008. 

I should like to avail myself of this opportunity to 
express My warm gratitude to B.S. Dzhelepov and K. Ya. 
Gromov for their kind counsel and discussion of the ex- 
perimental results, and to I. A. Yutlandov and Yu. V. 
Narseev for performing the experiments on the chemical 
isolation of the specimens. 
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FISSION CROSS SECTION OF Th?229 FOR MONOCHROMATIC 


NEUTRONS IN THE 0.02 -0.8 ev REGION 


Yu. Ya. Konakhovich and M. I. Pevzner 


Translated from Atomnaya Energiya, Vol. 8, No. 1, pp. 47-48, January, 1960 


Original article submitted June 4, 1959 


Until recently, studies of cross sections of fission- 
able substances had been largely limited to the three 
most easily available isotopes of greatest practical im- 
portance, viz., Tae Claes and Pu®9, Analysis of data 
obtained on the cross sections of fissionable isotopes [1], 
has revealed some interesting regularities in the distri- 
bution of level spacing between levels of fission and 
neutron widths [2-4]. 

One feature common to all the fissionable nuclei 
studies up to now, including Pu! and Am™! in addition 
to those enumerated above, is the large value of the 
fission cross sections in the thermal region, varying in 
obedience to the inverse law 1/VE, and not interpretable 
in terms of the resonances lying close by to the right, 
This "1NE background” is usually attributed to the effect 
of the "bound" level lying to the left of the neutron zero 
kinetic energy. The fact that a resonance appears in the 
energy range 0.2-0.3 ev in the cross sections of all iso- 
topes fissioned by thermal neutrons which have been 
studied to date merits attention, although there is no 
difficulty in demonstrating that the statistical probabili- 
ty of such a coincidence is very meager. The number 
of isotopes investigated must be increased in order to 
throw some light on the. underlying aspects of the behav- 
ior of the fission cross sections, However, an obstacle to 
measuring fission cross sections is the intense @ activity 
of most isotopes fissionable by thermal neutrons, result- 
ing in the effect of pulses triggered by alpha particles 
being superimposed on the fission counter, while measure-~ 
ments of total cross sections require a large number of 
different substances, inaccessible at present for most 
isotopes. 

One of the most convenient nuclei for research pur- 
poses is the Th?” nucleus, an &~active nuclide with a 
half-life of 7,300 years [5]. 

The fission cross section spectrum of Th?” over the 
neutron energy range 0.025-0.80 ev was studied, using a 
crystal neutron spectrometer [6] mounted on the horizon- 
tal beam of the IRT reactor. An ionization fission cham- 
ber containing the material under study and a "thin" pro- 
portional counter loaded with boron trifluoride were plac- 
ed in series in the path of the monochromatic beam of 
neutrons, The cross section for the B!°(n,«) reaction 
obeys the law 1/VE over a broad range of energies, and 


the relative values of the fission cross sections may be 
found from the formula 


04 = const ty ) < 
VE\ Ns 

where Np}, and Np are, respectively, the counting rate 
of the chamber and proportional counter for monochro- 
matic neutrons, after deducting background. In our 
arrangement, there was 6- 107 g Th?”? obtained by chem- 
ical separation from U?33, in which it has accumulated 
as aresult of a decay. 

The purity of the preparation we used is discussed 
in an earlier paper [7]. According to the data reported 
therein, uranium impurities did not exceed several hun- 
dredths of a percent. The preparation was coated in 
oxide form on platinum foil 7“ thick, The geometry 
of the layer was close to rectangular, with sides 20 X 26 
mm, The amount of Th?® present was determined by 
the total @ activity of the layer, with the four a decay 
steps following decay of Th229 accounted for. The foil 
bearing the preparation was placed on the high-voltage 
electrode of the plane ionization chamber, which was 
filled with a mixture of argon and carbon dioxide (8%) 
to a pressure of 500 mm Hg. Since the cross section of 
the diffracted beam (8X25 mm) was smaller than the 
dimensions of the layer, the electrodes were placed at 
a 30° angle to the beam path, favoring effective use of 
the material. On account of the very low counting rate 
for fission events, the time required to reach 10% statis- 
tical accuracy fell within the range of 50-100 hr for one 
point. Special measures were therefore utilized to 
cope with noise, and measurements of effect and back- 
ground were alternated, 

The diagram shows the variation of the fission cross 
section of Th29 over the energy interval 0,02-0,80 ev, 
with some indication as to the resolution of the instru- 
ment. The curve is normalized to the value of =(4521)) 
barns[1]. This normalization must be regarded as only 
tentative, since the value of 45 barns apparently has 
reference to the reactor neutron spectrum, while the 
variation of the cross section of Th?” does not follow the 
1/VE inverse law in the thermal region. The energy 
scale of the instrument was calibrated according to the 
level of iridium at 0.654 ev, by measuring the total 
cross section. A correction for higher-order reflection 
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Fission cross section of Th?”; 
O) crystal spectrometer; x) mechanical chopper. 


effects was introduced in the 0.025-0.07 ev region. 
Three resonance levels are manifest inthe cross section, 
at energies 0.240, 0.563, and 0.727 ev. 

The magnitude of the area S under the cross section 
curve was used to determine the product of)G = 2S/m for 
the 0.240 and 0,563 ev levels, which were 3447.5 and 
74 3.3 ev-barns, respectively. ~ Ane 

The resonance integral for Th?®, J-pp= \ 04 (E) Tp? 


0,4 
determined from the value of the cadmium ratio Rog 
in accord with the familiar formula 
Rev 4 

cad 
7 
; The ae , 
was 240 barns, since RCg= 2.34 and Rcg= 15.2 in the 
spectrum of the IRT reactor. 


ean Orn thermal 
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The resonance peaks at 0.563 and 0.727 ev provide 
an explanation for only a small part of the observed re- 
sonance integral, 

Measurements were made for control purposes near 
the resonance peak at 0.240 ev. The chamber counting 
rate found when the beam was masked by a cadmium 
layer 0.8 mm thick was equal to the counting rate for 
an incoherently scattered beam; this attests to the fact 
that the 0.240 ev resonance is due to first-order reflec- 
tion neutrons, 

In the 0,02-0,26 ev energy region, the fission cross 
section of Th? was measured by V. F. Gerasimov, at 
our request. The same fission chamber and target were 
used, with monochromatic neutrons selected out by a 


‘ mechanical monochromator. Results of the measure- 


ments appear in the diagram. 

The authors consider it a welcome obligation to ex- 
press their gratitude to I. V. Kurchatov for the interest 
kindly manifested in their research [8], to G. N. Yakov- 
lev and S, V. Pirozhkov for preparing a pure thorium 
specimen, and to Chang Huang-ch’ao for his aid in the 
measurements, 
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MEAN NUMBER OF PROMPT NEUTRONS PER SPONTANEOUS 


FISSION OF U238 


/ 
E. K. Gerling and Yu. A. Shukolyukov 


Translated from Atomnaya Energiya, Vol. 8, No. 1, p 49, January, 1960 


Original article submitted June 12, 1959 


To find V3g, the mean number of neutrons per 
single act of spontaneous fission of U238, knowledge of 
two quantities, determined independently of each other, 
is required: the number of prompt neutrons emitted by 
1 g of uranium in unit time An, and the number of spon- 
taneous fissions in 1 g of uranium in unit time ANoge. 
These quantities are related by the simple formula 

Vo39= ee . 
ANos, 
When the number of neutrons accompanying spontaneous 
fission in natural uranium is measured, the most reliable 


result is found to be the one reported in [1]: (1.65 4 0.09)- 


- 10°? neutrons per sec per g of U (cf [2]). 

The generally accepted value for the spontaneous 
fission half-life of U?* is 8.0-10"* years [3]. The corre- 
sponding value of ANog is 24.8 fission events per hr per 
gram of U[1]. Calculation of V2, based on data report- 
ed in [2, 3] yields 2,4 neutrons per spontaneous fission of 
u?*®; this value has been accepted as the most trustwor- 
thy [1,4]. 

By employing the new geochemical method for 
determining the spontaneous fission half-life of us, we 
obtained a value of (5.84 0.5)-1045 years [5] which, upon 
comparison with data reported by other authors in recent 
years [6-9], strongly suggests that the results in [3] were 
too high, and that the true value for the half-life of U238 
is 5.8-10% years, We then have ANpo3g= 34 fission events 
per hr per g of U, while the mean number of neutrons per 
act of spontaneous fission of U238 will be 1.7, instead of 
the previously accepted value of 2.4. 

As demonstrated in an earlier paper [10], the mean 
number of prompt neutrons may be calculated from the 


energy balance of the fission reaction. Results computed 
for PuS®, pua?) pu Pu>, Gm, Cm“; arts and 
cf%2 are found to correspond to empirical data. U?*8 is 
an'exception, calculations giving Jj,,=1.6. When this 
last value is compared to the one we found, there is evi- 
dence for practically complete agreement. 

It is apparent that the value (5.8 +0.5)-104 years 
should be accepted for the half-life of U?** resulting from 
spontaneous fission, while the mean number of neutrons 
per fission event should be assigned the value of 1.7. 


LITERATURE CITED 

1. D. Littler, Proc. Phys, Soc. 65A _, 203 (1952). 

2. K. A. Petrzhak, Physics of the Fission of Atomic 
Nuclei, Supplement No. 1 to Atomnaya Energiya 
(1957). 

3. E, Segre, Phys. Rev. 86, 21 (1952). 

4, J. Fraser, Proc. of the Symposium on the Phys. of 
Fission (Chalk River, Ontario, 1956) p. 238. 

5. E. K. Gerling, Yu. A. Shukolyukov and B. A. Makar- 
ochkin, Radiokhimiya 1, 223 (1959), 

6. P. Kuroda and R. Edwards, J. Chem. Phys. 22, 1940 


(1954). 

7, P. Kuroda and R, Edwards, J. Inorg. and Nucl. Chem. 
3, 345 (1957), 

8. G, Seaborg, The Transuranium Elements (New York, 
1958). 


9. B. D. Kuz'minovy, L. S. Kutsaeva, V. G. Nesterov, 
L. I. Prokhorova, and G. N. Smirenkin, Zhur. Eksp. 
i Teor, Fiz, 37, 406 (1959). 

10, B. D. Kuz'minov and G, N. Smirenkin, Zhur. Eksp. 


i Teor. Fiz. 34, 503 (1958). 


41 


THE EFFECT OF BORON-CONTAINING LAYERS ON THE YIELD 


OF SECONDARY GAMMA RADIATION 


D. L. Broder, A. P. Kondrashov, 
A. A. Kutuzov, and A. I. Lashuk 


Translated from Atomnaya Energiya, Vol, 8, No. 1, pp. 49-51, January, 1960 


Original article submitted August 3, 1959 


In most cases, the weight and size of biological re- 
actor shielding are decidedin terms of hard secondary y 
radiation originating in the reactor vessel and in the first 
shielding layers, as a result of neutron capture. The 
quest for methods of reducing this y-ray flux is therefore 
an extremely timely and urgent problem. 

The most intense secondary y radiation arises in re- 
sponse to radiative capture of thermal neutrons and of 
comparatively low-energy neutrons, formed as a result 
of slowing down in the lighter shielding materials. Ab- 
sorption of these neutronsin a boron-containing layer 
reduces the amount of neutron capture in structural ma- 
terials, particularly in steel which emits gammas at 


high energies [1, 2]. When neutrons are captured in boron, 


soft gammas of low penetrating power and ~0.5 Mev 
energy are generated, and shielding problems are much 
easier to handle. 

The present communication deals with an investiga- 
tion of the effect introduced by a layer of boron carbide 
sandwiched between steel and plexiglas to simulate 
water, upon the rate of production of gamma photons in 
steel. The geometry of the experiment is indicated in 
Fig. 1. Two different grades of steel were investigated: 
St-3 andiKh18N9T stainless. A prism of transverse di- 
mensions 710X710 mm was made from the steel-plexi- 
glas sandwich. 


Detectors 


& Water 

(CD Plexiglas 
@ Steel 

| B§ BgC 


4°20 mm 


Fig. 1. Geometry of the experiment. 
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A Po- &-Be source yielding 2° 10’ neutrons/sec was 
placed in a water shielding arrangement in front of the 
steel prism. Channels with plugs to hold indicators were 
made in the steel—plexiglas prism. Boron-containing 
filters for absorption of thermal and epithermal neutrons 
emerging from the water (or plexiglas) and impinging on 
the steel were inserted between the steel and plexiglas. 
Filters of three distinct types were employed (see vari- 
ants TL Mie Vain rica): 


10° 


10° 


Count, rel. units 


Fig. 2. Spatial distribution of neutrons in the St-3 
plus plexiglas prism: 1) Measurements by indium 
(ByC layer absent); 2) measurements by indium in 
cadmium (B,C layer absent); 3) measurements by in- 
dium (B,C layer 20 mm thick and 1.1 g/cm? density 
sandwiched in between steel and plexiglas); 4) meas- 
urements, indium in cadmium (BgC layer 20 mm thicl 


density 1.1 g/em*, sandwiched in between steel and 
plexiglas). 
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Fig. 3. Spectrum of y quanta formed in the St-3 
prism: 1) No B,C layer present; 2) ByC layer 20 

mm thick, density 1.1 g/em?, sandwiched in be- 
tween steel and plexiglas. 


The spatial distributions of thermal and resonance 
neutrons were measured, as well as the spectra of gammas 
formed in the steel and plexiglas in response to neutron 
capture. The neutron distributions were determined by 
means of round indium foils 20 mm in diameter, con- 
tained within cadmium cases (or without the latter). The 
spectrum of the gamma quanta was measured with a 
single crystal Nal (T1) scintillation y spectrometer 40 
mm in height and 40 mm in diameter. An 11% resolution 
was obtained for the Zn®™ line. A 128-channel pulse- 
height analyzer with ferrite-core memory was used for 
analysis of the pulses. Positioning of the scintillation 
crystal during the measurements is indicated in Fig. 1. 

Figure 2 shows the spatial distribution of neutrons 
measured in the St-3 plus plexiglas prism, both with and 
without a ByC layer sandwiched in. Figures 3 and 4 
give the measured gamma spectra reported herein (here 
also, with a BgC layer present, and in its absence). 

The spectrum of capture y radiation shown in Fig. 
3 contains y lines at energies 2.2 and 7.6 Mev, corre- 
sponding to the absorption of neutrons by hydrogen and 
iron. 

Insertion of a boron carbide layer between the steel 
and plexiglas has the effect of reducing y-radiation in- 
tensity at 7.6 Mev by a factor of 13.4. The spectrum of 
gammas emitted from 1Kh18N9T stainless steel has an- 
other set of y lines, in addition to those mentioned, 
forming as a result this time of capture of neutrons in 
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Fig. 4.Gamma spectrum formed in the prism of 
1Kh18N9T stainless steel: 1) No ByC layer present; 
2, 4) ByC layer sandwiched between steel and plexi- 
glas (variant II in Fig. 1) or layer composed of mix- 
ture of plexiglas, St-3 steel, and ByC (variant IV in 
Fig. 1); 3) BgC layer sandwiched between steel and 
plexiglas (variant III in Fig. 1). 


chromium and nickel, which contribute to the composi- 
tion of this steel. The absorption of thermal and epither- 
mal neutrons by the boron carbide layer, 20 mm thick, 
reduces the gamma flux at 7.6 Mev, emanating from 
1Kh18N9T stainless steel by 7.8-fold. 

The reduction in the yield of capture y radiation 
from St-3 steel when a boron carbide layer is sandwich- 
ed between the steel and plexiglas was computed theo- 
retically. The measured meutron distributions (cf. Fig. 
2) was used in the calculations. The spectrum of neu- 
trons thermalized in steel, used in this calculation, was 
determined in the age approximation for a semiinfi- 
nite steel slug with an absolutely "black" medium at 
the boundary. Corrections were introduced for self- 
shielding on the part of the indicators at the 1.44 ev 
resonance. The computed decrease in intensity of yield 
of secondary gammas in steel was 9.5, which is in ex- 
cellent agreement with the results of the direct measure- 
ments of y spectra. 

The authors would like to express their heartfelt 
acknowledgment to N. A. Alishin, V. S. Borisov, G. V. 
Rykov, and E. V. Shestopalov for their kind assistance 
in the technical work. 
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CRITICAL HEAT FLOWS IN THE FORCED FLOW OF 


LIQUIDS IN CHANNELS 
A. A. Ivashkevich 


Translated from Atomnaya Energiya, Vol, 8, No 
Original article submitted october 13, 1958 


In the design of nuclear reactors in which heat re- 
moval is effected by means of a boiling fluid, the value 
of the critical heat flow q,, at which nucleate boiling 
will be replaced by film boiling must be known. The 
purpose of the present note is to generalize on experi- 
mental data obtained in different investigations, and to 
derive an equation for determining q,.. 

Reference [1] provides a system of criteria for char- 
acterizing the boiling process of a liquid under conditions 
of natural convection. After transformation and discard- 
ing of secondary criteria, this system fakes the form: 

wey” wd, PY: Ul 

eV) oA sat esata, 3 
where w and u are the velocities of the liquid and vapor; 
g is the acceleration due to gravity; v' andr are the 
kinematical viscosity and heat of vaporization of the 
liquid; y’ and y" are specific weights of the liquid and 
vapor; i* andi, are the enthalpy of the liquid at t 
(temperature of saturation) and the enthalpy in the center 
of the heat flow. As the determining criterion, we may 


take the quantity F yee! 
oS ’ nu? 
Vea 


proportional to the collapse diameter of the bubble; o 
is the surface tension of the liquid. The system of crite- 
ria derived in [2] may also be used to derive the system 
in [1]. 

In the boiling of a liquid in channels with forced 
flow, system (1) must satisfy the criteria characterizing 
the channel geometry. By expressing u= Gerstyt [2] and 
substituting the expression for dy into system (1), we get 


q 
Ky =- ” = 7 7 
r* / gy” 9 (y'—¥") 
oa i Pes 
Fe eee een 5 — Sey? 
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where J, is the distance from the beginning of the heat- 
ed interval to the cross section in question; ly is the dis- 
tance from the channel cross section at which x=0 (if x> 
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>0 atthe channelentrance, J, is measured from the be- 
ginning of the heated interval) to the cross section in 
question; x is the vapor content of the flow; dp, is the 
hydraulic diameter of the channel. 

On the basis of treatment of the experimental data 
obtained from references [3-14], an equation was de- 
rived for q,, in the case of forced flow of liquid in chan- 
nels; 


Ker =1,9-10-5Re {1-+-1,8- 10-8 [Ky (4+ Ky)? 


+ Ay] Rey, (3) 
where 
fie aot ean 
Re= 2 Vo Tiss at dnV¥ =? > 1; 
SV AT aay 2V 5 
he wah dhV y' a y" : 
Res Dy" at Ee eee he 
aay eae 
mar(1+4 . ) for tiquid <‘sat ; 
beer eNOS. Ile TAOS OF Aa 
A=1-+0,065 ( > mé=r(i—ax) for «>0; 


Yiquid is the temperature of the center of the liquid 
flow. 

Also confirmed is Eq.(3) by data obtained in experi- 
ments with ethyl alcohol, Dowtherm A, liquid carbon 
dioxide, nitrogen, and hydrogen: qc, values obtained 
from Eq. (3) refer to the cross section for which the de- 
termining criteria are taken. The physical parameters 
of the liquid and vapor used in Eq. (3) are determined at 


Cat: 


Equation (3) is written in the general form both for 
surface boiling, when the temperature of the center of 
the liquid flow tijquid < tsat, and for boiling of the 
vapor- liquid mixture, when the vapor content of the 
flow x >0; for particular cases, it may be simplified. 
For example, for tliquid < tsat, we naturally have K4=0, 
and for x = 0, Ky=0. From the experiments described in 
references [3-14], we find that the stabilization of the 
velocity profile and vapor content sets in at Kg=50 and 
Kg= 125, so that, for Kg > 50 or Ky> 125, the values Kg= 
=50 or Kg=125 must be substituted in Eq. (3). 

Figures 1-3 represent experimental data from sever- 
al investigations [3-14] with a criterial treatment of the 
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data in the form of Ke, (1+Z Re) as a function of the 
function of Re, where Z = 1.8:10 ° [Kg(1+K>) 74+Ky]. For 
convenience, the experimental points are broken up in- 
to groups, with the ordinate scale expanded 5 and 25 
times. Straight lines are drawn on the basis of Eq. (3). 
Points obtained in experiments where water is boiled 
with a nucleus heated to below boiling point are design- 
ated as tliquid < t,,,» and points obtained when boiling 
a vapor-water mixture are designated as x > 0; d is the 
channel diameter; 6 is the slit width; 7 is the heated 
channel length. 

As is clear from Figs. 1-3, the data from different 
investigations, in criterial treatment, show agreement, 
since the initial data in [3] and in [11, 12], [4], and 
[5, 6] coincide; this is explained by the fact that com- 


parison is usually effected while ignoring channel length. 
The spread of experimental points relative to the straight 
lines drawn from Eq. (3) is nonsystematic, and is less 
than # 25-30% for the majority of points. This spread 
may be explained by the limited accuracy of the experi- 
ments and the simplifications assumed: the flow velocity 
w_was determined as the rate of flow divided by the 
channel area and by the specific weight, for average 
temperature or vapor content due to mixing; i, is deter- 
mined as the average value due to mixing. 

Equation (3) has been checked by the experimental 
data reported in 12 investigations carried out with boil- 
ing of water with the center heated to below the boiling 
point, and with a vapor-water mixture in channels of 
round, rectangular and annular cross sections over a broad 
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Fig. 3. Comparison of experimental data with Eq. (3), points 1-8 are taken 
from references [10, 8, 9, 11, 12, 13, 14, 3], resp. 
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range of parameter values: p = 1-220 atm; w= 0.3-24 
m/sec; tear— tiquid = 0- 160°C; x = 0-0.7; dp = 1-30 
mm;/ = 35-1800 mm, i.e., in the following range of 
variation of criteria: Re = 2:10°-4:10°; Ky = 0-0.8; Ky = 
=20-2000; Ky = 1-220; Kg = 0.4-48. 

Equation (3) is recommended for determining cr in 
forced flow of a liquid in channels with the center heat- 
ed to below the boiling point, and in the case of a vapor- 
water mixture over the indicated range of parameter var- 
iation. 
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INVESTIGATION OF HEAT TRANSFER IN THE TURBULENT 


FLOW OF LIQUID METALS IN TUBES 


M. Kh. Ibragimov, V. I. Subbotin, and P. A. Ushakov 


Translated from Atomnaya Fnergiya, Vol. 8, No 
Original article submitted January 6, 1959 


Experiments were carried out on two different pieces 
of apparatus. On one setup mercury was used, on the 
other lead and a eutectic lead—bismuth alloy (43.5% Pb+ 
+ 56.5% Bi). 

Experiments with Mercury. Before being put into the 
apparatus, P-3 type mercury (99.9% pure) was passed 
through a chamois filter. The experimental apparatus 
was a “tube-in-tube" type heat exchanger made from 
1Kh18N9T steel. The mercury flowed inside a tube with 
diameter of 17 X 12.5 mm. The heat given off in the 
pump and loop was carried off from the mercury by wa- 
ter flowing in an annular space 4.5 mm wide. Heat ex- 
change over a length of 760 mm between the mercury 
and water was realized on the counterflow principle, 
with the mercury moving from the bottom upward. Be- 
fore entry into the working part of the tube, allowance 
was made for a hydrodynamic stabilization interval, 
equal to 40 diameters. 

To reduce heat losses over the parts of the tube prior 
to the heat exchanger and after it, compensation heaters 
were mounted on those parts. The temperature of the 
heat-exchange surface between the mercury and water 
at the inlet and outlet were measured by thermocouples. 
The power of the heat exchanger was determined from 
the mass flow rate and cooling of the mercury and was 
controlled according to the mass flow rate and heating 
of the water. The arithmetic mean temperature, calcu- 
lated from the readings of thermocouples placed at the 
inlet and outlet of the heat exchanger, was assumed as 
the mean temperature of the mercury. The mean temp- 
erature of the heat-exchange surface was determined by 
planimetry of the temperature profile over the length of 
the heat exchanger, taking into account corrections for 
depth to which the thermocouples were built in. 


Experiments with Lead and Lead—Bismuth Alloy. 


The exposed liquid metal surfaces were protected against 
oxidation by a cushion of inert gas (argon). Before the 
experiment, the liquid metal was purified of oxides by 
hydrogen regeneration of the latter in the feed tank of 
the apparatus, Chemical analysis showed that the a- 
mount of oxygen in the liquid metal remained practical- 
ly invariant throughout the experiments and was about 
1.107% by weight. 
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. 1, pp. 54-56, January, 1960 


The test part (tube with diameter of 12% 9 mm, 
made from 1Kh18N9T steel) was set up vertically. The 
inner surface of the tube was reamed out. For measuring 
the temperature of the tube wall, a resistance thermo- 
meter made from 0.07 mm platinum wire, 1800 mm 
long and imbedded 118 mm, was used. The distance 
from the point where heating began to the resistance 
thermometer was 33 diameters, 

The temperature of the liquid metal at the inlet to 
the experimental interval was also measured by a plati- 
num resistance thermometer, the heating of the heat 
carrier by chromel-alumel thermocouples. The mass 
flow of the liquid metal was determined from the read- 
ings of a magnetic flowmeter, which was calibrated pre- 
liminarily by the volumetric method. The heat flow 
was created by an electric heater element. The length 
of the heat exchange interval was 520 mm. Heat losses 
were compensated for by an external heater interlocked 
with a thermometer. 

Experimental Results. The mean heat-transfer co- 
efficients for mercury were determined over a heat-ex- 
change interval equal to 60 diameters. Here the pre- 
liminary hydrodynamic stabilization interval was equal 
to 40 diameters. The data obtained in the experiments 
with lead and lead—bismuth alloy refer to an interval of 
the tube removed from the heat source at a distance of 
33 diameters, from the tube inlet at a distance of 41 
diameters. On the basis of [1-3] it may be assumed that 
the results obtained characterize the heat exchange in 
the stabilized zone. The limits of variation of some of 
the quantities during the course of the tests are shown in 
the accompanying table. 

In processing the experimental data,the physical 
parameters of the liquid metals were taken from data in 
[4], the thermal conductivity of the 1Kh18N9T steel 
from data in [5]. The results of the experiments, pro- 
cessed after similarity criteria, are illustrated in the 
figure. In the same graph are plotted curves calculated 
from an empirical formula given in [1}: 


Nu = 4,5+0.014 Pe”’® (1) 
and according to Lyon's semi-empirical formula [6] 


Nu = 740.025 Pe® , (2) 


Variations in the quantities during the tests. 


Parameters 


Temperature, Cs 
Heating, °C 


Liquid metal flow velocity. m/sec 
Heat-transfer rate, kcal/m? «hr . - 
Heat-transfer coefficient, keal/m? -hr:°C 


Reynolds number, Re: 107 
Peclet number, Pe - 10° 
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Results obtained after similarity criteria: 


4 4 


O) Mercury; @) lead; 0) eutec- 


tic lead-bismuth alloy; X) eutectic lead-bismuth alloy, with addition of 


0.1% Mg ; 


The experimental data that we obtained for heavy 
liquid metals are satisfactorily approximated by Eq. (1). 
An addition of 0.1% magnesium to the eutectic lead — 
bismuth alloy led to no significant change in the heat- 
transfer coefficient. At the same time, the experiment- 
al points are distributed far below the values calculated 
according to Eq. (2), with which the data obtained in the 
absence of a thermal contact resistance are found to a- 
SICCH 3, G13 0). 

The experiments whose results appear in the present 
paper were conducted under the usual industrial-experi- 
mental conditions. During the course of the tests the 
liquid metals were not purified of oxides, and the heat~- 


) according to[1] Nu=4.5+ 0.014 Pe®®, —--——) according 
to Lyon's formula Nu=7+ 0.025 pe’ 


exchange surface was not wetted. It may be assumed 
that under these conditions there existed a thermal con~ 
tact resistance, which led to the difference between the 
experimental data and Eq. (2). Our later experiments, 
in which the temperature field in the flow of heavy 
liquid metal was determined, called for this assumption. 
The heat-exchange coefficients found by extrapolation 
of the temperature profile to the wall agreed with the 
values calculated according to Eq. (2). 

Taking part in the experiments and setting up of the 
experimental equipment were Kh. A. Khachaturov, A. A. 
Sholokhov, V. I. Petrovichev, E. V. Nomfilov, and O. V. 
Remizov. 
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DETERMINATION OF MELTING POINTS OF BINARY MIXTURES 
OF URANIUM OXIDES WITH OTHER OXIDES IN AIR 


S. G. Tresvyatskii and V. I. Kushakovskii 


Translated from Atomnaya fnergiya, Vol.8, No, 1, pp. 56-58, January, 1960 


Original article submitted August 3, 1959 


In connection with the increased interest in uranium 
oxides and their reaction with other oxides, the authors 
determined the melting points of a series of mixtures in 
air, 


Starting materials and preparations of mixtures. The 


UO.,, BeO, MgO, CaCO 3, BaCO3, Al,03, La.Os, 
SiQs" 110; 0207). ThOs, GeOs, H,PO,, V.0;, Ta,0;; 
Bi.03, Cr203, (NH4)2MoO,, H,W0,, Fes03, MnCOs. 
Most of the compounds used were classified, according 
to their degree of purity, as pure, analytically pure, and 
chemically pure. The degree of purity of Ta,O;, pre- 
pared by the oxidation of tantalum, and of V2Os5 was un- 
known. The mixtures were prepared by mixing the sam- 
ples thoroughly in a mortar for one hour. The composi- 
tions of the mixtures corresponded to equimolar ratios of 
the oxides used, Additional compositions with molar 
ratios of 1: 2 and 2: 1 were prepared for mixtures of 
uranium dioxide and magnesium, calcium, strontium, 
and barium oxides. To oxidize the uranium di- 
oxide the mixtures were fired in air in porcelain cruci- 
bles at 800°C. Briquettes (25 mm diameter) were molded 
from the oxidized mixtures and fired at 800-900°C. The 
briquettes were ground after firing and the powders ob- 
tained again molded into briquettes, which were fired at 


1000-1100°C. In the next processing cycle the firing 
temperature was raised to 1200-1300°C. At tempera~ 
tures above 1000°C the briquettes were fired in platinum 
crucibles or on a platinum plate. The cycles of grind- 
ing, molding, and firing at 1300°C were repeated until 
the briquettes were strong and showed no signs of crack- 
ing or contracting after firing. The mixtures obtained 
after grinding of the briquettes thus obtained, were used 
for further investigation. 

Determination of melting points of mixtures. The 
melting points of the mixtures were determined by the 
cone deformation method.* The upper section of the 
cone was heated in an electric arc, passing between 
carbon electrodes, The temperature was measured with 
a microoptical pyrometer. The accuracy of the melting 
point determination was +30-50°C for samples melting 
with the formation of a drop, and +50-100°C for samples 
sweating on the surface and forming drops-with difficulty. 
The melting points of the mixtures and the nature of the 
melting are given in the table, 


* With the exception of those mixtures which melted 


during the preliminary firing. (P.P. Budnikov and S.G. 
Tresvyatskii, Coll: Physicochemical Basis of Ceramics) 
[in Russian] (Promstroiizdat, Moscow, 1959) pp. 301-315), 


ol 
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Melting Points of Binary Mixtures of Uranium Oxides with Other Oxides in Air. 


Composition, mol. %_ Melting ' 
In charge U/Me ratio point 
uo, other oxides chemical a °C Netite Cee iiae 
- analysis * 
50 50 BeO Avilere 2200-- 50 Melts with difficulty, almost no drop 
67 33 MgO 1/0,56 1900-+ 50 Drop forms with difficulty 
50 50 MgO EN pt 1750-+ 50 Melts with drop formation 
33 67 MgO 1/t,93 1850- 30 The same 
70 30 CaO 1/0,404 2000-+ 50 a Ke 
DO 50 CaO 1/1 2000-+- 50 iW 
30 70 CaO ROD 2200 £100 No drop forms, sweats on surface 
67 33 SrO 1/0,4 2000+ 90 The same 
50 50 SrO 1/0,93 2100+ 50 te " 
33 67 SrO ice oul 2200+ 50 Melts with drop formation 
67 33 BaO 1/0,54 1700-- 30 The same © 
50 50 BaO 1/1,03 LO ste OO) |e lau 
33 67 BaO 4/4.9 1680+ 00 Drop forms with greater difficulty than in pre- 
vious composition 
50 50 Al,Oy — 1940-+ 30 | Melts with drop formation 
50 50 La,Og 1/1,6 2350+ 10 Drop forms with difficulty 
50 50 SiO, 1/1 1770+ 30 Sweats on surface 
50 50 TiO, _ 14804- 30 | Melts readily 
50 50 ZrOg ~ 2600+ 50 Melts with difficulty, almost no drop 
5 50 ThO, == Does not melt 
50 50 CeO, PURE 2700+ 150 | Melts very little, sweats near electrodes 
n0 50 P2O5 — 12001350 |Sampie melts during preliminary firing 
50 50 V0; — < 800 The same 
50 50 MoOs — 900—1000 | ™ vs 
50 50 WO, = 1100-1200" 1 eae 
50 50 Ta.Os5 — 1850-+- 30 | Melts readily 
AND) 50 Bi,Og — 1800-- 50 Sweats slightly, evaporates considerably on 
heating 
an) 50 PbO, = 1850-- 90 | The same 
D0 50 SnOxg — = Cone does not melt in arc, considerable evapora 
: tion observed 
au 90 Gr,O3 on 2050+ 100 | Drop forms with difficulty 
20 90 Fe,Os eo 1370 30 | Melts readily 
50 D0 MnO, -—- 1650-+- 30 * ss 
100 None == 


*The change in chemical composition was practically insignificant in those cases when n 
carried out, with the exception of mixtures with bismuth, lead, and tin oxides. 


2700—2750-£100] Melts with difficulty 


o analysis was 


THE DISTRIBUTION OF IRON IN MICROVOLUMES OF 


ZIRCONIUM ALLOYS 


P. L. Gruzin, G. G. Ryabova and G. B. Fedorov 


Translated from Atomnaya Energiya,Vol. 8, No. 1, pp. 58-59, January, 1960 


Original article submitted August 5, 1959 


The principal factor limiting the use of zirconium in 
nuclear power reactors is its loss of strength and corrosion 
resistance at high temperatures. It is known that the pre- 
sence of hundredths or thousandths of one percent of ad- 
mixtures of nitrogen, carbon, aluminum, titanium and 
other elements in zirconium [1, 2] markedly reduces its 
resistance to corrosion in hot water or steam. The harm- 
ful effect of these admixtures can be neutralized by the 
presence of tin , iron, nickel, chromium and certain 
other [1, 3] admixtures of alloying elements which have 
a favorable action. With a specific combination of ad- 
mixtures, zirconium iodide smelted in an arc furnace 
can, therefore, possess a sufficient corrosion resistance in 
water and steam at elevated temperatures and pressures. 
The mechanism of the influence of admixtures and alloy- 
ing elements on the corrosion resistance of zirconium has 
been insufficiently investigated; the investigation of the 
distribution of admixtures (for example, carbon [4]) and 
alloying elements (for example, tin [5]) in zirconium by 
the autoradiograph method can be of real value for its 
determination. 

Iron is one of the elements which improve the corro- 
sion resistance of zirconium. Alloys of zirconium and 
tin obtain a maximum corrosion resistance only with 
additions of iron, nickel and chromium or one of these 
elements [3]. In some cases small admixtures of iron 
improve the corrosion resistance of zirconium iodide [2, 
6], particularly in superheated steam [2, 3]. Iron has a 
negligible effect on the strength properties of zirconium. 

The distribution of iron in an alloy of zirconium 
with 0.15 wt.% of iron and in an alloy of the Zircalloy 
type (also with a zirconium iodide basis) of the follow- 
ing composition (in wt.%): tin 1.1, iron 0.1, chromium 
0.1, and nickel 0.05, were investigated by the contact 
autoradiograph method. The radioactive isotope of iron 
Fe™ was added to the alloys in the form of a metallic 
powder when they were smelted in an arc furnace in an 
argon atmosphere. An NIKFI film of the MR type was 
used to obtain the autoradiograms. The Fe™ isotope has a 
y radiation with an energy of 1.295 Mev and B radia- 
tion with a maximum energy of 0.46 Mev. 

To ensure the necessary resolving power, the x-ray 
pictures were obtained from samples with a thickness of 
the order of several dozen microns. The intensity of the 


radiation of the section was 5-10 thousand counts/min- 
cm?. The exposure was selected by experiment. 

The investigation showed that the distribution of 
iron in a cast zirconium alloy and after various types of 
thermal-mechanical treatment remains uneven. In a 
cast alloy (Fig. 1) the major part of the iron is concen- 
trated at the edges of the blocks formed as a result of the 
Ba phase conversion, while part remains in solid solu- 
tion, unable to separate out from the latter because of 
the rapid cooling of the fused system. The distribution 
of iron in cast Zircalloy-2 has a similar character. 

Cold 5-10% deformation gives almost the same 
picture of the distribution of iron as in cast alloys. An- 
nealing of cold-deformed alloys leads to displacement 
of iron from the solid solution to the boundaries of the 
o-phase, Hot forging of alloys in air at 850-750°C re- 
sults in marked disintegration of the structure but the 
unevenness of the distribution of iron is retained. Ther- 
mal treatment of cast and hot-forged alloys was carried 
out in quartz tubes evacuated to ~10°*mm Hg. 

In the case of hardening of alloys from different 
temperatures from the 6 -region the distribution of iron 
in both alloys was similar to the distribution in cast 
alloys. Such distribution is due to the fact that at high 
temperatures (in the 8 -region) iron remains in the solid 
solution [7]. Rapid cooling leads only to partial separa- 
tion of iron at the boundaries formed as a result of B>a 
conversion. Slow cooling of alloys in a furnace from 
temperatures of the B-region leads to more complete 
separation of iron from solid solution at the boundaries 
and sub-boundaries of the @-phase (Fig. 2a and b). 

Figure 3 shows the microstructure and autoradio- 
gram of zirconium-iron cast alloy (not hot-forged) an- 
nealed at 800°C, i.e.,the eutectoid temperature of zir- 
conium-iron alloys [7]. During annealing, the iron 
separated out almost completely from the solution as 
the intermetallide ZrFe,, located at the boundaries and 
within the grains of the a-phase of zirconium. 

Annealing of zirconium-iron alloys in the a-region 
at 600°C (20 hours) and 500°C (40 hours) after chilling 
from 1200 and 900°C ,respectively, showed that the un- 
eveness of the iron distribution was retained. 

Rapid cooling from temperature of the 6 -region re- 
tains part of the iron in the o-solid solution, the other 
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In Zircalloy-2, redistribution of iron takes place 


Slow cooling from high temperatures and annealing in more slowly than in pure zirconium. This is probably 
due to the presence of chromium and nickel in the alloy, 


a-regi 3 ee ca 2 
the &-region assists more complete separation of iron at which — judging from preliminary data — are distribut= 
ed similarly to iron. It can be assumed that the reason 


part is distributed at the sub-boundaries of the a-phase, 


the boundaries of the grains and blocks of the a-phase, 


x70 x 56 


Fig. 1. Microstructure (left) and autoradiogram (right) of cast zirconium alloy 
with 0,15 wt.% Fe. 


x63 


x56 


Fig. 2. Microstructure (left) and autoradiogram (right) of alloys after annealing 
in the 6 -region (1000°C) for 3 hours and slow cooling: a) Zirconium-iron alloy; 
b) alloy of the Zircalloy-2 type. . 
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Fig. 3, Microstructure (left) and autoradiogram (right) of zirconium- 
iron alloy after annealing at 800°C for 3 hours, 


alloys of the Zircalloy group have a higher corrosion re- 

sistance than Zirconium-tin binary alloys is as follows: 

tin, being distributed uniformly in zirconium, increases 
the corrosion resistance of the solid solution, forming the 
main alloy. Iron, chromium and nickel, which are con- 
centrated at the boundaries of the grains and blocks, 
evidently delay boundary corrosion. 
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REACTIONS OF NITROGEN DISSOLVED IN WATER, 
UNDER THE ACTION OF IONIZING RADIATIONS 


M. T. Dmitriev and S. Ya. Pshezhetskii 


Translated from Atomnaya Energiya, Vol. 8, No. 1, pp. 59-62, January, 1960 


Original article submitted January 1, 1959 


The kinetics and mechanism of the reaction of ni- 
trogen with oxygen under the effect of electrons and y 
radiation were investigated in [1, 2]. It was also of in- 
terest to investigate reactions of nitrogen with water and 
its radiolysis products inder the action of radiation, Com- 
parison of these results would make it possible to obtain 
certain data essential for determining the reactivity of ni- 
trogen under irradiation conditions. These processes are 
important for the operation of nuclear reactors in which 
water containing dissolved gases is used as a heat carrier 
or neutron moderator. This paper gives certain charac- 
teristic results, 

Nitrogen dissolved in water is oxidized by ionizing 
radiations to nitrate and nitrite,and forms ammonia. 
Thus, after irradiation chemically active solutions are 
obtained. The formation of nitrate and nitrite was es- 
tablished [3, 4] when y radiation acted on aqueous solu- 
tions of potassium nitrate and nitrite. It was shown [5] 
that at a pressure of 1 atm the radiation reaction is due 
mainly to the absorption of energy in the gaseous phase. 
In the present work Co™ y radiation and fast electrons 
were used. 

From the accelerator the electron beam with a mean 
energy of 0.2 Mev passed through aluminum foil of thick- 


Concentration of NO3 and NOg ions, 
mole Aiter - 104 


1 2 3 
Absorbed energy, ev/em? «10 “19 


Fig. 1, Relation between the concentration of 
nitrate and nitrite ions and the absorbed y- 
radiant energy. Air pressure (in atm): 1) 100; 
2) 20; 3) 1. 
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Concentration of NOz and NO; ions, 


ness 15H and a 1.5 cm layer of air; it then fell on an 
aluminum membrane of thickness 10 closely covering 
the vessel containing the solution, To eliminate reac- 
tions in the gaseous phase above the solution the mem- 
brane was brought into intimate contact with the surface 
of the water. The gas residue was removed by pressure 
from the outer side of the membrane. The absorbed 
energy of the electrons was determined in a calorimeter 
with distilled water by comparison with a dc heater. The 
sensitivity of the calorimeter was 2.10 * deg/min- w. 

Co™ with an activity of 1.4 and 20 kC was used as 
the y-ray source, Irradiation was carried out in sealed 
glass or metal ampoules (of stainless steel). The gases 
were dissolved under pressure in the bidistillate. 

During irradiation the ion concentration in the solu- 
tions was determined by the conductivity. After irradia- 
tion the content of nitrate, nitrite and ammonium ions 
was determined calorimetrically. The nitrate ions were 
determined by reaction with disulfophenolic acid and 
ammonia, the nitrite ions by reaction with o-napthyl- 
amine and sulfanilic acid. The ammonium ions were 
determined by reacting with potassium mercuric iodide 
in alkali, with the formation of ammonium mercuric 
iodide. 


The minimum amounts of acid or ammonia 
Absorbed energy, ev/cem?. 10 “4 


15 


10 


mole Aiter - 1073 


a5 


50 


100 
Duration of irradiation, min 


Fig. 2. Relation between the concentration of 
nitrate and nitrite ions and the duration of ir- 
radiation by electrons with an energy of 200 

kev (gas pressure 1 atm); 1) Nitrogen; 2) air. 


TABLE 1. Relation Between the Reaction Velocity 
and the Duration of Irradiation (air pressure above 
the solution 1 atm, radiation intensity 1.8°10!7ey/em3 


sec) 


Duration of/Concentra- 
irradiation,tion of ni- 


hr 


(), 22 
0,67 
| OS) 
Bee 
2,67 
2,88 


Reaction Velocity 
velocity, .| Constant of 
trate and niqz 107! mole/ the first 
trite, 10°" Jem. sec order, 
mole /em? 104sec™} 
A085 1,20 io 
2,95 0,89 1,48 
6,92 Ofer ou 
6,98 0,15 4,53 
8,60 0,054 (pos 
8,85 0,014 1,56 


TABLE 2. Relation Between the Ratio of the Nitrate 
and Nitrite Concentrations and the Absorbed Energy 


Absorbed Ratio of the nitrate and nitrite 
energy, concentrations 
3 nee 

ev/em Nitrogen, atm Air, 
{ 2 50 1 atm 
1018 By Sh Re 7 as i) 
4()20 Ac Wan) 0,9 0,4 
1022 O25 0,4 0,3 ai 


TABLE 3. Relation Between the Yield of Ammonium 
Ions and the Total Nitrate and Nitrite Ions and the 
Absorbed Energy 


| Absorbed | Relative yield of ammonium ions 
energy, Nitrogen, atm Air, 
ev/cm? { 100 1 atm 
1018 1,6 ile? AW) 

1029 aa 0,9 0,8 

1022 0,9 0,8 Oe 


Concentration of NO3 and NO, ions, 


TABLE 4. Relation Between the Relative Yield of 
Nitrite and Ammonia and the Nitrogen Pressure Above 
the Solution (absorbed energy 2.2:10'* ev/cm) 


| : ! 2 
Pressure, za [NOgI/INOg] INIT; 1/1NO,)+[NOg] 
| ae 1), 85 

10 2,8 ey 
o0) 5a) dee 
100 (ee) 4) 22 


mole/liter - 107° 


Pressure, atm 


Fig. 3, Relation between the concentration of 
nitrate and nitrite formed by the action of y 
radiation and the gas pressure: radiation in- 
tensity 1. 79-10ewWsec> cm? , exposure 50 hr 
(1- nitrogen; 2- air). 


TABLE 5. Energy Yield of Reactions of Dissolved Nitrogen with Water 


| atm 


Pressure above 
the solution, 


Gas composition, % 


Nitrate yield, 
mole/100 ev 


Nitrite yield, 


Ammonia yield) Yield of com- 
mole/100 ev 


mole/100 ev |bined nitrogen, 


100 No 0,016 
80 No+-20 On 0,032 
50 No-1-50 Og 0,028 
80 Ny-+-20 Oo 0,068 

100 Ng 0,138 

100 No 0,238 
80 No +20 Oz 0,296 

100 No 0,267 


atom/100 ev 
0,054 0,108 0,18 
0,021 0,085 0.14 
0,008 0,052 0,09 
0,057 0,205 0,33 
0,216 0,472 0,83 
0,283 0,730 die2a 
OF it 0,615 1,08 
0,342 (0,760 iL aco 
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which could be reliably determined by these methods 
was ~10 *mg. 

The relation (Fig. 1) is shown between the amount of 
nitrate and nitrite formed at various gas pressures and 
the absorbed radiant energy. The amount of reaction 
products formed depends on the pressure above the sol- 
ution. At low degrees of conversion it is proportional to 
the absorbed energy. The amount of combined nitrogen 
is proportional to the absorbed radiant energy up to 10'9- 
10° ev/em®, At large amounts of absorbed energy no 
proportionality is observed. This is illustrated by data 
obtained during irradiation with fast electrons (Fig. 2). 

If gas is passed through the solution during the course of 
the experiment the reaction velocity remains constant, 
but the concentration of combined nitrogen increases 
proportionally to the amount of absorbed energy. The 
energy yield does not vary. The reaction obeys an equa- 
tion of the first order with respect to nitrogen. From 
Table 1 it is seen that the value of the velocity constant 
does not change. 

At a pressure of up to 1 atm the reaction velocity 
is directly proportional to the gas pressure above the 
solution. But at a pressure above 1 atm the proportion- 
ality between the reaction velocity and the pressure is 
disturbed, as may be seen from Fig. 3. 

Table 2 gives the relation between the ratio of the 
concentrations of nitrate and nitrite ions and the absorb- 
ed energy. With an increase in the latter the proportion 
of nitrate increases, 

Table 3 gives data on the relation between the ratio 
of ammonia yield/nitrate + nitrite yield and the ab- 
sorbed energy. With an increase inthe amount of absorbed 
energy the relative yield of ammonia decreases, 

In addition to the amount of absorbed energy, other 
factors influence the yield of various products. 

Table 4 gives data from which it follows that the 
relative yield of the reaction products varies with an 
_ increase in the pressure above the solution. 


58 


Table 5 gives certain values of the initial yields of 
the reaction products during conversion up to 10%. 

Thus, under the effect of ionizing radiations nitro- 
gen dissolved in water is combined in the form of nitrate, 
nitrite and ammonia, The relative yield of each of the 
products depends on the composition of the dissolved 
gas and the amount of absorbed energy. The energy 
yield is 0.09-0.18 atoms of combined nitrogen at a pres- 
sure of 1 atm and 1-1.4 atoms at 150 atm per 100 ev of 
absorbed radiant energy. 

The probable mechanism of the process will be ex- 
amined in the next paper. 
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METHOD OF CALCULATING DOSAGE FIELD 


OF POWERFUL ISOTOPIC UNITS 
N. I. Leshchinskii 


Translated from Atomnaya Energiya, Vol, 8, No. 1, pp. 62-64, January, 1960 


Original article submitted April 4, 1959 


The calculation of the dosage fields of radiators for 
powerful isotopic units is connected with considerable 
mathematical difficulties resulting from the complexity 
of the formulas for the calculations [1-6]. These cal- 
culations can be speeded up considerably and simplified 
by the use of the method presented below. 

To determine dose rate P from a point source in an 
infinite homogeneous medium, 


KM 
Re 


pe (1) 


e-""Bi(hv, pz, Z), 


where M is the activity of the radiation source; R is the 
distance from the radiation source to the point at which 
the dose rate is determined; K is a coefficient taking in- 
to account the dimensions of the quantities and total y 
constant of the radiation source in the formula. Use of 
this formula, however, is difficult, since the absorption 
e"HX and the age factor B depend on the linear absorp- 


tion coefficient of radiation in the substance p, the thick- 


ness of the absorbing layer x, the energy of the radiation 
hv and the effective atomic number of the medium Z. 
Analysis of the product e“/ * B(hv, px, Z) indicates 
that for large thicknesses of the absorbing media it may 
be approximated by an equation of the form f (x)= e"#* 
Brhy, px, Z)= ax?+ bx+c and for small thicknesses, by 


an equation of the form f (x): e HX phy,px, Z)=bx+ 1. 
The results of the approximation and the values of the 
coefficients are shown in the table. 

The approximation permits formula (1) to be repres- 
ented in the form 


P= = (ax?-+-bx-+-c) for large thicknesses (2) 
P= oe (bx +1) for small thicknesses. . (3) 


Formulas (2) and (3) are used for determining the 
dose rate from an elementary source which is part of an 
extended or volume source of radiation. Using these for-, 
mulas for a linear source of activity M and length L at a 
distance R between the point at which the dose rate is to 
be determined and the linear source, and with an absorb- 
ing layer thickness x,, one finds that the dose rate Py is 

® 


il 
7 5M ax? bx, ) 
Ln hl \ aa Gee meg Rare 
0 
¢ 2 (4) 
axy bx, j 
= d 
+\ Costs coo te ) ° | » 
where 
L—y, ; 


Qi =arctg 7 > Pg=arctg 


TABLE, Values of f., and Sizes of the Largest Errors in the Calculation of the Dose Rate 


Thickness of Co” radiation Cs!8 radiation 
Material absorbing | from the formula | greatest from the formula} greatest 
layer f(x) error, % f (x) error, % 
x, cm : 
Ae tee ee Poke ed 0100 1 (b=0) 0,5 1 (b=0) 0,5 
0—10 1 (b=0) 0,0 1 (b=0) 0,0 
Water (biological objects). , D=Ih0) 4,00-10-4,2— 3,0 = — 
; —0,025a-+1 ,0755 
7—70 — — 2+ 10-4a22— 4,4 
—0,0382-+4,1665 
0—10 —0,0135+1 1,8 —0,0135-+-4 fins 
: 0—10 —0,05a-+1 4,7 —0,046x-+1 1,9 
AMON mG Oo OOOO OOO | 05 —070550-+4 ee — 07042214 077 
0—2 —0,056a-+4 0,9 —0,04lr-+1 0,0 
G25 —0,1345a-+ I 3,0 = ae 
Steel oereeoe es ee oe eo © © Oo eo 0—3 —0,1417-+1 0,7 —(),1663a--4 ik 
0—2 —(),1415~-+1 0,4 —(),19452+1 s beef 
Gio MR See SUP AAT CET CieoeR Peis Pies eres? —0,149%-1-4 1,0 —(), 160xe-+1 2,0 
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Fig. 1. Dependence of the coefficient B= arctg 
yi/R + arctg (L-y,)/R on R/L for the following 
values of yy: 1) 1L/2; 2) 3L/4; 3) TL/8; 4) 151/16}. 
5) 4LA; 6) 17L/16; 7) 9LA; 8) 5LA. 


Integrating and inserting the values of ¢, and gp, 
we obtain 


KM at > 
(Py = LR axy R |-ba,A ; cB , (5) 
where 
arctgt 
{= In tg malas 5 tg “ws 
elie 
| arcle R ! )] 
ree 2 ote 
is 


Felli Be Se ; 
KM (ari the +e )( 9h) +-batA-+epB ) ake aan RT ontnA + enB ) 


Fig. 2. Dependence of the coefficient A=In> 


(= pate) =a 
x mw .arcts R sce are R 
ees = al 4 = } 


on RAL for the following values of y;: 1)1L/2; 2)5L/4; 
3) 9L/8; 4) 3L/4 and 1L/2; 5) 7L/8. 


B=arctgt + arctg oats ° 

The quantity B characterizes the dose rate from a 
linear source of radiation without taking into account the 
absorption and the age factor; the quantity A character- 
izes the influence of thin layers of absorbing media; and 
ax’,LA, the influence of thick layers. The values of the 
quantifies B and A are shown in Figs. 1 and 2. 

For the calculation of the absorption of the radia- 
tion in a multi-layer medium, the following expression 
holds; 


to 
att 


where n is the amount of absorbing layers. 

For the calculation of the dose field of a linear 
source along whose length the activity is not uniformly 
distributed, the linear source can be replaced by equi- 
valent point sources, and the dose can be determined 
from the formula 

=n 
P=K 
i= 

Comparison of the results of the calculations with 
experimental data[1, 4, 5] showed that the disparity 


Mj 
pr (aint + biti +e). (7) 


L 


does not exceed 10%. Such accuracy may be considered 
sufficient, since it is within the accuracy of direct meas- 


urements. 
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br-1 ’ (6) 


In designing radiators of any configuration, the ra- 
diators can be represented as consisting of linear or point 
sources, and the dose rate can be determined as the arith- 
metic sum of the dose rates from each of the radiation 
sources. 
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INTEGRATING DETECTOR OF PENETRATING RADIATION 


O. A. Myazdrikov 


Translated from Atomnaya Energiya, Vol. 8, No. 1, pp. 64-65, January, 1960 


Original article submitted October 17, 1959 


The electret state of a dielectric, according to the 
theory of Gross [1, 2], is due to the presence in its vol- 
ume of bound and real charges of density pF (2) = -div Piz) 
and p,,. Surface densities TE + Py and T, correspond 
to these charges. 

Since the methods of measurement of the surface 
charge densities based on the principle of electrostatic 
induction give the total value of this density, we can take 
Tf+Tr= + Pp the normal component of the equivalent 
electric moment of a unit volume of dielectric. 

In considering the influence on the electret state of 
the dielectric (frozen polarization) of one or another 
form of penetrating radiation in order to explain the law 
of time variation of the surface density, we shall start 
from a model of the detector shown in Fig. 1. We assume 
that the electret is placed on a grounded base, and above 
its upper surface, which bears a charge density T, is an 
enclosed volume of a gaseous medium of height Z_ at 
a pressure P_. Of course, a change of T in the irradia- 
tion process will occur, owing to the absorption of energy 
in both the gas volume V¢ and in the dielectric volume 
Vp: The field intensity in the gaseous region, accord- 
ing to the condenser model suggested in [2], is determin- 
ed by the expression 

as: ATT Zr : 
tfatty 
where €,, is the dielectric constant of the dielectric. 
The quantity E reaches a numerical value of hundreds 
of volts and higher and, consequently, the ions produced 
in the gas volume give rise to a saturated ion current. 
The influence of the gas volume is taken into account 
by the density of the leakage current in the gas volume 


JG: 


(1) 


Since the electret charge is bound, the neutral- 
ization of the ions is not possible, and the field of the 
latter decreases the value of T. 

The energy absorption in the dielectric volume leads 
to a more complicated process; the formation of space 


10 


ig 
docamprasers ERTL SS 
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DADA TES EG LA IDET EET 


Figs 1; Model of the detector. 
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charges and hole— electron pairs, i.e., pseudo-dipoles. 
Ionization lies at the basis of these processes, too, and 
hence their effectiveness depends, in particular, on the 
form of radiation. In the end result, these processes de- 
pend on the change of P}. In this connection, the effect 
on the dielectric volume is taken into account by means 
of the leakage current density in the dielectric volume jp - 

Thieisethe radiation intensity, then the general equa- 
tion for the behavior of the frozen polarization in the 
irradiation process can be written in the form 

— Ft migtip =f \ div rav+hs ( div ray. 

‘G Lerner: 

For complete electron equilibrium and a saturated 
condition in the region V, for the y -radiation field, 
formula (2) takes the form 


dr 


PG 4 
—=-_ = eee c Vise a he 3 
c2 Py op ese M2) 


wat 
where c is the electric roentgen equivalent; P, is the 
y “ray dose, Hp is the linear absoprtion coefficient of 
the dielectric (for hu>1 Mev, it is determined mainly 
by Compton scattering); Hof is the pseudo-dipole mo- 
ment; Po is the atmospheric pressure. The first term of 
formula (3) expresses the value of the ion current pass- 
ing through 1 cm? of gas volume of height ZG. The 
meaning of the second term is somewhat more compli- 
cated. As is known, the scattering angles y and @, for 
the case of Compton scattering are related by the expres- 
sion 


(4) 


i.e., in each act of scattering the recoil electron under- 


ctg += (14a) te 0, 


goes a displacement of positive sign in the direction of 
the primary y quanta flux m,. 

Assuming that there is a small probability of re- 
combination of the holes and electrons resulting from 
their capture by trapping centers, we may state that each 
act of scattering leads to the creation of a long-lived 
electron—hole pair (pseudo-dipole), where, according 
to formula (4), these pseudo-dipoles are oriented. If 
we neglect the reverse emission of recoil electrons from 
the trapping centers, if we assume that their concentra- 
tion in these centers is considerably less than the con- 
centration of the centers themselves, and if we take in- 
to account the attenuation law in differential form, then 
the meaning of the second term of formula (3) becomes 
obvious. For this, it is sufficient to sum the number of 
flop-overs in a dielectric column per 1 cm’Aec and 


multiply the result by po¢- The coefficient k takes into 
account the processes of pseudo-dipole multiplication, 
since the recoil electron energy is large. 

The production mechanism of pseudo-dipoles per- 
mits the interpretation of the experimentally observed 
polarization of a dielectric by radiation [3, 4]. A num- 
erical estimate indicates that ig >> j_, even when the 

space-charge field of the dielectric produced by elec- 
trons knocked out of the walls bounding the volume fe 


is taken into account. Therefore, 
dy Jb 5 
—FZ= cZqP, Py ‘ (9) 


If the electret consists of organic dielectrics and con- 
tains hydrogen, then the process of interaction of the 
electret structure with a neutron flux will depend on the 


recoil proton. In this case,formula (2) takes the form 
d 
—=IZP no (6) 
i is the neutron flux density; F is a function 
whose numerical value is determined by the percentage 


content of hydrogen nuclei. 


where P 


Surface charge density T, cm72 


ANT OKe, Ti 


Fig. 2. Dosimetric characteristic for y radia- 
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Integrating expression (5), we have 


“eg Fe o 
e., the change in the surface density Av is a measure 
of the recorded dose D Similar conclusions may be reached 
in the case of neutron fluxes. 

The results of the experimental investigation of 
formula (7) for y radiation is shown in Fig. 2 for pate 
different values of Z 

Thus the ehaheecl in the frozen polarization can serve 
as a measure of the recorded dose, and the electret it- 
self is an integrating radiation detector [5]. The sensi- 
tivity of this detector for y and neutron fluxes, deter- 
mined respectively by the expression 

Sys hg (8) 


and 
Sno=FZp, (9) 


can be adjusted within sufficiently wide limits, owing to 
the terms on the right-hand side of formulas (8) and (9). 
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MEASUREMENT OF Co® y -RAY DOSE CLOSE TO THE 


BOUNDARY BETWEEN TWO BODIES 


V. I. Kukhtevich, B. P. Shemetenko, 


and B. I. Sinitsyn 


Translated from Atomnaya Energiya, Vol, 8, No. 1, pp. 66-68, January, 1960 


Original article submitted August 10, 1959 


In the present work, the passage of Co™ y rays 
through water (Medium I) close to a plane boundary with 
a second medium (Medium II) was investigated, Air, 
lead, nickel, and aluminum were used as Medium II. 
The source-to-detector distance Pp was varied from 0.7 
to 5.0 mean free paths of y rays in water, and the height 
of the source and detector h was from 0,05 to 2,0 mean 
free paths. A diagram of the experiment is shown in 
Fig. 1. The measurements were carried out in a water 
tank of dimension 2.0 X 2.2X 1.6 m. Lead, nickel, and 
aluminum plates of area 90 X 150 cm and thickness equal 
to 2.5 mean free paths of Co”y rays in each of the ma- 
terials of Medium II were placed at the bottom of the 
tank, The measurements at the boundary between water 
and air were carried out in a water tank with an open 
surface. A scintillation y dosimeter with an anthracene 
crystal [1] was used as the y-ray detector. The Co” 
source was made in the form of a sphere of 0.5 cm di- 
ameter; the source activity was 0,153 + 0.005 curie. 

The purpose of the measurements was to determine 
the dose rate D,(p, h) in water close to the boundary 
between the two media. The influence of the boundary 


Ke ee and! U7 


§ 2 3 


on the passage of rays can be characterized by the co- 
efficient L from the ratio Dp(p, he) D,(p, h)=L, where 
Do(p» hea) is the dose rate in a medium of infinite ex~ 
tent measured at the same distance p used for the deter- 
mination of D,(p, h). 

The results of the measurements are shown in Fig. 3, 
where the values of the coefficients L are given for dif- 
ferent distances p and heights h. The errors in the deter- 
mination of the coefficient L did not exceed ¥, The 
obtained results are compared with data [2], in which the 
dissipation energy of y rays of initial energy 1.28 Mev 


Detector 


Medium I 
¢ tf 
SN 

Medium II 


Fig. 1, Diagram of the experiment. 
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Fig. 2. Comparison of the experimental values of the coefficient 1/L (as a function of the distance 
between the source and detector for igh = 0.5 for different materials of Medium Il) with the theoretical 


coefficient K and k’ of [2]. 
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in Medium I (with Z close to that of H,O) at the plane but is several times more dense; then the influence of the 


boundary with a second medium (Medium II) was cal- boundary of separation is taken into account by the cor- 
culated by a Monte Carlo method for two extreme cases; rection coefficient K dependent on source—detector dis- 
1) Medium II has the same atomic number as Medium I, tance and the height of the source and detector above the 
L 
1,40 0 
a a6 
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Fig. 3. Dependence of the experimental values of the coefficient L on p and h for Medium II: 

a) Air (h in cm: 1-2.8; 2-4,8; 3-6.8; 4—8.8); b) lead (h in cm: 1-0.8; 2—2.1; 3-4.4; eA 
5=10.4; 6-16.45 7-20.4); c) nickel (h in cm: 1-0,.8; 2-2.4; 3-4.4; 4-6.4; 5-10.4); d) aluminum 
(h in cm: 1-0,8; 2-2.2; 3-4.2; 46.2; 5- 10.2). 
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level of the boundary of separation; 2) Medium II repre- 
sents a vacuum or material with a very low albedo; in 
this case the correction factor is denoted by K'. Coeffi- 
cients K and K’ are the inverse values of the quantity L 
which we defined earlier, i.e., K=1/L. 

The results of the comparison for pigh= 0.5 are 
shown in Fig. 2. The values of 1/L were obtained by 
interpolation of the experimental data. The values of 
the quantity 1/L for the media water—air and water— 
lead coincide, and within the limits of error of calcula- 
tion and experiment, are in agreement with the values of 
K', The experimental point 1/L for 9p = 1.5 was ob- 
tained by linear extrapolation of this quantity to unity. 
Strictly speaking, such an operation is not entirely just- 
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ified, since*the quantity L tends asymptotically to unity. 
The correction factors 1/L for the media water— nickel 
and water—aluminum lie between the values of K and 
K', as is also the case for the calculations of the author 
of [2]. 

The authors thank S. G. Tsypin for discussion of the 
results of the measurements. 
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ON THE EFFICIENCY OF GAS-DISCHARGE COUNTERS 


V. P. Bovin 


Translated from Atomnaya Energiya,Vol, 8, No. 1, pp. 68-70, January, 1960 


Original article submitted February 13, 1959 


The efficiency * of gas-discharge counters has been 
investigated by a number of authors [1-6]. The formulas 
and experimental curves obtained clearly explain the 
main physical processes occurring during the recording 
of y radiation, but they cannot be used for practical 
purposes, since they do not reflect the peculiarities of 
the design of modern counters and the conditions of their 
use, 

We determined the variation of efficiency with the 
radiation energy in gas-discharge counters of various 
types produced by domestic industry. These counters 
differed from one another regarding the material and 
the design of the cathode and the filling gas. Counters 
of the MS and VS type, filled with an argon-alcohol 
mixture have, respectively, copper and tungsten cathodes 
deposited on the interior of a glass envelope of wall 
thickness 1 mm. Counters of the GS type have a graphite 
cathode and halogen counters of the STS type, a steel 
cathode. 

Figure 1 shows the averaged experimental character- 
istics of five counters of each type. For the individual 
counters of each type, there was observed a spread in 
the efficiency within the limits of 10-20%, but the shape 
of the characteristics remained unchanged. The great- 


est spread in parameters was observed for counters of the 
VS and STS type, and also for all counters filled with an 
argon=alcohol mixture after 10°10" pulses had been re- 
corded, 

Analysis of the curves shown in Fig. 1 indicates that, 
in the quantitative measurements, counters of the MS 
and GS type have better characteristics, but they have a 
high operating potential (1300-1350 volts), Counters of 
the VS type have the highest integral efficiency owing 
to the high coefficient of photoabsorption of y rays in 
the tungsten cathode. The maximum observed on the 
curve is explained by the fact that, for a given region, 
the efficiency is determined primarily by the product Rr, 
where R is the range of photoelectrons in the cathode 
material of the counter andTt is the absorption coeffi- 
cient for y rays. The product Rt has a maximum in the 
energy region close to the K jump of the photoabsorption 
coefficient. 


"By efficiency of y-ray receivers, we willmean the ratio 
of the number of recorded y quanta to the total number 
of y quanta of a given energy entering the radiation re- 
ceiver volume during the same time. 
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Fig. 1. Dependence of the efficiency of gas-discharge counters of 
various types on the radiation energy. 
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The shape of the curves for counters of the STS type 
in the energy region above 1.5 Mev has not yet been 
sufficiently investigated. 

When gas- discharge counters are used, they are, as 
a rule, placed in aluminum or steel sheaths of thickness 
1-3 mm, and sometimes are specially shielded by lead 
filters of different thicknesses for equalizing their spec- 
tral efficiency or are shielded from soft, scattered radia- 
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tion. The use of filters and shields has a basic influence 
on the efficiency of y-ray registration. 

Figure 2 shows the experimentally obtained depend- 
ence of the efficiency on the radiation energy for count- 
ers of type STS-6, STS-2, MS-4, and VS-4 placed in- 
side cylindrical filters of various materials (aluminum, 
iron, lead). In the given case, the efficiency was deter- 
mined to a considerable degree by the interaction be- 
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Fig. 2. Dependence of the efficiency of counters with various filters on the radiation energy: a) STS-6 
counter [(1) without filter; (2) with iron filter of thickness 3 mm around counter; (3) with lead filter of 
thickness 2 mm around counter; (4) with lead filter of thickness 2 mm behind counter; (5) with lead 

filter of thickness 2 mm between counter and radiation source]; b) STS-2 counter [(1) without filter; 

(2) with iron filter of thickness 3 mm around counter; (3) with lead filter of thickness 2 mm around ohne 
er; (4) with aluminum filter of thickness 2 mm around counter};c) MS-4 counter [(1) without filter; (2) 
with iron filter of thickness 3 mm around counter; (3) with lead filter of thickness 2 mm around eCuntery 
d) VS-4 counter [(1) withour filter; (2) with iron filter of thickness 3 mm around counter; (3) with lead 


filter of thickness 2 mm around counter]. 
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tween the wall material of the counter and scattered y 
radiation, Compton electrons, and photoelectrons prod- 
uced in the filter. 

The influence of secondary electrons can be traced 
via the characteristics of the B counter STS-6 with a 
fine steel cathode (Fig. 2a). A lead filter of density 2.3 
g/cm* gave almost no change in the efficiency of reg- 
istration for energies of 0.3-1.5 Mev, since it produces 
a large quantity of secondary electrons which are record- 
ed by the counter. The strong absorption of soft radia- 
tion in lead affects the low-energy region of y quanta, 
and the efficiency sharply decreases. The glass envelope 
used in counters intended for recording y radiation 
shields the cathode material from the action of second- 
ary electrons, and therefore the above phenomena are 
not observed in these counters, 

In iron filters, the predominant process in this en- 
ergy region is Compton scattering. The soft, scattered 
radiation is recorded by counters whose efficiency in- 
creases with the atomic number of the cathode materi- 
al; this is very clearly seen from the characteristics of 
the counter type VS-4 with an iron filter of density 
2.4 g/cm’ (Fig. 2d), This is also confirmed by the fact 
that the efficiency of registration drops sharply after a 


thin filter of thickness 0.1 mm is placed between the 
counter and the iron filter. 

Aluminum filters of small thickness (2-3 mm) in- 
troduce practically no change in the dependence of the 
efficiency on the radiation (Fig. 2b). 

The best filters for gas-discharge counters with glass 
envelopes are lead shields, since they cause the forms 
of the efficiency characteristics of counters of various 
types to approach one another. 

During work with thin-walled counters, it is neces- 
sary to use combinations of filters; iron—lead or alum- 
inum=lead, where the metal with.the smaller atomic 
number should be placed closer to the counter, 

LITERATURE CITED 
1. H. Baeyer, Z. Physik. 95, 417 (1935) 
2. G. Droste, Z. Physik. 100, 529 (1936) 
8. G.Sizoo and D. Coumou, Physik. 3, 921 (1936) 
4, H, Jukawa and S, Sakata, Sci. Papers Inst. Chem. 

Physik. Res. 31, 187 (1937) 

5. E. Norling, Arch. Math, Astr., Och. Fisik. 27a, 27 

(1941) 

6. H. Bradt, P. Gugelot, O. Huber, H. Medikus, P. 

Preiswerk, andP, Scherer, Helv. Physik. Acta. 19, 

77 (1946). 


69 


AIRBORNE RADIOMETER-ANALYZER 
V. V. Matveev and A. D. Sokolov 


Translated from Atomnaya Energiya, Vol. 8, No. 1, pp. 70-72, January, 1960 


Original article submitted December, 23, 1958 


Airborne radiometric searches for deposits of useful 
minerals are being utilized more frequently in geophys- 
cal work, Even with the great achievements in apply- 
ing aeroradiometry to this type of search, the possibility 
of further development and perfection of the present 
methods still exists. As our initial goal, we must refine 
the method of airborne determination of the qualitative 
composition of radioactive ores. 

Though the idea has long been [1] used in other geo- 
physical methods [2, 3], aeronautically this method 
has, up to the present time, hardly been used; this is ex- 
plained basically by the absence of appropriate airborne 
instruments, 

In the present paper we describe a highly sensitive 
scintillation airborne threshold y spectrometer to be 
used in determining the intensity and relative hardness 
of the earth's y radiation. 

The radiometer-analyzer consists of a scintillation 
counter which sends signals to a two-channel threshold 
amplitude pulse analyzer (see block diagram, Fig. 1). 


Fig. 1. Block diagram of the device;.I) counting 
head [(1) scintillation counter; (2) cathode repeat- 
ers]; II) high-voltage supply for the F EU; ITI) elec- 
tronic circuit [(3, 6) cathode repeaters; (4,7, 8) 
amplitude discriminators; (5, 11) measurers of count- 
ing speed; (9) anticoincidence circuit; (10) form- 
ing monovibrator]; IV) automatic indicating and 
voltage supply circuit [(12, 13) recording galvano- 
meters; (14) voltage transformer and rectifiers (15) 
high-voltage battery]. 
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The counting head of the device (the schematic 
cross section shown in Fig, 2) consists of a hermetically 
sealed aluminum cylinder (20 cm in diameter and 60 

cm long) containing a liquid scintillator 2, and two 
photoelectric multipliers 5 of the type FEU-24 or FEU- 
19M. 

A series of experiments using specially constructed 
apparatus [4] was performed in order to prepare the liq- 
uid scintillator which would have the optimum para~ 
meters for a counter with the given geometry (transpar- 
ency to its own radiation, conversion effectiveness, radia- 
tion spectrum, etc.). 

A solution of paraterphenol in toluene was used as 
the scintillator, and this, as is known,is one of the best 
liquid scintillators. Experiments in the choice of opti- 
mum concentration of paraterphenol in toluene were 
based on the fact (this follows from [5]) that conclusions 
obtained from studies of liquid scintillators with small 
volumes may be applied toscintillators with large volumes 
(if in the process there are no basic changes in the spec- 
tral distribution of the scintillator light and the condi- 
tions of light reception at the photocathode of the multi- 
plier are preserved), As a result of the studies we obtain- 
ed a value of 4g/liter as the optimum concentration of 
paraterphenol in toluene; this is in good agreement with 
the data given [6, 7]. 

To the above solution we added napthaphenyloxazole 
(ONPO) (of a concentration of 20 mg/liter), This in- 
creased the sensitivity of the device by approximately 
10%, 

It is well known that scintillators are subject to ox= 
ygen poisoning. It has been shown [8] that this poison- 
ing is reversible and that the properties of the scintillator 
may be restored if inert gases are passed through the solu- 
tion, The mean oxygen extinguishing factor is approxi- 
mately 20%, Actually, after chemically pure argon was 
passed through the solution, the conversion efficiency of 
the scintillator was increased by 10%, In addition, the 
argon atmosphere which was created in the reservoir 1 
(See Fig. 2) served to protect the scintillator from possi- 
ble poisoning during the further use of the device. 

The voltage supply circuit for the photoelectric 
multiplier consists of a high-voltage sectioned battery 
which is located outside of the counting head of the de- 
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Fig. 2. Schematic cross section of the counting head of the device; 1) reservoir; 
2) liquid scintillator; 3) packing, 4) window; 5) FEU; 6) high-voltage divider, 


7) cathode repeater, 


Fig. 3. Exterior view of the electronic circuit and the counting head of the device. 


vice. This method of voltage supply simplifies the de- 
vice significantly and increases the multiplier stability. 

It is necessary to note that the method of feeding 
the photoelectric multiplier from a sectioned battery 
has the disadvantage that the total voltage and the po- 
tential distribution between the electrodes of the photo- 
electric multiplier is fixed. However, it was shown [4, 
9] that for FEU-19M we can choose a mode of operation 
for which the signal-to-noise ratio is high and the opera- 
tion is linear. 

The electronic circuit of the radiometer- analyzer 
is shown in the block diagram of the device, Fig. 1. 

In order to ensure operation within the necessary 
energy limits, switches are provided which permit us to 
change the limits of both channels over a wide energy 
range: from 50 to 500 kev in the channel for the total 
count and from 1.5 to 2.5 Mev in the channel for the 
separate count. 

The limits of the discriminators in the channel for 
the separate count are chosen in such a manner that only 


the intensity of the hard y radiation which has an energy 
exceeding a given energy is registered. Cosmic radia- 
tion with an energy exceeding 3 Mev is excluded with 
the aid of an anticoincidence circuit. Such a decrease 
in the background noise noticeably improves the accu- 
racy obtained in determining the nature of the layers of 
radioactive ores, 

The parameters of the airborne radiometer- analyzer 
are as follows: device sensitivity using the FEU-19M in 
the total count channel is 500+50 impulses/ sec for 1 
mC/hr of radium radiation; the coefficients characteriz- 
ing the relative hardness of the y radiation are equal 
respectively (in percent deviation from the indications 
shown in the total count channel) to 4 for thorium radia- 
tion and to 2 for radium radiation, The exterior view of 
the electronic circuit of the device and counting head 
with the cover of the thermostat raised is shown in 
Figvios 

Trials of the device in 1956 and its use in 1957 have 
shown that its high sensitivity and the possibility of the 
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threshold spectrometry of y radiation distinguish the de- 
vice and permit us to hope that the radiometer-analyzer 
will be useful in aeronautic searches for minerals, 
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INVESTIGATION OF THE PRODUCTION OF AN ELECTROMOTIVE 
FORCE IN A SYSTEM OF SEMICONDUCTORS WITH URANIUM 
DURING IRRADIATION IN A REACTOR 


Yu. K. Gus'kov, A. V. Zvonarev, and V. P. Klychkova 


Translated from Atomnaya Energiya, Vol. 8, No. 1, pp. 72-75, January , 1960 


Original article submitted August 3, 1959 


It is known that under the action of ionizing radia- 
tion on a system of p-type and n-type semiconductors in 
the presence of a contact field an electromotive force is 
produced in the system. Information on this question, 
available in the literature, refers primarily to irradiation 
by light, x rays, y rays, B and a particles[1]. A 
system of p-type and n-type semiconductors irradiated 
by uranium fission fragments was investigated. Some 
compounds of uranium (for example, U3Og), as is known 
[2], are semiconductors. If a system of semiconductors 
with uranium is irradiated by a neutron flux, then the 
fission of uranium leads to the formation of uranium fis- 
sion fragments. These fission fragments, apart from ion- 
ization, produce a large radiation disturbance in the ir- 
radiated materials and thus strongly change the electric 
properties of the latter [3, 4]; therefore, the system of 
semiconductors with uranium should consist of substances 
whose electric properties are not affected very much by 
radiation. Such substances may be semiconductors with 
a large number of initial lattice defects, for example, 
polycrystalline semiconductors or single crystals with a 
large number of impurities. 

In the investigated systems, crystalline oxide semi- 
conductors were used, Uranous-uranic oxide U303, which 
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Fig. 1. Dependence of the voltage 
Vo O) and current I,,(@) of the 
U30,-MgO sample on the neutron 
flux density nv. 


has a large work function, was used as the p-type semi- 
conductor in all experiments, and the oxides BaO, TiO,, 
MgO, and Al,O3, which have small work functions [5], 
were used as the n-type semiconductors. As materials 
for the electrodes we used ‘gold or copper for the p-type 
semiconductors and magnesium or titanium for the n- 
type semiconductors. The samples of U3Os,—BaO and 
U303— TiO, were prepared by successive thermal sputter- 
ing of a system of layers of semiconductors and electrodes 
on a base in vacuo by the method described in[6]. In 
the samples of U3Og—A1,03, the Al,O3 layers were pre- 
pared by metallization on a titanium base; then U3O, 
and a gold or copper electrode were sputtered in vacuo 
on the Al,O3 layer. In the UsOg—MgO samples, it was 
impossible to obtain the MgO layers by sputtering in vacuo, 
and metallization; therefore, we used a ceramic plate of 
MgO of thickness 0.5 mm as the MgO layer, A magnesi- 
um electrode was sputtered on one side of this plate and 
a layer of U3O, and a gold or copper electrode on the 
other, The area of the working surface of the samples 
was 6 and 2.8 cm’, the thickness of the layers of UgOg, 
BaO, TiO, was 5-20 and of Al,O3, 100-200p. 

The prepared samples were fixed in holders with 
clamped bronze electrodes. In order to control the temp- 
erature of the sample surfaces, a copper—constantan 
thermocouple was soldered to one of the electrodes, The 
holders with the samples were placed in aluminum cylindri- 
cal containers and irradiated in the experimental, cooled 
channel of the reactor of the atomic-electric power plant, 
where the density of the neutron flux and y rays was 
10°-10% cm “sec * depending on the power level of the 
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Fig. 2. Load characteristic of the UsOg-MgO sample. 
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reactor, The sample temper’ture during irradiation was 
a}:20 Ce 

During the irradiation, the following characteristics 
were studied for each sample: the no-load voltage Vio, 
the short-circuit current I __, the load characteristic, 
sample surface temperature, and the sample resistance 
R at 1.4 volts in the forward and inverse directions. The 
volt-ampere characteristics for the samples were taken 
before and after irradiation in the reactor. 

During the irradiation in the reactor, an electro- 
motive force was observed to arise in all the samples. With 
an increase in the neutron flux density nv the voltage 
V,. and current Isc increased. The character of the de- 
pendence of the voltage V,, and current Is, on the flux 
density for the UsO,-MgO sample is shown in Fig. 1. It 
is seen from the graph that V,, reaches saturation at 
nv@6X 10" cm™?sec *, and Isc increases linearly with 
an increase in the neutron flux density. The load char- 
acteristic of the U303-MgO sample is shown in Fig 2; 
it is similar to the characteristic of a barrier-layer photo- 
cell. The neutron flux density in this case was 8 X 10” 
cm sec}, 

Figure 3 shows typical volt-ampere characteristics 
of the U30,-Al,03 sample before and after irradiation by 
an integral neutron flux of 5X 10cm™*, From this 
figure it is seen that, first, the sample resistance, after 
irradiation, decreases by an order of magnitude and, 
second, after irradiation, for smal] voltages (0.01-0.1 
volt) a barrier-layer effect, which was absent before ir- 
radiation, was observed, A rather marked rectification 
effect was observed in samples during irradiation; the 
rectification factor here reached 2-10 at a voltage of 
1.4 volts. 

The behavior of the UsO,-MgO samples as a func~- 
tion of the integral irradiation flux is shown in Fig. 4. 
Upon irradiation all characteristics varied greatly, The 
sample resistance before irradiationwas 10 megohm °; after 
irradiation, it fell to 2 kohm, as a result of which I, 
increased sharply; with an increase in the irradiation 
flux, Veo also increased. At an integral irradiation flux 
of ~10'7 cm™2, there was observed a tendency towards 
stabilization of the effect. After the end of the irradia- 
tion of the sample in the reactor, the resistance of the 
sample increased to 1 megohm. 

In order to explain the role of uranium fission 
fragments in the production of the electromotive force 
a special experiment was set up. Several samples were 
given a predose in a y -ray flux before being irradiated 
in the reactor. In this case, the effect in the samples 


was 3 orders less than in the case of irradiation by the 
same neutron flux. This is in agreement with the amount 
of specific ionization produced by y rays and the urani- 
um fission fragments, Double samples of UgO,- Al,O3 
were prepared jn which uranium with 10% enriched U”™® 
was deposited on one half of the sample and natural 
uranium on the other half. During irradiation, the 

10% enriched sample gave an effect 15 times as great 

as the sample with natural uranium. The results of these 
experiments indicate that the electromotive force pro- 
duced in the samples is due to the uranium fission frag- 
ments, 

In addition to the above~mentioned semiconductors, 
oxides and sulfides of beryllium, nickel, molybdenum, 
tungsten, zinc, and copper were examined qualitatively. 
In, all samples, there was observed the production of an 
electromotive force, but the best combination of semi- 
conductors proved to be the combination of samples of 
U303-MgoO. 

Our calculations indicate that in the U30,-MgO 
samples 0.01% of the fission fragment energy is convert- 
ed into electric energy. This low efficiency is apparent- 
ly explained by the small lifetimes of the current cartiers 
and the poor ratio of the diffusion length of the current 
carriers to the thickness of the semiconductor layers in 
the samples. From an analysis of the cotained results it 
is possible to conclude that the production of an electro- 
motive force in the samples depends basically on the 
barrier-layer effect, although volume and thermal elec- 
tromotive force play a certain role. 

In conclusion the authors thank Professor A. K. Kras~ 
in for his interest in this work, G. N. Ushakov for assist- 
ance in performing the experiments, R. G. Bulychev, V. 
A. Shalin, and G, V. Rykov for taking part in the experi- 
ments. 
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INTERNATIONAL SYMPOSIUM ON THE METROLOGY OF 


RADIOACTIVE ISOTOPES 


K. K. Aglintsev and V. V. Bochkarev 


Translated from Atomnaya Energiya, Vol. 8, No. 1, pp. 76-78,January, 1960 


An International Symposium on the Metrology of 
Radioisotopes, organized under the auspices of the Inter- 
national Atomic Energy Agency, was convened in 
Vienna on October 14-16, 1959. About 100 scientists 
from 27 countries, including about 20 representatives of 
various international bodies (CERN, Euratom, Interna- 
tional Standardization Organization, IAEA ) participated. 
The symposium was opened by the assistant general di- 
rector of the IAEA, H. Laboulier, A total of 37 reports 
and papers were read at the various sessions of the sym- 
posium, with lively discussions ensuing. 


Seven nations, Great Britain, Canada, USSR, USA, 
France, West Germany, and Sweden, presented review 
papers summarizing the officially adopted methods for 
the metrology of radioactive isotopes in their countries. 
In these reviews, information was made available on the 
organization of standardization work, and on the methods 
and facilities used in the standardization of isotopes and 
neutron sources, 

Some of the papers contained descriptions of the 
techniques used in fabricating the specimens and meas- 
uring their activity, and some data on the manufacture 
of standard emitters were made available. 

Canada manufactures (at the Ottawa laboratory) up 
to 100 radium preparations annually, contained in plati- 
num-iridium needles, 

The number of standard emitters manufactured in 
France during the years 1957, 1958, and 1959, were, in 
the form of solutions,9, 118, and 129 specimens, respec- 
tively,and 1, 111, and 177 specimens in the form of so- 
lid sources, 


Standard emitters with the long-lived isotopes Co™ 
and Sr” have been manufactured in West Germany, and 
1960 plans call for the manufacture of emitters with the 
short-lived isotopes Na?4, pees and Au, 

Great Britain is producing a large quantity of stan- 
dard emitters; the National Physics Laboratory prepared 
142 standard emitters in 1958, with Na‘, p®, Co, 
Sr? + vo ps1 and Aut (specific activities varying 
from 14 C/g to 1 mC/g); the Harwell Nuclear Research 
Center has put out about 500 sources with 60 different 
isotopes in 1958; the Radiochemical Center at Amer- 
sham produced about 800 emitters during the first nine 
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months of 1959, including 100 a-, 200 y-, and 500 B- 
emitters (of the last, 100 liquid and 400 solid sources). 
The National Bureau of Standards of the USA is making 
available standard emitters with isotopes Hi.Ch Na‘ 
p88 142 C69 785 5p 4 790, p51 187 py t8? 
Au; Hg?®, tu basa as well as standard emitters of Po!” 
U3Og, Ra (D+E) solutions and solid preparations of radium, 

A survey paper on the methods used in the metro- 
logy of radioactive isotopes in the USSR was presented 
by a group of authors (K. K, Aglintsev, V. V. Bochkarev, 
V. M. Grablevskii, F. M. Karavaev). The report made 
a deep impression on the participants at the gathering. 

A report on comparisons made of standard radio - 
active emitters produced in various countries was pre- 
sented by L. Cavallo and W, Mann(USA). The comp - 
arisons, effected during recent years at a number of lab- 
oratories, showed that discrepancies in the determination 
of specific activity were limited: in particular, for the 
standard emitters (in solution form) of Na’ and pe the 
error was 3%; for Na” and Au!® sources, 5%; for p®? and 
K® sources, 1%. This attests to the reliability of the 
methods used in standardizing radioactive preparations. 

There were 21 papers devoted to original research 
methods and measurement techniques for measuring the 
activity of w-, B-, and y -emitters and emitters decay- 
ing by electron capture. H, Robinson from the University 
of California (USA) described a novel design for a pre- 
cision alpha counter based on 1 # geometry, which was 
achieved by means of a special collimating attachment. 
A report by J, Stein and F, Strelov from the National 
Physics Laboratory (Union of South Africa) was devoted 
to the absolute count of q particles in a liquid scintilla- 
tion 4m a-counter, which was successfully used to deter- 
mine the half-life of U?**, 

B, Pate (USA) gave a review of the development of 
the method used in measuring activity with 4m coun- 
ters, and pointed out the advantages of such counters 
Over counters with small solid angles, emphasizing the 
significant influence of 6 radiation on the results of 
self-attenuation measurements, especially for emitters 
with low 6 -particle energies, and the importance, stem- 
ming therefrom of careful manufacture and quality con- 
trol of specimen thickness and regularity for measure- 
ments of specimens on thin organic film. 


The problem of the study of self-attenuation of 8 
radiation in the test specimen was dealt with in two re- 
ports delivered by Canadian scientists. L. Jaffe and J. 
Fishman presented ample experimental data on the de- 
pendence of self-attenuation on the thickness of the 
source for different energies in the spectrum ranging 
from 0,018 Mev (T) to 1.7 Mey (P*?), Self-attenuation 
curves obtained on counters with narrow solid angles and 
on 4m counters were analyzed, and an estimate was 
made of the effect of the geometrical shape factor on 
self-attenuation, P. Campion reported on original re- 
search on self-attenuation correction terms, in particu- 
lar,on new approaches to the extrapolation of self-atten- 
uation curves with the aid of mixedB- and B ~y- 
sources, developed jointly through the efforts of J. 
Taylor and J. Merritt. Calculations of corrections for 
absorption and self-attenuation of B radiation in sources 
with 4m geometry were alsotaken up in a report by A. 
Baptista (Portugal). 

Some of the reports discussed problems involved in 
measuring extremely low-level sources, which are of 
particularly great importance for some problems en- 
countered in biology, medicine, and geology, in re- 
search using labeled atoms, etc. The three basic path- 
ways for increasing the reliability of such measurements 
were considered, viz; devising improved counters, de- 
vising more reliable noise-free electronic networks, and 
reducing background noise,G. Manoff (USA) reported, in 
particular, on work involving a study of the effect of 
cosmic radiation, natural radioactive impurities in the 
structural materials of the counting equipment and hot- 
cell materials, etc, on background; the radiochemical 
purity of reagents was studied, different counting cir- 
cuits (in an anticoincidence circuit for reduction of 
background) were considered, and an estimate was made 
of the relative quality of the circuits for measuring ac- 
tivity levels of the order of 0-10 decay events per min, 
the feasibility of using specially prepared long-lived 
standards, i.e., low-activity standard emitters, for this 
purpose was demonstrated, A number of recommenda- 
tions were advanced for reducing background, choosing 
materials with a low content of radioactive substances, 
etc., by N, Trott (Great Britain), A team of French 
scientists (W. Gallique, B, Grinberg, and M. Tenard) de- 
veloped a low-background device consisting of three 
plane counters which makes it possible to measure the 
specific activity of liquids containing B emitters in 
quantities corresponding to the amounts currently ac- 
cepted as limiting allowable concentrations, This me- 
thod also makes it possible to detect small radioactive 
impurities in solutions, e.g., a 0.027% P** impurity in an 
ov preparation, 

A paper by R, Hasstrom (Finland) on improving li- 
quid scintillators for measuring low activity levels and 
on the use of transistors for fast coincidence counters, 
and reports by M, Forte and A. Ansani (Italy) on a meth- 


od of scintillation counting of very low-energy particles 
(below 1 kev) stimulated much interest, 


J. Stein (Union of South Africa) discussed the re- 
lative efficiency of a liquid scintillation 4 ™ counter 
compared to a proportional 4m™6-counter and a 4m B y- 
coincidence counter, finding that the efficiency drops 
as the energy of B radiation decreases, with a value of 
98.5% reported for Cone 

G, Engelman (France) told of an electrostatic me- 
thod for measuring the activity of B emitters, making 
it possible to standardize p™ and Sr” preparations to 
within 2% accuracy. 


Three papers were devoted to measurements of ra- 
dioisotopes in gas form, K. Hogreve (West Germany) 
discussed the process of absolute measurements of gases 
from the standpoint of the introduction of various cor- 
rective terms (including terms for drift of part of the 
electrons to the counter walls and absorption of the gas 
to be measured on the walls) reaching 2-37 accuracy, 

W. Mann (USA) submitted results of several research ef- 
forts with internally filled compensation counters in use 
at the National Bureau of Standards, particularly on the 
role of inhomogeneity of the electric field, the depen~ 
dence of “wall effect" on the pressure of the gass filling 
the counter, and the effect of pulse discrimination and 
operating voltage on the plateau slope. 


Several techniques for absolute measurements of 
B-active gases with the aid of special devices and in- 
ternally filled compensation counters and some methods 
for absolute and relative measurements of the activity 
of volatile liquids (including liquid organic prepara- 
tions) through their vapors were described in a report de- 
livered by V. V. Bochkarev and V. A. Bazhenov (USSR), 

The British scientists,R. Allen and S, Curran,de- 
livered reports on methods for measuring isotopes decay- 
ing by electron capture, Measurements of low-Z iso- 
topes by means of a 4m(x —y) coincidence counter and 
proportional 4m counter under high pressure (up to 40 
atm), and of high-Z isotopes with high-pressure counter 
for x rays and a scintillation 4m counter were described, 
Several problems concerning the mechanism of the op- 
eration of the counting circuit designed to measure iso- 
topes decaying by electron capture were also discussed, 
particularly the effect of the diameter of the anode lead 
of the proportional counter on its energy resolution. 


Problems concerning calibration of neutron flux 
and standardization of neutron sources were dealt with 
in papers presented by J, Thomas (Denmark) and K. 
Geiger (Canada) calorimetric measurements were dis- 
cussed by M, Lecoin (France) and A. Brinjolfson (Den- 
mark), a Compton dosimeter for measuring x rays and 
high-intensity gammas was discussed by B. Gross (Brazil). 
G. Kainz, K. Scharf, and R. Patselt (Austria) proposed 
a special device to determine fatigue effects in photo- 
multipliers. 
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J. Domanus and B, Osuczewski (Poland) advanced a 
suggestion to redefine the concept of the "gram-equiva- 
lent of radium” as a radiological unit, taking into ac~- 
count the absence of any uniquely determined relation 
between the activity of the source and the dose rate due 
to the source, since y emitters usually are not point 
sources and the dose rate is a factor of self-absorption of 
radiation within the source, 

A report by G. Manoff (USA) on a technique for pro- 
cessing measurement data with the aid of various types 
of electronic computers was met with interest. 


Several other reports bearing on the use of absolute 
counting techniques for solving such problems as varia- 
tion in crystal sensitivity upon boring (K. Bechter, West 
Germany), particle detection in free-neutron decay (T. 


‘Novay, USA), detection and determination of the activi- 


ty of B emitters distributed in specimens, by brems- 
strahlung measurements (J. Patmen, Great Britain) were 
also discussed. 

The proceedings of the symposium will be published 
by the IAEA. 


INTERNATIONAL CONFERENCE ON ACCELERATORS 


A.N. Lebedev 


An International Conference on High-Energy Ac- 
celerators and Scientific Research Equipment, organized 
by the European center for nuclear research (CERN) with 
the support of the International Union of Pure and App- 
lied Physics, met in Geneva, September 14-19, 1959, 
About 350 physicists and engineers from 29 countries 
participated in the conference, 

The basic task confronting the conference was to 
get a picture of the progress made in theoretical and ex- 
perimental research on the development of new acceler- 
ation techniques, leading to greatly increased energies 
and particle beam intensities. The conference also dis- 
cussed ways of perfecting already existing facilities or 
those in the process of construction. 

The conference was opened with a report given by 
a well-known specialist in the field of high-energy phy- 
sics, the leader of the project working on the design of 
the huge 45 Bev linear accelerator, W. Panovsky 
(Stanford U., USA), This report discussed phenomena 
in the field of elementary particles which must be stu- 
died, in Panovsky's view, by means of high-energy ac- 
celerators, as well as the requirements set by the needs 
of the physical experiments onthe energy and intensity 
specifications for the new accelerators, 

Reports by Johns (Univ. of Mich., USA) and O'Neill, 
(Princeton, USA) presented the requirements for asso- 
ciated equipment and particle beam parameters, as well 
as the features of the experimental work to be pursued 
in colliding-beam facilities. These questions are cur- 
rently being much discussed, O'Neill gave a review 
of the work done at Stanford University on storage rings 
of various geometries incorporating injection from lin- 
ear and cyclic accelerators, According to O'Neill, the 
installation of a storage magnet with colliding beams 
in the Stanford linear accelerator will be completed 
sometime in 1960, The system for introducing particles 
with a quick-acting magnetic shutter (operating time 
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lating current values required, 


10%— 10° sec) already installed and tested in the 
Stanford facility was found highly interesting. The 
mechanism for trapping 500-Mev electrons in the stor- 
age system is closely related to the radiation effects ac- 
companying the motion of relativistic electrons, and 
leading to intense damping of transverse oscillations, 
Since it was not possible to simulate these effects at 
low energies, the full-scale facility was built from the 
start. It must be particularly noted that O'Neill ex- 
pects to produce an extremely high vacuum (10% 10 
mm Hg). This vacuum will make it possible to extend 
the beam lifetime in the accelerating chamber to sever- 
al hours, and will consequently provide the high circu- 
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The task of setting up a colliding-beam system has 
another solution which is not related to the need to have 
a separate injector-accelerator, This involves a sym- 
metrical fixed-field ring synchrocyclotron suggested 
about three years ago by both A. A, Kolomenskii (USSR) 
and Okawa (Midwestern Universities Research Associa- 
tion). Facilities of this type make it possible to effect 
simultaneous acceleration and storage of two colliding 
beams of identical particles. As reported at the con- 
ference, this trend is being pursued currently by MURA 


(papers delivered by Cole and Simon), the Physics In- 


stitute of the USSR Academy of Sciences (papers by 
Kolomenskii and Lebedev), and CERN (paper by Penz), 
The fundamental approach common to the three 
teams is, in addition to theoretical research, the de- 
velopment of electron ring synchrocyclotrons of relative- 
ly low energy; e.g., a colliding-beam machine acceler- 
ating to 38 Mev has already been built (but not com- 
missioned) in the USA, and is designed primarily for re- 
search in techniques, CERN is engaged in the design of 
a similar accelerator and is aiming at 100 Mey, the idea 
being to use the machine for physics research. Both 
groups are in addition considering building a stacked- 


beam ring synchrocyclotron as one of the possible ways 
of obtaining a stabilized Budker beam, Reports presen- 
ted by American scientists informed us that work is in 
progress on a 10-15 Bev proton accelerator, which is ap- 
parently to consist of radial sectors with a spiral-modu- 
lated field, 

Soviet papers on the same topic told of theoretical 
and simulated research carried out at the Physics Insti- 
tute of the Academy of Sciences in connection with the 
development of electron ring synchrocyclotrons, 

A large part of the conference reports was devoted 
to acceleration involving the use of plasma, four of the 
eight reports in this area being authored by Soviet 
scientists. Work on so-called plasma betatrons, in 
which electrons stripped from a discharge ignited at the 
beginning of the acceleration cycle in a chamber are 
accelerated, has been moving apace in several coun- 
tries in recent years. CERN, for instance, is building 
five such machines, four of which are already in opera- 
tion (report by Meusonnier), Measurements have shown 
that the value of the current of accelerated electrons in 
all these machines is extremely low, in every case sev- 
eral orders of magnitude lower than what was expected, 
At the present time, research is being pursued to probe 
into the reasons underlying such large losses, 

Another interesting machine of the same type is 
the so-called "megatron,” a plasma betatron designed 
by the Stevens Institute of Technology (USA), which has 
not yet been put into operation, It is designed to ac- 
celerate to 150 Mey, with the current of the order of 
1000 amp. Its most striking feature is the very high 
value of the magnetic field strength (as high as 100 koe) 
and its correspondingly small size. 

Two reports on plasma method were made by Okawa 
(Japan), speaking in the name of two collaborating 
groups. The first report gave an account of the theory 
and development of the plasma betatron, in which the 
magnetic field bears a strong resemblance to that in a 
stellarator, The second report told of a small linear ac- 
celerator designed to accelerate plasma bunches by the 
pressure of radio waves propagated in waveguides, This 
last technique is based on the idea of coherent accelera- 
tion as set forth by V. I. Veksler. 

Several modifications of this basic idea were also 
discussed in papers presented by Soviet scientists Veksler 
and Tsytovich, on some problems dealing with the shock 
mechanism of acceleration in the case of collision of 
rings with relativistic current; Askar'yan, Levin, and 
Rabinovich delivered a paper on the stability of some 
plasma formations, 

As was previously the case, a great deal of interest 
persists abroad in the idea of a stabilized electron beam 
as advanced by Budker. Emphasizing the difficulties 
lying in the way of realizing this idea in practice, 
Finkelstein ( Stevens Institute, USA) told of theoreti- 
cal research performed on the stability of a stabilized 
beam, and in particular,related the view that, from 


this standpoint, the most promising approach is the use 
of an external strong-focusing guide magnetic field. 

Of several new methods proposed for acceleration, 
we might single out the stochastic mode of operation, in 
which interest has grown measurably during the recent 
years, although the idea of the stochatron had been 
voiced by Burshtein, Veksler, and Kolomenskii as long 
ago as 1948, Some of the general problems relating to 
stochastic modes of acceleration were dealt with in the 
report by Kolomenskii and Lebedev, in which a number 
of considerations were voiced on the choice of the spec- 
trum of accelerating voltages, and other topics. A re- 
port on the first working machine embodying these 
principles was made by Keller (CERN), who told of ex- 
periments with a small 4.7 Mev proton stochatron. With 
the aid of simple techniques for creating a “noise” vol- 
tage on the dees, the authors successfully obtained an 
intensity of the same order as in the synchrocyclotron 
mode of operation, In addition, Keller told of the pos- 
sibility of employing the stochastic mode of accelera- 
tion to increase the intensity obtainable on the large 
CERN synchrocyclotron. This same group of topics 
should include the theoretical investigations by Robin- 
son (Cambridge, USA) on an accelerator with simul- 
taneously operating accelerating voltages of different 
frequencies. 

An interesting project for designing an ironless 
strong -focusing synchrotron was related by Christophilos 
(Berkeley, USA). An interesting feature of the magnet 
is that it is to be cooled to a temperature of 80°K or 
even below. The report gave the machine parameters 
and some structural considerations regarding a 24 Bev 
machine of that type. 

Some other reports by Soviet and foreign scientists 
were devoted to cyclotrons with spatial variation of the 
magnetic field. Of those, a report by Zamolodchikov 
and Dmitrievskii (Joint Institute for Nuclear Research, 
Dubna) on a spiral-sector field machine now in opera- 
tion excited particular interest, since it is the first oper- 
ating machine of this type. 

Papers devoted to limitations and perturbing fac- 
tors in accelerators were read at a special session. 

Reports by Simon (MURA, USA), Orlov and Tara- 
sov (Physics Institute of the Armenian Academy of 
Sciences, and Institute of Theoretical and Experimental 
Physics of the USSR Academy of Sciences), Robuche 
and Bariber (CERN) gave further developments in the 
nonlinear theory of betatron oscillations; in particular, 
the last-mentioned paper reported on investigations in- 
to particle losses in passing through a nonlinear reson- 
ance, 

Simon's report further dealt with a theoretical ana- 
lysis of the instabilities arising in a circulating charged 
beam, similar to the discussion in Kolomenskii's and 
Lebedev's report. These effects may be highly impor- 
tant in storage systems, imposing additional limitations 
on the parameters of the system, 
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In a report by Vladimirskii and Kapchinskii (Ins- 
titute of Theoretical and Experimental Physics, USSR 
Academy of Sciences), a detailed exposition of the 
theory of particle motion in a proton linear accelerator 
was given, with the space charge effect taken into ac- 
count. This problem has assumed an important practi- 
cal stature at the present conjuncture, in view of the 
need for higher currents from linear accelerators, par- 
ticularly for injection into large cyclic accelerators. 

Large electron accelerators did not get much at- 
tention at the conference, It could be gathered from 
the report by Robinson (Cambridge, USA) that the con- 
cept formulated on radiation effects as a result of the 
work of Soviet and foreign authors is held to be 
valid at the present time. Most of the papers on elec- 
tron synchrotrons were devoted to theoretical analyses 
of individual effects: choice of damping systems, cal- 
culations of losses with attenuation accounted for, etc., 
(papers by Orlov and Kheifets, Orlov, Kheifets, and 
Tarasov, USSR). 

Of the experimental papers, we note the report by 
Sands Cal Tech, USA) in which the quantitative agree- 
ment between the experimental data obtained from the 
1,2 Bev electron synchrotron and the theoretical calcu- 
lation of particle losses due to quantum fluctuations in 
radiation is pointed out. 

As the session devoted to the state of the art in the 
design of high-energy accelerators, 10 short reports were 
heard on machines being designed or built in various 
countries. Plans for large synchrotrons were presented 
in a report by Livingston on the Cambridge 6 Bev ac- 
celerator, and another by Entschke on the 6 Bev Hamburg: 
accelerator, which is still in the drafting stage. 

Several weak-focusing proton synchrotrons are being 
built in various countries at present, the largest of them 
being the 7 Bev Harwell machine, the 12.5 Bev ANL 
machine, and the 3 Bev machine at Princeton, 

Several papers dealt with the construction of strong- 
focusing proton synchrotrons. Thus, Green (Brookhaven, 


USA) told of a 30 Bev accelerator on which pre-opera- 
tional work has been in progress for one-half year; 
Walton (ORNL, USA) told of a 12 Bev machine, and 
Vladimirskii (Institute of Theoretical and Experimental 
Physics of the USSR Academy of Sciences) gave an ac- 
count of a 7 Bev accelerator. Neal (Stanford, USA) also 
stimulated much interest with his report on the design 
of a linear electron accelerator (10-20 Bev in the first 
cycle, 20-45 Bev in the second). 

No special report was presented on the CERN 25 Bev 
strong-focusing synchrotron, but it became the subject 
of any number of unofficial discussions, in view of the 
fact that a huge success was obtained on that machine 
while the CERN conference was in session, namely the 
achievement of the first acceleration cycle on the 
machine. 

The limited scope of this survey unfortunately does 
not afford us any opportunity to go into the enumerated 
papers in any detail, Furthermore, space considerations 
do not permit mention of many interesting papers, such 
as the contributions of the Khar’kov Physics and Engin- 
eering Institute on linear accelerators and plasma tech- 
niques, as well as quite a few papers on chambers, 
counters, and the use of beams in existing accelerators. 
A detailed account of the proceedings will be contained 
in the full proceedings of the conference, to be pub- 
lished by CERN in the near future. 

The conference demonstrated with full clarity that 
Soviet science has rendered a sizable contribution to 
accelerator physics and engineering; this is attested to 
by the large number of Soviet papers (which comprised 
some one-third of the entire number, and about half of 
the papers on topics concerned solely with accelerators), 
This stands out even more clearly with respect to the 
new ideas propounded by Soviet scientists, which have 
provided the groundwork for further work in many coun- 
tries all over the world, 


AT THE INSTITUTE FOR PHYSICAL METHODS OF SEPARATION 


(GERMAN DEMOCRATIC REPUBLIC) 


N. M. Zhavoronkov and K. I. Sakodynskii 


Responding to the invitation of the German Acade- 
my of Sciences in Berlin, the authors visited Leipzig for 
several days, in which time we became familiar with 
the work accomplished at the Institute for Physical Me- 
thods of Separation (Director: J, Mullenpfordt), which 
is the main center for the dissemination of information in 
work with stable isotopes in the German Democratic Re- 
public (East Germany), 
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The Institute was organized in 1955, Its basic task 
is to promote the use of stable isotopes in scientific re- 
search and to render assistance to all interested organi- 
zations in arranging and performing work using stable 
isotopes, 

The Institute is engaged in intensive work on the 
systematization of basic information on the uses of 
stable isotopes in chemistry, biology, medicine, geo- 


logy, etc., on techniques for handling stable isotopes, 
on the theory underlying separation processes. Several 
review articles written by workers at the Institute have 
appeared in print on these topics, 


Despite the fact that the Institute was organized 
comparatively recently, work of no mean interest has 
been accomplished inthis time. But the chief merit of 
the Institute is that it is the mainspring behind the de- 
velopment of work with stable isotopes in the entire 
German Democratic Republic, making available labeled 
compounds and carrying out isotopic analysis for all in- 
terested organizations, 

The Institute is made up of six sections; a section on 
experimental separation by the fractional distillation 
method (headed by E, Krehl), a section on experimental 
separation by chemical exchange techniques (headed by 
K, Vetzel), a section on the theory of separation methods 
(headed by G, Focht), a theoretical section (headed by 
G. Voit), an analytical section (headed by G. Birkenfeld), 
and a section on the use of stable isotopes (headed by 
_H. Hitbner), 


The tasks of the various separation sections are to 
produce concentrates of stable isotopes of hydrogen, bo- 
ron, carbon, nitrogen, and oxygen. The Institute has 
for this purpose set up facilities for the separation of 
stable isotopes by the methods of fractional distillation 
and chemical isotope exchange. The oxygen 08 iso- 
tope and deuterium are produced by fractionation of 
water in a packed column; the concentration of o'® thus 
obtained is 4%, and deuterium, 7%, To obtain concen- 
trates of the N® isotope, the method of isotopic 
exchange between nitrogen oxides and those of nitric 
acid is resorted to. These techniques yield up to 10% 
enrichment in N®, Experiments have been performed 
on obtaining the carbon isotope oe by fractional distil- 
lation of carbon tetrachloride, It was shown that carbon 
tetrachloride of high purity does not decompose when 
stored in steel equipment, 


Plans for the immediate future call for the expan- 
sion of production of stable isotopes, through the build- 
ing of new and larger facilities, including facilities for 
producing B”’ concentrates by exchange distillation be- 
tween BF; and its anisole complex, 


In the section on experimental separation by frac- 
tional distillation techniques, interesting research has 
also been performed on the effect of surface roughness 
of the packing material on separation efficiency. Op- 
timum values of surface roughness have been ascertained, 
increasing packing efficiency by ~20%, In the same 
section, a proposal has been put forth to use the isobutyl 


ester of benzoic acid as a medium for carrying out the 
droplet technique of isotopic analysis of water. 

The work done in the theoretical sections, whose 
task is to elaborate the theory underlying the separatory 
processes and to calculate the values of the separation 
factors by quantum-statistical methods, is also of con- 
siderable interest. Work has been completed on calcu- 
lation of the separation factors in the chemical exchange 
of several carbon-containing compounds, in the search 
for systems possessing the highest separation factors, 
Some interesting research has also been conducted in 
calculations of the work performed by a nonadiabatic 
fractionating column, A schedule of work on the uses 
of stable isotopes and other methods for separating stable 
isotopes has been drawn up at the Institute, 

Over 30 labeled compounds have been synthesized 
by the section on the uses of stable isotopes, Many of 
the procedures evolved have already been adopted by 
industry, Labeled compounds put on sale include, e.g., 
ammonium nitrate and trimethylamine hydrochloride, 
labeled with N35, acetic acid labeled with of; and 
other compounds, 

The section for the uses of stable isotopes is con- 
ducting independent research on the uses of stable iso- 


‘topes, in addition to providing consultations to various 


organizations, This section is engaged in a study of the 
process of biological transport of water in plants, using 
OV; oxygen exchange between water and glass, ex- 
change between the water of crystallization in sulfates 
and the water of the solution, Interesting results have 
been obtained in a study of the mechanism of oxidation 
of paraffins to fatty acids. Work has also been initiated 
in the field of the uses of isotope methods in geochemis- 
try. 

The analytical section has three mass spectrometers 
at its disposal, One of them, an MS-1, is specially 
equipped for carrying out determinations of deuterium 
content, the second, an MI-1305, is used for determin- 
ing isotope ratios in solid compounds, while the third 
serves for conducting mass isotope analyses, Special 
purpose mass spectrometers are being developed in the 
analytical section, e.g., a mass spectrometer with a 
high-frequency filter, and tests are being performed on 
ion sources of various types. Work on the uses of elec- 
tronic photomultipliers to determine the intensity of 
spectral lines is also a center of interest, 

The analytical section is performing isotope analy- 
sis not only for research being conducted at the Institute 
itself, but is also engaged in analytical research for all 
interested organizations in the German Democratic Re- 
public, 


81 


BUILDING AND DESIGNING OF ATOMIC -POWERED VESSELS 
IN WESTERN AND EASTERN COUNTRIES 


A. V. Klement'ev 


Maritime nuclear power has found its greatest ap- 
plication and development abroad in the area of naval 
and military shipbuilding, where the tactical and tech- 
nical possibilities inherent in atomic- powered vessels 
are paramount and economic performance is a secondary 
factor, 

The use of atomic energy on cargo vessels faces 
great difficulties of an economic and technical na- 
ture. The atomic ships under construction or in the de- 
sign stage at present do not measure up to conventional 
cargo ships as do operating costs, 

According to data published in the foreign press, 
capital expenditures for the building of an atomic ship 
exceed the cost of a conventional vessel by 10-30%, 
Operational expenses due mainly to high fuel costs and 
the presence of special] auxiliary maintenance facilities 
run 15-20% higher. However, the increased cargo ton- 
nage capacity and the unhampered sailing range of 
atomic ships compensates somewhat for the high operat- 
ing expenses, 

It has been established from technical and econo- 
mic surveys that the use of atomic propulsion plants is 
most feasible on large ships for long hauls of liquid 
cargo or bulk cargo. Modern large-tonnage tankers 
spend 80-90% of their operating time on the high seas, 
with the weight of the fuel tank alone accounting for 
10% of the weight on board, in the case of steam tur- 
bine-propelled ships on a long trip. Allowable capital 
outlays on tankers are therefore at a higher level than 
on dry-cargo ships, which are in port about half of their 
operating time. 

All of the atomic power units under construction at 
present for installation in cargo and naval vessels are 
based on water-moderated, water-cooled reactors, 

The most promising power plant approach is the 
single-loop heat removal system in a block with the gas 
turbine. 

A great emphasis in research and experimental 
work is laid on designing a more economical and re- 
liable nuclear reactor, Conventional steam turbines are 
used in the atomic power plants, having proved their 
use fulness over many years of use, 

USA. Paramount attention is given to warships in 
the development of atomic vessels, Work on atomic- 
propulsion passenger and cargo ships is pursued with 
much less intensity, The first vessel of the merchant 
ship variety to be propelled by atomic energy is the 
Savannah (cf, Table), The power plant to be installed 
aboard that ship (similar to the power plant used in the 
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American atomic submarine Nautilus) will not possess 
optimum economic and design parameters. On this 
point, we might mention the long-term program an- 
nounced in the USA in the beginning of 1957 toward the 
development of atomic installations for merchant ships. 
Some machinery manufacturers have been attracted to 
this program. 

In the field of experimental investigation, a good 
deal of attention is being given to atomic gas-turbine 
installations, now being studied by several large com- 
panies, Under AEC auspices, General Dynamics Corp. 
is engaged in the design of a land-based prototype of a 
20, 000 hp maritime atomic power unit, A single-loop 
system with helium coolant and helium as the working 
fluid in the gas turbines is used. The prototype is ex- 
pected to be ready for operation by early 1963, Tests 
of the setup will determine to what extent a single-loop 
system with helium coolant will prove feasible for pro- 
pulsion of sea-going vessels [4], 

The problem of installing an atomic propulsion 
unit on one of the tankers currently under construction 
is being studied, If the decision is positive, an atomic 
tanker of 22, 500 tons deadweight and capable of 20 
knots speed will be launched sometime in 1961 [5]. 

Great Britain. During the past three to four years 
Great Britain has been mani festing no mean interest in 
the use of atomic energy in her merchant and naval fleets, 

The British firms Comel Lord and Babcock and 
Wilcox, Ltd, have drawn up a mutual agreement on re- 
search and design projects aiming at the development 
of an atomic power plant to power a tanker of 65, 000 
tons deadweight .” 

Government scientific research organizations are 
engaged in the development of other projects for ato- 
mic power plants on tankers, in particular 25, 000 and 
50, 000 hp units for tankers of 45,000 and 100, 000 tons 
deadweight, respectively, In contrast to the US, Great 
Britain gives priority to economic feasibility in its 
specifications for atomic propulsion plants, As a re- 
sult of the studies conducted by the British Shipbuilding 
Association in 1957, it was found that the use of ship- 
board atomic power plants will not be profitable for 
some years to come [6], 

The problem of creating a submarine atomic- 


*Deadweight is a term for the total load-carrying capa- 
city of the ship, incl, weight of cargo, passengers with 


baggage, crew with baggage, fuel, lubricants, water, 
etc. 


Essential Data on Ships Under Construction or in Design Stage in Foreign Countries [1-3, 8] 
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powered tanker is getting a good deal of attention, 

For example, on the "Cavess" docks of the Mitchell 
Engineering corporation, tests on various models of un- 
dersea atomic vessels are in progress [7]. 


Japan. In early 1956, a scientific research body 
was set up under the auspices of the Ministry of Trans- 
portation to study atomic shipbuilding possibilities, At 
the present time, this association includes over 50 
scientific organizations and industrial firms throughout 
the country, A five-year plan elaborated by the asso- 
ciation envisages completion of work on the power plant 
for the first atomic tanker sometime in 1961, In 1957, 
the body drew up plans for four types of power reactors 
for tanker propulsion units of 20,000 and 40, 000 hp, 
The schedule of atomic shipbuilding approved by the 
body is based on the use of reactors acquired from the 
USA, in the first stage. Enriched uranium will also be 
imported, 


In addition to this association, several firms are 
involved in atomic ship construction problems, Mit- 
subishi Heavy Industry has worked out plans for two 
tankers with atomic power plants, including one under- 
sea tanker; Hitashi Shipbuilding and Engineering is 
working on a project of an atomic tanker of 65,000 tons 
deadweight capable of making 17 knots, This firm is 
also working on projects for atomic submarine tankers 
and atomic power plants for surface tankers with gas- 
cooled reactors, Plans for an atomic-powered passen- 
ger ship to ply the Japan-South America line have 
reached completion [8]. 


France, In late 1957, the Commissariat de I’ Ener- 
gie Atomique announced a contest for a draft project 


of an atomic power plant of 20, 000 hp to power a tank- 
er of 40, 000 tons deadweight. Several Franch compan- 
ies have developed plans, but practical realization has 
not been achieved to date [9]. 


Sweden. The Hetaverken company has completed 
plans for an atomic tanker of 65, 000 tons deadweight 
(displacement; 84, 000 tons), capable of making 17,75 
knots, A boiling water reactor has been proposed for 
the propulsion unit, The authors of the project feel, 
however, that the vessel will be unable to compete 
with conventional tankers [10].. 

Reports of plans for atomic tankers of smaller ton- 
nages have also come to light in Sweden, in particular 
plans for a tanker of 45,000tons deadweight with an 
atomic power plant of 20, 000 hp. 
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BRIEF COMMUNICATIONS 


USSR. The VII session of the Learned Council of 
the Joint Institute for Nuclear Research at Dubna met 
November 20-25, 1959, at Dubna, Leading special- 
ists in the field of nuclear physics from the 12 member- 
nations of the Institute took part. Reports were deli- 
vered by the Director of the Joint Institute, Correspond- 
ing Member of the Academy of Sciences of the USSR, 
D. I. Blokhintsev, Vice Director Emil Dzhakov (Bul- 
garia), Prof. Wang Hang-ch'ang (China), and the direc- 
tors of five laboratories, giving an account of the re- 
sults of the scientific activity pursued by the staffs un- 
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der their supervision. 

USSR. The atomic icebreaker Lenin has been of- 
ficially incorporated into the Soviet icebreaker fleet by 
the State Commission, Tests completed on December 
5, 1959 showed that all of the mechanisms in the ice- 
breaker are functioning normally, 

The tests were carried out under autumn and winter 
conditions in the Baltic Sea. 

Rumania. Work nas peen completed at the Physics 
Institute of the Rumanian Academy of Sciences on a 
laboratory for work with radioactive isotopes, 
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